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SUMMARY 


This  report  is  piincipnlly  concerned  with  the  use  of  the  triple  conrelatlon  function  and  its 
Fourier  trantfonn,  the  blspectium,  for  the  imaging  of  randomly  translating  objects  or 
objecu  viewed  through  atmospheric  nirbulence.  The  triple  correlation  has  significant 
advanuges  over  the  second  order  correlation,  although  it  does  involve  additional 
computational  complexity.  Specific  conslderadon  is  given  to  algorithms  for  computing 
the  triple  correlation,  the  use  of  a  moving  time  window  and  centroiding  procedures. 

The  Appendix  contains  the  nine  publications  and  conference  papers  that  resulted  tiom 
this  research. 
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1  INTRODUCTION 


This  report  deicribss  the  application  of  triple  correlation  signal  processing  techniques  to 
photon-litnited  data.  Such  dau  will  typically  arise  at  low  light  levels  associated  with 
night  dme  viewing,  but  it  should  be  borne  In  mind  that,  as  detectors  get  faster  in  their 
response,  photon-limited  dau  will  also  arise  at  high  light  levels  (and  very  short 
exposures).  The  type  of  dau  we  are  considering  is  distorted  in  some  way,  either  by 
random  translations  (ax  In  a  tracking  system]  or  by  atmospheric  turbulence  and  thus 
simple  long-exposure  Imaging  is  not  applicable,  as  a  blurred  Image  would  result. 

The  triple  correlation  l(^)(xi,xi)  and  itr  Fcurler  transform,  the  blspectrum  I(5)(ui,u2),  of 
a  one-dimensional  funcdon  l(x)  are  deflnec  by 

i(-)(xi,X2)  ■  j  i*(x)  i(x+xi)  i(x+X2)  dx  (1) 

•  •• 


end 

I<%l.U2)  «  Km)  KU2)  l*(ui+U2),  (2) 


where  I(u)  Is  the  Fourier  transform  of  l(x).  The  triple  correlation  (and  blspectrum)  have 
a  number  of  usefbl  properties; 

•  they  are  Invariant  on  transladon  (this  is  clear  hrom  the  definitions  above) 

•  they  can  be  propagated  through  turbulence  with  llnle  loss  of  information 

•  the  original  signal  can  be  recovered  uniquely  from  either  function  with  no 
ambiguity. 

The  triple  correlation  and  blspectrum  have  recently  been  re  discovered  in  optics  by 
Lehmann's  group  in  Erlangen'0-13  and  this  present  study  wet  motivated  by  their 
results.  This  report  describes  four  aspects  of  the  triple  correlation 

•  application  to  randomly  moving  objects  (}2) 

•  application  to  imaging  through  turbulence  ({3) 

•  use  of  a  moving  time  window  to  improve  signal-to-noise  ratio  (§4) 

•  algorithms  (1.1). 
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In  Kldition,  we  have  Inveidgaiecl  the  um  of  centrolding  for  photon-limited  imeging  and 
demonatreted  a  new  technique  of  Imaging  for  randomly  translating  objects  (#6), 

The  achlevemenu  on  the  above  five  topic',  are  summarised  in  to  J6  on  the  foUowlng 
pages.  References  1-9  are  the  publicadons  resulting  ftom  this  reseamh  and  are 
reproduced  in  the  Appendix,  For  tluther  details  on  each  topic,  the  references  should  be 
consulted. 


2  TRIPLE  CORRELATION  OF  RANDOMLY  TRANSLATING 
OBJECTS 

It  is  clear  from  Bq  (1)  that  the  triple  correlation  of  a  function  l(x)  Is  Invariant  on 
translation  (random  or  otherwise)  and  thus  the  triple  correlation  of  a  randomly 
translating  object  is  simply  equal  to  that  of  the  stationary  object.  Thus,  provided  one 
can  recover  l(x)  from  Its  average  triple  correlation,  we  have  a  simple  way  of  Imaging 
which  is  particularly  relevant  to  optical  systems  involved  in  tracking^'*. 

We  have  applied  tills  to  both  real  and  computer-simulated  data  in  one  and  two 
dimenslons^'^.  Typical  results  for  the  one-dlmenslonal  case  are  shown  in  Pigs  2  At  3  of 
Ref  3;  for  an  object  function  consisting  of  128  sample  points,  unambiguous  imaging  at 
photon  levels  as  low  as  <4  detected  photons  per  ftame  has  been  demonstrated  using  real 
data  gathered  by  a  resistive  anode  imaging  photon  detector  (IPD).  In  this  case  the  dark 
count  was  -0.3  detected  photons  per  frame  and  ••30,(XX)  frame  i  were  lequired  (this 
would  represent  a  data  collection  time  of  "5  mins  at  KX)  frames/s).  Simple-minded 
centroiding  does  not  work  weU  at  this  photon  level  (but  see  (6  below).  The  presence  of 
a  uniform  background  in  the  one-dimensional  esse  wa.«  also  studied;  signal-to- 
background  ratios  >1  did  not  significantly  affect  the  reconstructed  image  quality. 

The  two-dimensional  case  was  studied  in  Ref  4.  For  real  photon  data,  a  simple  object 
consisting  of  >5x3  resoludon  cells  and  -2000  frames,  a  photon  rate  of  -  IS  per  frame 
was  sufficient  to  give  an  adequate  reconstruedon  (Pig  2  of  Ref  4),  although  computer 
simulated  data  of  the  some  object  showed  that  under  these  condldons  the  reconstruction 
should  be  almost  perfect  It  Is  believed  that  imperfecdons  in  the  IPD  were  responsible 
for  the  reduced  quality  in  the  teal  dan  case.  For  computed  simulated  data,  a  more 
complex  object  consisting  of  -8x13  resoludon  cells  and  -2000  frames,  a  photon  rate  of 
-13  per  frame  gave  an  excellent  quality  reconstruedon,  comparable  to,  and  arguably 
better  than,  the  original  stationary  object  imaged  with  -30,000  photons. 

The  conclusion  of  this  study  was  that  the  triple  correladon  is  an  extremely  powerful 
technique  for  high  resoludon  imaging  of  randomly  translating  objects  at  very  low  light 
levels  - 1-20  photons  per  frame.  We  are  considering  extensions  of  the  technique  to 
randomly  rotated  and/or  scaled  objects;  some  progress  on  this  has  been  made  by 
Lohmann’3  in  the  high  light  level  case,  but  the  pracdcal  feasibility  in  the  photon  limited 
case  has  not  yet  been  explored. 
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3  TRIPLE  CORRELATION  IMAGING  THROUGH  TURBULENCE 


Since  (he  invention  of  stellar  speckle  interferometry  by  Labeyrie  in  1970,  a  variety  of 
techniques  for  speckle  Imaging  have  been  suggested^^^'^.  Of  these,  triple  correlation 
imaging  appean  to  provide  the  highest  quality  reconstructions.  Our  research  has 
focussed  on  a  detailed  study  of  the  signal-to-noise  ratio  of  the  photon- limited  triple 
correlation  (or  bispectrum)  and  a  comparison  of  the  Knox-Thotnpson  and  triple 
correlation  methods^^. 

When  imaging  through  turbulence,  the  average  bispectnim  of  the  image,  <I(^Hui,U2)>. 
is  related  to  the  object  bispectrum  0(3)(ui,U2)  by 

<I(3)(ui,U2)>  ■  OW(ui,U2)  <P(^i(ui,U2)>  >  (3) 

where  <P<3)(ui,i<2)>  i*  (be  average  bispectrum  of  the  atmosphere/telescope.  The  phase 
of  <P(3)(ui,u2)>  Is  zero  and  thus  the  phase  of  the  average  image  bispectrum  (which  is 
crucial  for  the  reconstructloi,  process)  simply  equals  that  of  the  object  bispectrum. 

Our  calculadons  and  computer  simulations  show  that,  although  the  whole  atmospheric 
bispectrum  contains  diffhtction-llmlted  Informadon,  the  only  useful  parts  In  the  photon 
limited  case  ate  those  for  which  ui  ■  U2<  In  this  region,  the  signal-to-noise  is  dominated 
by  power-spectrum-like  temts,  yielding  an  SNR  ••  n ,  the  average  number  of  photons 
per  frame.  The  behaviour  at  intermediate  and  high  light  levels,  and  at  other  parts  of  the 
(ui,U2)  space,  is  considerably  more  complex^.  As  a  simple  guide,  fairly  complicated 
objecu  with  aSx  13  resolution  cells  can  be  imaged  with  near-dlffraction-limited 
resoludon  at  a  level  of  ••200  detected  photons  per  frame  and  ••20,000  frames,  in  good 
seeing  (-1  "). 

A  comparison  of  the  Knox-Thompson  and  triple  correiadon  techniques  has  shown  them 
to  be  remarkably  similar,  provided  each  method  is  implemented  in  the  optimum  way. 
The  triple  correiadon  method  is  always  superlot  to  the  Knox-Thompson  method 
because  centroiding  Is  allowed  for  exactly  in  the  triple  correlation,  whereas  It  is  only 
approximately  corrected  for  in  the  K-T  algorithm.  The  triple  correlation  does  not 
Involve  any  more  complicated  or  time-consuming  processing  than  tlje  K-T  method  - 
both  are  equally  computer-intensive. 
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4  USE  OF  A  MOVING  TIME  WINDOW 


Historically,  spatio-temporal  imaging  and  image  processing  has  Mitd  frames  of  data, 
with  the  implicit  assumption  that  the  exposure  time  of  each  fhune  is  sufficiently  shon  to 
"freeze"  any  motion  within  the  scene.  When  dealing  with  photon-limited  data  of  a 
dynamic  nature.,  the  concept  of  data  frames  becomes  less  useful  and  a  more  valuable 
approach  is  to  consider  the  data  as  a  set  of  space-time  (x,y,t)  coordinates  of  photon 
events.  A  time  "frame "  may  then  move  through  this  data  set  conrinuously  and  the 
window  Itself  may  be  of  variable  length. 

As  a  fust  step  in  this  direction,  we  have  evaluated  the  improvement  in  the  slgnal-to- 
noise  ratio  of  an  Nth-order  correlation  when  a  continuously  moving,  but  fixed  width, 
time  window  is  used,  compared  to  the  normal  discrett  frame  window  of  the  same 
length^.  The  detailed  argument  is  given  in  Appendix  A  of  Ref  8.  To  illustrate  the 
method,  consider  the  photons  to  arrive  at  a  uniform  rate  of  N  per  time  window  and 
evaluate  the  number  of  Independent  photon  difference  vectors  possible  in  the  estimation 
of  the  autocorrelation  function  (2nd  order  correlation).  When  the  window  moves 
discretely  to  its  next  position,  the  total  number  of  new  difference  vectors  is  N(N- 1  )/2 
(forward  and  reverse  vectors  are  redundant).  Moving  the  window  continuously  we 
pick  up  (N-1)  new  differences  for  every  new  photon  Included  In  the  window,  thus 
accumulating  a  total  of  N(N- 1)  new  differences  in  all.  In  this  case,  there  Is  thus  an 
Increase  In  the  SNR  of  V2  compared  to  the  discrete  case.  For  the  Nth-order  correlation, 
the  Improvement  Is  -  Nn. 

The  concept  of  treating  the  Image  information  tu  a  set  of  photon  coordinates  is  a  very 
important  one.  We  propose  to  extend  triple  correlation  analysis  to  the  space-time 
domain  using  this  approach.  For  example,  the  length  o)  the  time  window  may  be 
different  for  different  spado)  frequencies,  reflecting  the  knowledge  that,  in  many  cases, 
lilgher  spatial  frequencies  decoirelate  more  quickly  in  time  than  lower  frequencies. 
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5  ALGORITHMS 


Calculation  of  the  complete  tiiple  coneladon  or  bispectnim  is  compuuiionally  intensive 
and  requires  a  large  computer  memory.  For  example,  a  512x512  image  has  a 
bispectrum  consisting  of  (512)^  elements,  each  element  being  the  pnxiuct  of  three 
Fourier  transforms.  For  applicadons  such  as  real-dme  imaging  through  turbulence,  one 
may  wish  to  calculate  and  average  xIOO  bispectra  per  second. 

We  have  Implemented  three  algorithms  for  computing  the  bispecmim  of  two- 
dimensional  photon-limited  dau.  Although  all  three  give  similar  results  and  have 
similar  computational  complexity^,  algorithm  (iii)  (below)  it  preferred,  as  it  enables  the 
moving  dme  window  concept  described  in  g4  to  be  Implemented  efficiently.  In  all  three 
algorithms,  only  a  frecdon  >  t  the  total  bispectrum  is  c<»nputed,  either  that  close  to  the 
axes  as  In  (1)  and  (iii),  or  using  a  limited  number  of  secdons  as  in  (il). 

Algorithm  (i)  involves  direct  Implemenution  of  Eq(2)  and  it  thus  not  llnvted  to  photon- 
limited  data.  A  key  step  in  this  algorithm  (and  (iii))  is  deciding  what  bispectrum 
locadons  to  compute;  this  is  based  on  consideradons  of  signal-to-noise  rado.  Algorithm 
(li)  uses  the  Radon  transform  or  projecdon-slice  theorem  of  Fourier  transforms:  the  two 
dimensional  data  Is  projected  onto  a  series  of  one  dimensional  vectors  06  in  our  case^) 
and  for  each  vector  the  average  triple  coiTeladon  or  bispectium  is  calculated.  The 
corresponding  slice  in  the  object  spectrum  is  found  from  each  bispectrum  and  the  (36) 
slices  are  combined  to  form  the  object  spectrum  and  the  object  map,  as  In  computed 
tomography.  This  approach  may  be  suited  to  highly  parallel  processors.  Algorithm  (iii) 
uses  a  photon  differencing  approach  to  calculate  only  those  pans  of  the  bispectrum  that 
are  necessary  for  the  reconstruedon  and  is  similar  to  the  photon  differencing 
implementadon  of  the  Knox-Thompson  algorithm. 

All  three  algorithms  have  a  similar  compudng  time  on  our  hardware  (SUN  3/160  with 
Sky  'Warrior'  array  processor,  purchased  under  this  Contrret),  -12  s  per  frame  of 
128x128  pixels  at  -200  photons  per  frame  (algorithm  (i)  uses  the  array  processor,  (ii) 
and  (ill)  do  not).  The  algorithms  are  programmed  in  'C  and  could  probably  be 
improved  in  speed  by  a  factor  of  up  to  2. 
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6  CENTROIDING  OF  RANDOMLY  MOVING  OBJECTS 


During  the  study  of  triple  correlation  imaging  of  randomly  moving  objects,  a 
comparison  was  made  with  centroideri  photon  coordinate  data^  The  simplest  type  of 
centroid  is  that  of  the  photon  coordinates  of  each  hame,  without  regard  to  the  number 
of  photons  in  each  frame;  this  Is  not  very  useful,  as  the  results  reported  in  Ref  3  show. 
We  have  studied  another  type  of  centroid,  that  of  the  N-photon  image.  To  compute 
this,  the  data  is  divided  into  sets,  each  set  containing  only  those  images  with  exactly  N 
detected  photons.  A  centroided  image'  is  feund/or  each  set  and,  denoting  its  spectrum 
by  Qn(u).  we  have  shown  that*'** 

Qn(u)  -  I(u(l-l/N])(I*(u;N)lN-t  (4) 

where  I(u)  is  the  spectrum  of  the  original  object.  Note  that  the  spectrum  of  the  N- 
photon  centroided  image'  is  non>linearly  related  to  that  of  the  original  object.  For 
example,  if  N  ■  2,  Qn(u)  i$  to  the  power  spectrum  of  the  object,  and  if  N  ■  3.  it 
is  equal  to  a  particular  section  of  the  blspectrum. 

A  high  quality  reconstruction  of  the  original  eon  be  found  by  combining  all  of  the  Qn'(u) 
using  a  recursive  algorithm*-®,  Using  computer-simulated  one-dlinenslonal  data,  we 
have  reconstructed  a  two-point  object  at  an  average  photon  level  of  -3  photons  per 
frame  and  ■80,000  frames  (see  Fig  3  of  Ref  6).  The  algorithm  has  the  interesting 
property  that  to  reconstruct  the  nth  discrete  spatial  frequency,  one  requires  an  image 
with  at  least  n  detected  photons.  This  method  is  currently  being  studied  for  the 
applicadon  of  imaging  through  tuibulence.  It  has  the  advantage  that  it  is  extremely  fast 
to  compute. 
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7  CONCLUSIONS 


This  research  has  demonstrated  the  potential  usefulness  of  the  triple  correlation  and  its 
Fourier  transform,  the  bispectrum,  in  the  analysis  of  photon-limited  images.  The 
invariance  of  tl.e  triple  correlation  upon  translation  or  propagation  through  turbulence, 
coupled  with  the  fact  that  a  signal  can  be  reconstructed  uniquely  from  its  triple 
correlation,  make  it  a  powerful  tool  for  imaging  in  unfavourable  circumstances. 

We  have  demostrated  the  application  of  the  triple  correlation  to  the  imaging  of  a 
randomly  translating  object  and  objects  viewed  through  atmospheric  turbulence.  In  both 
cases,  dlffracdon-limited  imaging  is  possible  at  very  low  photon  rates,  <10  per  f^e 
for  certain  cases  of  random  transladon  and  "200  per  frame  for  typical  situations  of 
imaging  through  turbulence.  The  signal-to-noise  ratio  of  the  triple  correlation  has  been 
evaluated  in  detail  and  three  algorithms  for  computing  the  bispecirum  have  been  written 
and  compared.  A  photon-differencing  algorithm  is  perferred  as  use  can  be  made  of  a 
moving  time  window,  yielding  a  higher  signal-to-nolse  ratio. 

Centroldlng  Ims  been  studied  as  a  means  of  diffiaction-llmited  imaging  of  randomly 
translating  objects  and  a  new  algorithm  fOT  objea  reconstruction  based  on  centroldlng 
has  been  developed. 
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TuDl-1 


IHAQLNQ  A  RANDOHUt  TRARAATINO  OIJICT  OSINO 
TMt  TRini  CORRB.ATION 

J  C  Dainty  and  H  j  Northoett 
BiaoKatt  Laboratory, 
mpartai  Collaia, 

London  SWT  SRI,  UK. 


Introduatlon 

Tha  trlpla  oorralatlon  function  and  ita  Feuriar  tranafora. 

tha  biapaotruB  of  a  ona  diaanaional  intanaity  diatribution 

I<x)  art  dafinad  by  til.  . 

I^”(Xj,Xj)  •  j  I(x)  Itx+Xj)  I(x+Xj)  dx  (1) 

and 

!^*’(uj.uj)  -  l(Uj>  I(u,)  l(-uj-uj)  .  (3) 

wharo  7(u)  ia  tha  Fouriar  tranafora  of  I(x).  Lehaann  and  Wirnitaar  [11 
hava  ahown  that  it  ia  poaaibla  to  raoovar  a  aignal  I(x)  free  ita  trlpla 
oerralation  or  biapaotrua.  Tha  raoevarad  aifnal  ia  unique  axoapt  in 
oartain  oaaaa  of  no  praotloal  aignifloanoa  in  tha  praaant  atudy. 

It  la  olaar  froa  Iqa.(l)  and  <3>  that  tha  trlpla  oorralatlon  and 
biapaotrua  ara  invariant  with  objaot  dlaplaoaaant  and  Bartalt  and 
Wirnitaar  (11  daaonatratad.  in  a  ooaputor  a.^Bulatlon.  that  a  aharp 
raoonatruotlon  of  a  randoaly  aoving  objaot  on  a  dark  baokground  oould 
ba  obtalnad  at  vary  low  light  lavala  ualng  tha  triple  oorralatlon.  In 
thia  paper  wa  praaant  a  oraotioal  daaonatratlon  of  thia  uaing  phetona 
raoordad  by  an  laaging  Photon  Dataotor  (IPD)  [I]  and  hava  Inoludad  t 
unlfora  background  in  aeaa  of  tha  axparlaanti. 
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TuDl-2 


toaaplMHf 

RaMM  dlaplaoMant  of  a  ona  uiaanaienal  objaot  waa  aohlavad  by  vlavlnc 
tha  ebjaat  paflaotad  ott  a  airror  attaohad  to  a  randoaly  dpivan 
laadapaakar  oeil.  Tha  in  waa  aaakad  by  a  praallt  in  tha  ona 
dlMnalenal  axparlianta  and  photon  avanta  wara  dtcltlaad  into  an  array 
Ij.  J  ■  durlnt  aaoh  tlaa  window  (fraaa)  of  lancth 

apppoxlaataly  (  aa.  Tha  avaraca  dataetad  photon  rataa  In  tha 
axparlaanta  wara  in  tha  ranaa  1-40  par  fraaa  or  llO-dOOO  a~^>  Tha 
lowar  rataa  wara  ooaparabla  to  tha  dark  oount  rata  (-75  t*^)  and  tha 
hichaat  rataa  wara  probably  Influanoad  by  tha  daad-tlaa  affaota  that 
ooour  with  th.a  kind  of  dataotor. 

In  a  typical  runi  a  x  10*  -  5  x  10*  photon  ooordlnataa  wara  raeerdad  in 
10*  -  10*  fraaaa  and  atorad  on  dlao.  In  HP  Pljf  oonputar  waa 
proiraaaad  to  oaloulata  tha  triple  oorralatlon  of  tha  photon  avanta 
ualni  a  aaohlna-oodad  ooordlnata  diffaranolnt  algorltha.  whloh  oould 
oaloulata  tha  trlpia  oorralatlon  at  a  rata  of  100  a**^  for  an  avarata  of 
10  dataotad  photona  par  fraaa. 

A  lone  axpeaura  lu|a  (-5  x  10*  photona)  of  tha  atatlonary  ebjaot  (a 
doubla  allt)  la  ahown  In  Pig.  l<a>.  kith  tho  objaot  aevlng  randoaly  In 
tha  frana.  a  long  axpoaura  ppoduoaa  tha  axtranaly  blurrad  raault  ahown 
In  Pig.  Kb).  Cantroldlng  of  aaoh  fraaa  prior  to  avaraglng 
< 'ahlft-and-add‘ )  glvaa  an  laprovad  raault.  aa  ahown  In  Plga.  Ko),  (d) 
and  (a)  for  avaraga  photon  rataa  of  T.d.  4.4  and  1.1  par  fraaa 
raapaotlvaly.  but  tha  praaanoa  of  tha  oaoond  allt  la  only  apparant  at 
tha  hlghaat  photon  rataa  and  avan  than  la  highly  dlatortad.  Ualng  thi 
trlpia  oorralatlon  nathed  following  tha  raeonatruotlnn  algorltha  of 
Bartalt  at  al  (4)  (which  uaaa  ona  ooaplata  eotant  of  tha  blapaotria). 
tha  raoenatruotlona  ahown  In  Pigs.  Kf).  (g)  and  (h)  ara  obtalnad. 
koto  that  both  silts  ara  olaarly  dlatlngulahabla  in  all  oasas.  In  tha 
oasa  of  Pig.  1(h).  tha  avaraga  dataotad  photon  rata  of  1.1  par  fraaa 
(of  which  0.1  la  dark  oount)  aaant  that  only  7«  of  frtass  oontalnsd  tha 
naoassary  I  or  aora  photons  for  ovaluation  of  tha  trlpia  oorralatlon 
and  tha  raoonstruotlen  was  obtalnad  using  only  20IP  uaabls  fraaoa. 
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TuDl-3 


Pl(ur«  3  fhowt  th«  afftot  of  •  unlfora  baokcrouod  auporlapoMd  upon  tba 
doubla  allt  laafa,  for  baak|reund  to  alinal  ratlea  of  apprexlMtaly  lil 
(ooluan  1)  and  Oil  (ooluan  3)  and  a  al(oal  rata  of  apprexlaataly  4.4 
phetona  a~^,  Flcuraa  3(a)  and  (f)  ahow  tha  raoonatruotlen  ebtalnad 
ualns  a  aingla  llna  of  tha  biapaetrua.  wbaraaa  3(f)  and  <h)  art  found 
froa  a  oeaplata  eotaot.  SaUafaotory  raoonatruatlon  la  ebtalnad  whan 
tha  baotcfreund  rata  doaa  net  aubatantlally  axeaad  tha  alfnal  rata. 


a  b 


g 

jUfPct-i 
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IMAGING  A  FANDOMLY  TRANSUTINC  OBJECT  AT  LOW  LIGHT  LEVELS 
USING  THE  TRIPLE  CORRELATION 

J.C,  DAINTY  ind  M.J.  NORTHCOTT 
Blatktii  tmptnt!  Colltgt,  LmBm  V/C 

Ritiivcd  6  Jinuify  19t6 


An  lttta|in|  photon  doit'tlor  hoi  bMn  UMd  to  rfcofliirt<cl  on  iino|0  ol  o  rondomly  tronilolini  ono'dimoniionol  objtci  oi 
oviroio  photon  rotoo  el  looo  thon  ono  oiinol  photon  por  lioiito.  Thi  mothod  uud  rohoi  on  tho  ihilt-invorionco  ond 
Founor-phoH  prtu rvin|  propottioi  of  tho  Iripio  oorrolotiait  function 


1.  Inlrodticllon 

The  irlpio  coffclillon  function /^V.(|,X})ind  Iti 
Foufler  ttinsform.  the  blipecirum7'’'(U|,U]),of  i 
ond-dimtnilonil  inisniity  dlitrlbutlon  }(x)  irt  di> 
llndd  by  [I), 

and 

7<*’(u,.uj)o7(u,)7(uj)7(-ui  -uj),  (:) 

white  7(ii)  In  thi  Fourlir  trinifoim  of  /(*).  It  hit 
bun  ihown  [1-3)  thst  i  ilfnil/U)  dm  b«  ricovtrid 
from  Its  tripli  corrilitlon  oi  bliptctrum  ind  thil  thi 
ticovirid  ii|ntl  li  unlqut  ixcipt  for  cirltln  cuii  of 
no  prictlcil  liiniflcmci  In  ihi  pritint  iiudy. 

It  li  cliir  from iqi.  (I) ind  (3)  ihil  thi  iilpli  cot- 
filitlon  tnd  blipicirum  in  invirimt  with  objici  dli- 
plicimint.  Btnilt  end  Wlmltzer  [4]  dimonilriiid, 
tn  1  computii  ilmulitlon,  thit  •  ihirp  riconciructlon 
of  1  iindomly  ihlftid  object  on  i  diik  bickgtound 
could  bi  obtained  ei  very  Inw  tight  liveli  ttilng  the  bl' 
ipeptnun.  In  this  paper  we  preient  a  ptacileal  demon- 
•nation  of  thli  uilng  photoni  recorded  by  m  Imaging 
Photon  Diteetor  (IPD)  [3]  md  have  included  o  uni¬ 
form  background  in  eoma  of  the  experlmenti. 


0  030-40ld/86/S03.30  O  Elievler  Science  Publliheri  B.V. 
(Notth-Hollind  Phyilca  Publiihing  Dlvlilon) 


1.  Method 

A  large  number  (typically  3  X  10*)  of  framii  of  a 
randomly  tranilatlng  object  were  recorded  and  their 
average  triple  correlation  computed  uilng  photon  dif¬ 
ferencing  algorithm.  The  average  blipectrum  wai 
found  from  the  avirage  triple  correlation  by  Fourier 
tranaformatlon  and  an  eatimate  of  the  modulus  and 
phaie  of  the  object  spectrum  was  extracted  from  the 
average  blipectrum.  Finally  an  inverie  Fourier  tram- 
form  yielded  m  eitlmata  of  the  original  object. 

Random  diaplacemont  of  the  one  dimensional  ob¬ 
ject  wai  achieved  by  viewing  It  reflected  off  a  minor 
attached  to  a  randomly  driven  louvlipeaKer  coll.  The 
IPD  li  a  microchannel  plate/reilatlve  anode  device 
|5]  with  236  X  236  pixel  coordinates,  which  In  thli 
experiment  wai  masked  by  i  presllt  and  iDltwari  so 
that  the  photon  events  were  digitised  Into  a  onedl- 
menslonal  array  D/,/  •  0,  1,  ...  I3T  during  esch  time- 
window  (frame)  of  approxlmalily  6  ms  duration.  The 
iverige  dirk  count  rite  corre-ponjed  to  0.3  detected 
photons  pet  frame. 

In  1  typical  run,  3  X  10*  frames  are  recorded  and 
stored  on  disc.  An  HP  9836  computer  was  program¬ 
med  to  calculate  the  triple  correlation  ualng  a  machine. 
coded  coordinate  differencing  algorithm.  It  Is  stralght- 
forward  to  show  [6,7]  that  a  histogram  of  vector  dif 
ferinces  of  photon  coordinates  exictly  equals  the  re¬ 
quired  average  triple  corrflatlon  of  the  date. 

The  srerage  blipectrum  of  the  dita  Is  formal- 
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|>  rtlittd  to  thit  of  thf  umplid  Lnt«nilty, 
by  (8), 

♦  (3) 

whttf  ff  l!  tht  avirigi  nutnbir  or  photoni  p«r  rrim«. 

In  fict,  tha  rorractlon  tarmi  In  iq.  (3)  that  hivi  lo 
bt  appllid  to  tha  tripla  eorralallon  to  aiti- 

mata  alntply  cotiaipond  io  Mlf<omtallon  of 
photon  a'^anti.  Thui  an  unbliMd  aitimata  of 
la  obtalnad  If  all  lalf-cotralation  tarmi  ara  Ignoiad  In 
tha  dlraet  calculation  of  thara  li  no  naad  to 

maka  tha  corractlon  via  tha  blipactrum  of  aq.  (3).  Ui- 
ing  our  axlitlni  machlna-codad  loftwara,  tha  unblaiad 
triple  eorralallon  can  be  compulad  at  a  rata  of  appro't- 
Imataly  100 1~‘  for  an  ivaraia  of  10  photoni  per 
frama, 


FIf.  t.  Fhita  riitlavil  rrom  itia  Oupicuum.  Ttii  kanii  tipti- 
lint  viluii  of  iha  bUpaairum,  for  lainipli,  tlii  boa  it  (].]) 
l•ptlHlllla||.  Thf  upper  vilut  or  a  In  nch  boa  li  ihtvilut 
or  thf  object  ipectrum  re  cove  iible  rtom  that  viiue  or  a  ina 
thf  two  fiiUfi  vflufi  of  a  liven  in  the  lower  pert  of  the  boa. 
For  eaemple,  from  In  we  can  Tina  ai  ir  we  alreaay  know  aj 
tnaa|. 

12 


Tho  blipecinim  la  obtained  by  two  dlmenilonal 
Fourier  iranifurmatlon  of  tha  triple  correlation.  Our 
maihod  of  modului  and  phaia  raconitructlon  of  the 
object  ipaciruffi  from  tha  blipactrum  awantlally  fol- 
Iowa  that  of  Bartalt  at  al.  [2]  except  that,  In  daiar- 
mining  tha  phaaa,  a  weighted  average  la  uiad  to  allow 
for  the  varying  modului  of  tha  blipactrum.  If  Of  da- 
notaa  tha  phaaa  of  the  object  ipacirum  a:.d  Si„  da. 
noiai  tha  phaaa  of  the  blipactrum,  than,  for  a  raul  ob- 
Jaci, 

8;m 

If  both  Oq  and  O;  are  let  to  zero,  then  tha  naultlni 

object  la  centred  at  tha  origin.  The  phiios  02,03 . 

ate,  can  be  found  by  lotting  m  ■  I  and  /  •  2, 3, .... 
ate.  Pig.  I  lUuitriiai  how  ivaragad  aitimitei  of  the 
phaiai  Oj  can  be  found  by  averaging  tha  viluoi  obtain, 
ad  eleng  diagonali  of  tha  blipactrum. 


3.  Raaulii 

A  long  axpoiura  Image  (--S  X  10^  photoni)  of  ilj 
itatlonary  object  (a  double  illt)  li  ihown  In  flg.  2(a). 
With  tha  object  moving  randomly,  a  long  axpoiura 
producai  tha  axtramaly  blurted  raiult  ahown  In  fig, 
2(b).  Canttoldlng  of  each  frame  prior  to  avaraglng 
C'ahlft'and-idd'')  glvai  an  Improved  raiult  11  ahown 
In  flgi.  2(c),  ,d)  and  (a)  for  ivariga  dataciad  photon 
ratal  of  7.6, 4.4  and  l.l  per  frame  reipactlvaly  (of 
which  0.5  par  frame  li  dark  count).  Canttoldlng  bare¬ 
ly  ravaali  tha  praienca  of  tha  lacond  ilit,  even  at  the 
highaat  photon  ratal.  Uilng  tha  triple  correlation  and 
blipactrum,  tha  raconitructloni  ahown  In  flgi.  2(fl, 

(g)  and  (h)ara  obtained.  Note  that  both  illti  ara  clear. 
ly  dletbtgulihad  In  all  caiai.  In  fig.  2(h),  the  average 
daiectad  photon  rata  of  l.l  pat  frame,  of  which  0.5  11 
dark  count,  meant  that  only  7'A  of  framai  contained 
iha  nacauiry  3  or  mote  photoni  taqulred  for  the  ava- 
luatlon  of  Iha  triple  correlation  and  the  raconiirucilon 
wae  obtained  uiing  only  2089  framai. 

Fig.  3  ehowe  the  affect  of  a  uniform  background 
lupartmpoiid  upon  tha  double  lilt  Image  for  back- 
ground-to-algnil  tatioi  of  approxunataly  1.4:  I  (col  I), 
2.8  :  1  (col  2)  and  5.4 :  I  (col  3)  and  a  ilgnal  photon 
rata  of  3.1  par  frame.  In  each  caia,  row  3  ihowi  Iha 
raconatruclloni  obtained  uilng  limply  a  ilngla  line  of 
Iha  blipactrum  and  row  4  uiai  tha  complete  octant. 
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Fm.  }.  ESMt  »r  •  wifiKin  tackfiavna  ta  tilpl*  uirtliUBa  iraailnHUm.  (•>-(o>  bap  npavt  taoMU  sf  Wt-  cmn  >bJMi 
alife  u  tvmp  of  } J  plwtooVCnm  hi  Ui<  iliHl  (dPiibb  UH)  ta<  4.),  1.1  tai  K-l  pboua/lbm  li  th>  lattof ■  h(ck|io«ad 
hi  (p).  (b)  Md  («)  iMpaodnlp,  (d)-(n  onholdM  ibstl  fxpoain  laaco,  ((>-0>  iripb  ooaibiba  tmuhiMtleu  ahb|  on*  Ihii 
ol  ihp  M^niua,  <n-(n  iitpb  eomhUap  tteeaitniMbai  utiai  Uw  aanpbia  bbppcuum. 
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(oiM<  or  two^iiMfulonil)  for  10  pheioiu  pit  fltm* 
ind  I  teltilvtly  unlihtforwtrd  hud-wirid  comb* 
lor  biMd  on  CMOS  luflo  could  tehlivt  the  uin«  itu 
tor  SO- 100  photons  p«r  toinM,  thus  tUowlni  on-line 
computitlon  of  tht  triple  eoritlttlon.  The  prinelpel 
problem  In  extendini  this  lo  two  dimensions  lies  In 
the  stortie  eepecliy  required  for  one  complete  ociint 
of  the  bispeetrum  (**32  Mwords  for  e  I2B  X 128 
ims|e).  It  msy  be  possible  to  use  s  selection  of  ptsnes 
In  the  bispeetrum  to  reconstruct  the  object  with  only 
t  imtU  lots  In  il|neMo<noiis  ritlo  compired  to  i  re- 
consiruetlon  using  the  complete  octent. 
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A  ktprtnt  from  iht 


Photon  limitod  imaging  through  turbultnct 
0,  R,  Ayara,  J.  C.  Dainty  and  M.  J.  Northeott 
Optica  Saotlon,  ilaehatt  laboratory,  Imporlal  Collaga,  london  211,  UK 


Tha  Imaging  of  oblaeta  through  atmoapharlo  turbulanea  hat  baan  atudlad  for  two 
particular  eaiati  (II  whan  tha  talaaoopa  dlamatar  D  la  muoh  anallar  than  tha  atmoapharlo 
ceharanea  tcala  r  and  (11)  whan  0  >>  r  .  Caaa  (1)  produoat  randomly  moving  Imagaa  and 
caaa  (11)  raaulta  In  tpaoKlad  Imagaa.  In  both  eaaaa,  triple  oorralatlon  taehnlquaa  aro 
Invaatlgatad  and  aomparad  with  poaalbla  altornaclva  mathedt,  tuoh  aa  eantreldlng  tor  (1) 
and  tha  Knox-Thompton  mathod  for  (11).  Tho  triple  oorralatlon  taohnlqua  appaara  to  bo  the 
moat  aatlafaetory  rna  for  both  problamt  at  low  light  lavala. 

Introduotlon 

Tho  Imago  of  a  email  object  viewed  through  atmoapharle  turbulanoa  ahowa,  In  general, 
two  eharaetorlatle  faaturoai  It  randomly  tranalatoa  and  la  brolton  up  Into  a  apoelda 
pattern.  A  datallad  doacrlptlon  of  tha  affoet  of  turbulonua  In  Imaging  la  given  In 
Rat.  1,  but  tha  two  eharactarlatloa  can  be  undorttood  in  almpla  tarma  by  conaidarlng  tha 
eaaaa  where  the  talaaoopa  dlamatar  0  la  vary  much  amallor  than  tha  atmoapharle  eoharanea 
aoala  r  or  vary  much  greater  than  It.  in  tho  former  eaaa  (D  <<  r  )  a  point  oouroa  object 
appaara  aa  a  randomly  tilted  wavefront  In  tha  taleaoopa  pupil  and  thua  glvaa  rlaa  to  a 
randomly  tranalating  image  of  tha  objaot.  In  loot,  randomly  tranelatlng  Imagaa  oan  alao 
arlaa  If  tha  object  Itaalf  la  randomly  tranalating  or  If  tha  talaaeope  la  doing  ao,  for 
example,  baeauaa  of  guiding  arrora  or  wind  buffeting  tha  aupport  atruotura.  Tha  method 
of  objaot  reeonatruotlon  daaerlbad  In  thla  paper  doea  net  depend  on  tha  cauaa  of  tha 
random  ahift  In  tho  Image. 

In  the  aacond  oaaa  (0  >>  r.).  It  la  wall>linown  that  a  apeeklaJ  Image  la  produeed  by  a 
point  aourea  object  aa  a  raauit  of  Inter feronea~et  randomly  phaaed  aiaaa  In  the  tojaaeopa 
pupil.  Tha  technique  of  tpeoltla  Intbrforometry  '  rooovara  a  dlffraetlen  limited  auto- 
oorralatlon  of  tha  object  Intenalty  and  It  worka  at  very  low  light  levela,  at  lew  at  a 
few  detected  photont  per  frame  for  timple  atrugturea  aueh  aa  binary  atart.  A.multltude  of 
apaeltle  imaging  taehnlquaa  have  been  auggeeted*.  The  Knex>Thompaon  teohnlque’’’  oalculatea 
phata  dlffarenooa  In  tha  Fourier  plane,  the  phate  itaalf  being, found,  In  effeot,  by 
integration  of  tha  dltferoneaa. .  Tho  triple  oorralatlon  method  which  la  related  to 
tha  technique  of  phaaa  cloaure  >  la  apparently  more, difficult  to  implement  than  tha  Knox* 
Thompaon  taohnlqua  but  poaaeaaea  Inherent  advantagaa  ’  dleeueaad  below. 

In  thla  paper  wo  praaont  axparlmantal  reaulta  and  computer  almulatlcna  for  tha  two- 
dlmanalonal  Imaging  of  a  randomly  tranalating  objaot  and  oomputor  almulatlon  raaulta  for 
apaclda  imaging,  In  bot)-'  oaaaa  at  vary  low  light  lavala.  Ttao  datalla  of  tha  computational 
tachnlquoa  and  detail  d  comparlaona  of  dlffara.it  methoda'’  are  praaantad  olaawhara,  Wa 
atart  with  a  *  laf  aummary  of  tha  triple  correlation  and  Ita  relation  to  othar  correlation 
functlona. 


Tha  triple  correlation 


Tna  triple  oorralatlon  l*^’(fl|,Kjl  of  a  real  function  l(x)  la  daflnad  by° 


i(3)( 


Sl'iSj*  •  /  *  Hi  •  ♦  Hj* 


(1 ) 


and  Ita  four  dlmunalonal  Fourlar  tranaform,  tha  blapactrum  I’^Nu^.Uj)  by 
l’”(u,,u-)  ■  K^,)  I*(Ui  +  Ug)  • 


(?) 
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where  Ilu)  le  the  two  dimeneional  Pourlar  trenefortn  of  Kxl.  Note  theti  for  two 
dimenelonel  functlone  the  triple  correletion  end  blipeetr^in  ere  foor  dlmenelonel  end 

thle  ceueet  en  l|iiiedlete  problem  regerdinp  oemputer  etorege.  tven  taking  account  of  the 
known  aynmetriee*,  one  would  require  500  Hworda  of  memory  for  the  biipectrum  of  e 
25(  X  35(  function. 

Pot  compariaon.jthe  autocorrelation  function  i*^Mx.)  and  Ita  Fourier  treneform,  the 
energy  epectrum  l'  are  defined  ter  real  i(x)  by'* 


I  •  J  1  If)  ilu  *  i<i  > 


(3t 


and 


i*^’(u, ) 


1(U,  I  1  (Ji,  ) 


(41 


( KT  ( 

Ttaa.lfnox-Thompeon  correletion  i  (x,, 
defined  by 


Au)  and  Ite  Fourier  tranaform  (with  reepcet  to  x, ) 


i**’'>(x,,Au) 


i ( X  *  X, )  axp( IxAu.x)  ^S. 


(5) 


and 


,Aji) 


Kji, )  I*(l,  ♦  du) 


(t) 


where  Ai^  ia  a  fixed,  amal'  vector. 

If  Cq.(3)  le  rewritten  with  Uj  ■  Au,  then 

l’*’(u,,Au)  -  KAu)  I  (2,,  Au)  . 


(7) 


Thue,  a  tingle  plena  of  the  biepeetrum  for  fi  email  fixed  value  of  Au,  ^e  equal  to  the 
Knox-Thempaon  croae>epwetrum  multiplied  by  an  ob)act-dapandent  quenTity  l(Au). 

Thera  are  aevaral  different  tachnlquea  that  can  bn  uced^to  compute  the  varioui  correlation 
functione  and  their  Fourier  traniforme  at  low  light  levele  For  direct  eelculation 
correlation,  a  vector  differencing  algorithm  ie  the  moat  efficient  at  lew  light  ievela 
However,  roconatruotion  algorithme  that  require  the  biepeetrum  would  neeeaaitata  taking  the 
four-dimanaional  Fourier  traneform  of  thia  triple  correlation  and  thia  it  net  practical 
on  grounde  of  epaed  or  computer  memory  for  large  arraye.  In  auch  caeoa  we  either  compute 
eelected  aub-plenae  of  the  biepeetrum  or  ute  tomographic  tachnlquea.  Selected  planet  of 
the  biapeotrum,  Bq.(7l,  are  found  by  celoulatlDC  their  inverte  Fourier  traneform  directly 
using  a  weighted  photon  differencing  algorithm  ”  and  then  Fourier  tranafermlng  the  renult, 
or  by  direct  Fourier  tranaforming  using  Eq.(3l.  in  the  tomographic  method,  the  two 
dimenelonel  photon  data  ia  projected  onto  a  line  at  angle  S,  thue  reducing  tna 
dimensionality  of  the  calculation.  Thle  ie  repeated  for  tevnrtl  angles  (16  in  tnr  results 
below).  In  all  c..ett  of  photon  differencing,  bias  terms  due  to  aelf-oorrelition  of  photons 
are  avoided  by  excluding  from  tna  calculatlont  all  vector  dlffttencee  between  a  photon 
event  and  ittelf. 

The  method  of  raconetructirn  of  the  object  epectrum  from  the  biepeetrum  sseantia'.ly 
followed  that  daecribed  in  Ref*  9  and  10,  with  eomj  speed  imprevement.  The  object  intanaity 
wee  found  by  aimpla  invaraa  Fourier  traneformatlun.  both  of  theta  raconetructlon  atopa  have 
acopo  for  improvement)  a.g.  the  inverse  Fourier  transformation  le  not  poeitlvo'Conatralnad 
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Ptnaanlv  movtna  objtctt 

(••beratery  axparlmanta  and  ootnputar  almulatlona  have  baan  oarilad  out  on  Imaalnf  a,, 
randomly  movinq  two  dimanalentl  objaet  at  vary  lew  llqht  lavala.  Carllar  axparlmanta  on  a 
ona  dlmanaional  objaot  indlaatad  that  objact  raconatruotlon  oan  ba  aohiavad  at  avaraqa 
alqnal  lavala  on  tha  ordar  of  i  dacaotad  photon  par  frame  and  30,000  framaa. 

rigura  i  ehnwa  a  aehematie  of  tha  axparlmantal  apparatua.  Tha  random  motion  waa 
aohiavad  by  tilting  a  mirror  linked  by  a  lavar  to  randomly  driven  leudapauKar.eeila.  Tha 
Imaging  Photon  Oetaoter  (IPDl  ia  a  miereehannal  plata/raalativa  aneda  device  with 
336  X  3S(  pixel  ooordlnataa,  although  In  tha  axparimant  the  pixel  area  waa  raetrletad  to 
33  X  33. 

To  Imago  a  randomly  tranalatlng  obiaot  wa  clearly  muat  average  aoma  quantity  that  ia 
tranilation*invariant.  Tha  autcenrraiation,  tq.d),  (or  energy  apaetrut.O  tatiafiaa  thla 
eendltion,  but  It  la  not,  in  ganaLal,  peanlbla  to  raeenatruet  an  ebjaet  unamblgueualy 
from  ita  auteoerrelation  funotlen  .  The  trlola  correlation  la  alee  tranalaticn-invariant 
and  la  theretcra  ideally  aulted  to  thia  problem  aa  flret  pointed  out  by  tartalt  and 
Wlrnltiar*°. 

All  caloulatlona  ware  carried  out  on  an  HPgtSt  daaktop  computari  tha  total  time  taken 
to  oompleta  tha  biapaetrum  of  3000  framaa  with  an  avaraga  of  13  dateotad  phetona  par  frame 
waa  approximately  30  houra  and  tha  raconatructien  from  tha  biapaetrum  waa  aeccmpllahed 
in  a  taw  minutaa.  The  raaulta  aru  ahown  in  Plgure  2.  Ilmpla  oentroidinn  prior  to  averaging 
yielda  a  blurred  raault,  ai  in  Plg.llbl  -  in  feet,  cantroldlng  can  peaaibly  ba  dona  in  a 
more  intelligent  way  by  grouping  togathar  all  framaa  with  tha  aama  number  of  datactad 
phetona  I  but  thia  hat  not  bean  attempted  yat.  Tha  raault  derived  from  tha  triple 
cerralation.  Figure  3(e),  uaing  3000  framaa  and  an  average  of  13  detected  photoni  par 
frame  (approxlmitaly  0.1  photona  per  eblaet  pixel)  oomparaa  favourably  with  the  original. 
However,  tha  raconatructien  ia  not  perfect  and  it  ia  believed  that  thia  waa  due  to 
imparfactlena  in  tha  fltfflald  raaponae  of  tha  IPO.  A  computer  aimulatlen  for  thia  eaaa 
gave  almoat  perfect  reeonatruction  with  3000  framaa. 

A  computer  cimulation  atudy  waa  carried  out  for  a  aaeond  obiaot,  oomprlalng  approximately 
130  plxala  with  a  pixel  area  of  133  x  133  (aaa  Figure  3).  Again,  3000  framaa  at  an  avaraga 
photon  rate  of  IS  photona  per  frame  (*■  o. i  photona  par  object  pixel)  ware  uaed.  A  aun3/lt0 
computer  and  tky  "Harrier”  array  preeaeaer  ware  uaed  to  generate  tha  data,  calculate  the 
biapaetrum  and  reeenttruet  the  object  from  the  biapeetrum.  Two  mathoda  of  calculating  tha 
biapaetrum  ware  uaadi  in  tha  firct  ( racenattuetion  ahown  in  Fig.Kc) )  .aalactod  planaa  ware 
computad  and  in  the  aaeond  (Fig  3(d))  a  tomographic  technique  uaing  33  anglea  waa  uaed. 
Daapita  the  very  low  photon  flux  of  approxlmatoly  0.1  par  objact  pixel  on  average,  the 
raconatructient  are  of  high  quality.  Tha  computing  time  for  Fig  3(e)  waa  approximately 
30  min  (compared  to  approx  30  hr  for  the  33  x  33  array  in  Fig.2),  i.a.  approx  I.Sa  per 
frame. 


ioeckla  imaaina 

apaekla  interferometry  givaa  a  good  aatimata  of  the  Fourier  modulua  |0(u)|  of  an  object, 
but  doaa  not  provide  tha  Fourier  phaae.  The  Knox-Thompion  method  glvea  aa^imatai  of  phaae 
dift'erancei  in  the  Fourier  plena  of  the  Image  aa  can  ba  aaan  by  taking  the  argument  of 
Fq .  ( 1 3 ) 1 

(u,  ,41^)  •  ^1*1*1*  *  *** 


where  g.  la  tha  phaae  of  the  image  apeetrum  and  Ajj  it  a  fixad  (email)  vaetor.  Tha  Knox- 
Thempaen  tranafer  function  ia  purely  real  (to  firat  order  )  and  thua  tha  phaae  of  tha 
object  apactrum  a^lji)  ia  almply  equal  to  that  of  tha  image  apeetrum  «{(ii). 

.Taking  tha  argument  of  Cq.(7),  it  it  clear  that,  for  a  tingle  image,  the  pnaaa 
Vt  '  ^ •  3u)  of  tha  biapaetrum  la  equal  to  that  given  by  the  Knox-Thompaon  method  plua 


(3) 


For  a  tingle  image,  f,(4u)  ia  a  eonatant  and  may  ba  taken  to  ba  lero.  Thua,  tha  Knox* 
Thompaon  quantity  la  iquTvalant  to  a  tingle  plane  of  tha  biapaetrum, 

Conaldar  now  tha  maaturamant  of  the  antambla  avaraga  of  Eqa.(3)  and  (3)  for  tha  triple 
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corraUtlon  and  Knoa-Thonpion  matheda  raapaetlvaly,  whtra  tha  lina9aa  hava  random  diaplaoa- 
mantl  dua  to  random  tllta  in  tha  ineldant  wavafronc.  Tha  tripla  eerralation  ia  Invariant 
to  trantlatlon,  wharaai  It  la  naeaasary  to  axaotly  oantroid  tha  Imaga  Intanalty  to  obtain 
tranilatlon  Invarlanoa  In  tha  Knox«Thompaon  taehnlqua  In  prlnolpla,  thla  can  ba  dona 
at  hlph  light  lavala  but  net  whan  thara  ara  only  a  taw  dataotad  photona  par  frama.  Howavar, 
tha  raaulta  e<  eemputar  almulationa  Indieata  that  >  <00  dataotad  photona  par  frama  ara 
probably  raqulrad  ter  imaging  all  but  tha  vary  almplaat  objaeta  and  thua,  In  praotioa, 
tha  KnoN-Thempaon  mathod  la  affaetlvaly  trantlatlon  Invariant  provided  tha  data  la  eantroldad 
prior  to  preeat ting. 

Tha  algnal-to-nolaa  ratio  of  a  maaauramant  of  a  alngla  Plaaa  of  tha  biapactrum,  Bq.(7), 
la  lata  than  that  on  tha  Knea^Thempton  aroaa>apaotrumi  tq.(t)  i  lewar  ordar  momenta  ean 
alwayt  ba  aatlmatad  mera  aeourataly  than  higher  order  enaa.  On  tha  othar  hand,  additional 
phaaa  aatlmataa  ean  ba  obtained  from  othar  planet  in  the  bltpaotrum.  Thaaa  phaea  aatlmataa 
ara  not  all  atatlatleally  Indapandant  and  thalr  error  Ineraaaaa  at  one  mevet  off  tha  main 
diagonal  u.  ■  -Uj  of  the  blapaetrumt  thla  faet  ahould  ba  taken  Into  aeeount  In  tha 
raeenatrueilen  algorithm,  Raaulta  of  eemputar  almulationa  Indieata  that  tha  two  taehnlquaa 
give  approximately  the  tame  algnal-to-nolaa  ratio  for  polnt-llka  objaeta. 

Tha  principal  ndvantaga  of  raeonatruetlon  from  tha  bltpaotrum  la  that  It  avolda  tha 
problamt  etuaad  by  email  valuta  of  tha  nodulua  in  the  Knox-Thempaon  mathod,  Thaaa  ooeur  at 
pelntt  In  tha  Fourier  plana  where  tha  medulua  la  tare  or  amail,  and  tha  fixed  value  iu 
doaa  net  allow  almpla  waya  around  thla  problem,  Tha  affect  In  praotioa  la  to  limit  tno 
upper  apatial  traquanoy  of  phaaa  raoovary  and  la  objaet-dapandant. 

Computer  almulationa  of  tpaokle  imaging  hava  bean  carried  out  to  compare  tha  Knox- 
Thompaun  and  triple  eerralation  taehnlquaa.  In  a  typical  alm'ilatlon,  20,000  photon-llmltad 
121  X  <21  apaokla  framaa  ara  ganaratad  In  tha  following  wayi  (a)  ganarata  a  oorralatad 
Qauaalan  phaaa  preeaaa  of  appropriate  varlanoa  and  eerralation  length  In  tha  pupil,  lb) 
Fourier  tranaform  and  taka  aquara  modulua  to  ganarata  Intanalty  point  apraad  function 
(apaekla  Image  of  point  aourea),  (e<  convolve  with  object  Intanalty,  Id)  ganarata  photona 
via  a  Folaaon  random  number  generator,  Ualng  tha  Sun}/<(0  and  array  prooaaaer.  It  takaa 
approximately  <a  to  ganarata  ona  aat  of  photon  eoordlnataa  (l.a,  one  image).  Figure  4 
ahowa  aoAia  eemputar-ganaratad  etallar  apaekla  pattarDf,  tor  a  D  •  2m  talaaeopa,  t.  ■  <0  om, 

1  ■  SOO  nm,  apatial  fraquanoy  eut-eff  ■  <2.4  are  aao  . 

Figure  $  oomparaa  raeonatruetlena  of  an  aatarold-llka  objaet  via  the  two  taehnlquaa 
ter  <S,000  framaa  with  an  average  of  <00  dataotad  photona  par  framai  aeaumlng  a  wavelength 
bandwidth  41  •  20  nm,  axpeaura  time  of  <0  *a  and  a  D  ■  2m  talaaeopa,  thla  corroapendi  to  an 
integrated  magnitude  of  m^  -  1$  for  tha  aatarold,  whoea  extant  la  approximately  O.S  are  aae 
by  < ,0  are  aao. 
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Plaura  i  lohanatie  of  aqulpntnt  tor 
attpartmanta  ualng  a  randomiy  moving  ebjaet. 


n 


u 


Randomly 
tilting 
,  mirror. 


OllfuMly 


A  B  C 


PlQura  2  CHparlmantal  ratulta  on  tha  imaging  of  a  randomly  moving  objact  at  vary  low  light 
lavaia  lavaraga  of  19  dataetad  photons  par  trams,  2000  frtmsa),  (A)  •  Imtga  of  stationary 
object  (11  •  csntroidad  imaga  (C)  •  imaga  raconstruetad  from  tha  tripla  cerralation 
msthod. 
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Flauty  a  CompuMr  •iinui«tlon 
raiuTtt  (or  Imaplng  •  randomly 
moving  ebjaet  (averaga  of  iS 
dataetad  phetena  par  trama. 
3000  franaa) . 

A  •  imaga  of  atatlonary  ebjaet 
•  -  oantreidad  Imagt 
C  -  raeonatruetlen  from  bl* 
apaetrum  using  •alactad 

pianaa 

D  <•  tiaing  radon  tranaform 
and  biapaetrum 


Fieura  4  Computar-gansratad  atallar  apaokls  pattarna  of  a  point  aouroa  I  left)  and  a 
uniform  disc  iright). 


24  /  tm  Vtl  10$  Inrtrit  $nttnni  In  Opilii  ll$$7l 


Cemputar  titnulatlon  o(  intpinp  in  ••t«roid-llk«  objiot  through  turbulanei. 
Tnoien  oiioi  1]|000  fTr.-.nat  with  in  ivirigi  at  100  dotactid  photona  par  frima  wai  uaad. 
(A)  •  dl  ff  not  ion-1  Iml  tad  Imgi  at  objaet 
(•)  -  long-anpoaura  imioa  through  turbulanea 

(C)  -  raoenatruotion  uilng  tha  Knox-Thonpien  itathod 

(D)  -  raoenatruetion  uiing  tha  biipaetrum  mathod. 

At  thia  light  loval,  mora  Crama  ara  raguirad  to  yiald  an  aeeaptabla  raconatruetion. 
Bnamination  of  tha  apaetruit  a(  tha  raaonatruetiona  C  and  0  ahewa  that  tha  biapactruit 
raeonatruction  la  auparler  to  that  of  tha  KV  aathod. 
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REF  5 


Spacklt  imiflng  taehnlquat 
J  C.  Dainty 

Blacktcc  l.aboratory.  Imparlal  CollaB#.  London  SW7  2B7i  Ux 


A|fZS£CX 

lavaral  prepoaaj  taennlquat  of  fooal  plana  imaqing  through  turbulanea  ara  daaorlbad  and 
comparad.  Tha  taehnlquaa  arai  ahift-and-add,  axponantial  (Jltar,  Knox-Thonpaon  (oroaa- 
aptrtrum)  and  trlpla  eetralatlon  Particular  attantion  la  glvan  to  tha  parfornanea  of  tha 
mathoda  at  axtramaly  low  light  lavall  of  a  few  dataetad  photona  par  frair.a. 

INrHOPWCTIQN 

bpoexla  intarfaroRiatry,  aa  invantad  by  Labayrla'.  prnvldai  an  aatlmata  of  tha  diffraction 
lliTiltod  auteeerralatlon  or  anargy  apactrum  of  an  objaet.  in  a  llnaar.  laoplanatlc, 
ineeharant  imaging  ayatam,  tha  Imaga  intanalty  K]^)  la  ralatad  to  tha  ob^act  Intanaity 
o(x)  by  a  convolution  ralatlonahlp, 

1 1 X)  •  oljl  •  plj)  ,  ( 1 ) 

whara  *  danctaa  convolution  and  p()^)  la  tha  Inatantxnaoua  point  apraad  function  of  tho 
optical  ayatam.  Whan  imaging  through  turbulanea,  p(x)  ia>a  random  apaokla  pattern  whoaa 
atatlatica  ara  almllar  to  thoaa  nf  laaar  apackla  pattarna^.  In  Foutlar  apaoa,  Cq.P)  la 

I(y)  -  0(u)  .  P(ji)  . 

whara  tha  uppar>caaa  functiona  danota  tha  Pourlar  tranafor.ta  of  tha  lowar*caaa  onaai  tha 
quantity  P(uI  la  tha  inatantanaoua  taIaacopa«atiroaphara  tranafac  funetlon. 

In  conventional  long-axpoaura  Imaging,  one  forma  tho  avaraga, 

<  Kill  )  ■  0(21  .  (  p<2t  '  (3) 

and,  aa  la  wall>known,  tha  avaraga  tal9aeopa>atmoaphara  tranafar  function  CP(n)>  dropa  off 
rapidly  for  angular  apatial  fraquanelaa  greater  than  r  /l,  where  r.,  tha  Priae  paramatar, 
la  typically  lO  cm  at  a  good  obaarvlng  aita. 

In  apackla  Intarfaromatry,  one  forma  an  avaraga  of  tha  anargy  apactrum, 

(  l‘^’(ul  >  ■  (  Ku)  I*(u)  ) 


-  o‘**(u)  . (  P*^’(u)  ^  .  (41 

III  Eq.(4!,  tha  auparacrlpt  danotaa  that  a  aocond  ordar  correlation  la  being  parformad, 
which  in  image  apaca  coulo  be  written 


f  l  lx,  )  -  <  J  iljjl  Kji  *  >ti  >  ^  •  •*  > 


At  originally  thown  by  Xorff  ^  ,  tha  apackla  tranafar  function  (  P  ^  (y)  ^  eontalna  a 
dlf fraction-limited  component  and  thua  allowa  tha  aatlmatlon  of  tho  anargy  apactrum  of  the 
Object  up  to  an  angular  apatial  fraquonoy  of  0/1,  whara  D  la  the  taleacopa  diameter. 

At  a  raault  of  fq  (41,  apackla  Intarfaromatry  haa  thrao  propartlaa  that  have  contributed 
algnl f Icantly  to  Ita  auccata  In  aatronomyi 

(11  It  provldaa  dif fraction-limited  Information  about  tha  object  atiuctura  , 
lUl  that  information  la  parturbad  by  a  linear  tranafar  function,  and 
'll; I  tha  image  anargy  apactrum  can  be  avaluatad  at  vary  low  light  lavala  ualng 


DA  I  I 


32 


CiLO  Ui 


phovon  corralition  tpohnlquap  *  . 

Ra^trding  point  (il)i  It  it  not  •tttntlil  that  th«  traniftr  proeitt  bt  llrttr.  but  it  it 
utjtUy  dotiraola  and  in  any  cttt  th*  imtgino  procott  nutt  bt  wtll-undt-ttood  and  hava  an 
itiuartt  loiution  ao  that  ob]aet  data  oan  ba  found  quantitativaly , 

Clairiy,  naithar  rq,(4)  or  IS),  provtdaa  an  aatiiaata  of  tha.objaot  .nap  o(fl)  dlroetly. 
Hathodt  of  trying  to  attimata  o(ji)  freffl  ita  anargy  tpaotrum  0  ’(u)  ara  raviawad  in  thaaa 

Procaadingt  by  rianup  and  altcwhara?.  Tha  prebloA  it  that  of  tatlmating  tha  phata  a(u)  of 
Ol^ii  glvan  tha  aquarad  nodulua  ■  |0(u)|'.  Although  a  oontldarabla  aaiount  of~affort 

hat  gona  into  tolvlng  thit  problaoi  no  oanaraT  aolution  axlttt.  rurtharmora,  apaekla 
ptttarnt  contain  information  that  it  additional  to  that  givan  by  tha  anargy  tpaotrum  and 
thit  thould  ba  axploitad  whan  finding  tha  objaot  map. 

In  tha  following  taetion  wa  raviaw  briafly  four  prepoaad  mathoda  of  apaekla  imaoino.  aach 
of  which  procattat  tha  raw  apaokla  data  in  a  diffaront  way  to  l.abayria't  apaekla 
intarfaromatry ,  Hewavari  it  ia  daairabla  that  any  propoaad  taehniqua  of  apaokla  Imaging 
rataint  tha  tdvantagaoua  propartiaa  li)  •  !iii)  of  tha  original  apaekla  taehniquai  at  wall 
at  giving  an  attimata  of  tha  objaet  map.  Property  (ill),  tha  ability  to  work  nt  vary  low 
light  lavalt.  it  particuiarly  important,  tinea  jutt  about  any  mathod  workt  whan  photon  i.oiaa 
it  naglaetad. 


POUR  PROPOirP  TICMWIQUEI  OP  IPtCXli!  IMAaiNO 


A.  $nift -and-add 

In  thit  mathod,  firtt  propotad  by  l.yndt  at  al‘  and  eonaidarably  axtandad  by  latat^'*  and 
othara  ,  ona  ivortgat  eantreidad  varalona  of  aalaotad  bright  tpaoklat.  At  high  light  Iwvalt, 
tha  mathod  would  htva  llnatr  trant'ar  lunation  it  ciaoh  "apaekla”  ware  a  randomly  ttanalitad, 
linoarly  daoradad  imapa  of  tha  orieinti  object.  Thia  ia  not  tha  eaaa,  at  thown  by 
Hunt  at  il  ,  ind  tha  thi f t-anj*tda  mathod  it  not  linear  at  high  light  lavalt.  Tnara  it 
alto  tha  problem  of  irianuifyinj  bright  tpteklaa  in  photon-limitad  datt, 

Raeantly,  wa  have  anaiytad.thif t-and'tdd  imagine  ;t  very  low  light  lavalt  for  tha  otta  of 
E  randomly  trantltting  objaet  (whota  imaging  la  linear  at  high  light  lavalt).  Tha 
photon-iimittd  etta  it  alto  fundamanully  non-iinaar.  if  <  Ou(u)  )  danetat  tha  Peuriar 
tpaotrum  of  tha  tvaraga  of  all  oantroidad  framaa  that  oontiin"aKactly  M  dttaeted  photont  and 
Kul  danotaa  tha  normalitad  Pouriar  apaetrum  of  tha  high  ligi.t  level  intanaity  of  tha 
original  (unihiftad)  imaga,  than 


/  O^iul  'I 


I  u) 


U) 


?njt,  for  axampla.  for  N  •  1 


^  Oj(u)  )  *11  I(j) 


|7) 


that  it,  tha  tpaotrum  of  tha  ahl f t-and-idd  image  ia  proportionui  to  tha  energy. apaetrum  of 
tha  true  image  Fquation  (6)  hia  intaraatlng  implicationa  for  phaaa  retrieval '  . 

B.  FxDonantial  filtatinq 

For  modulua-only  raconatruetior.  in  ona  dimanalon.  tha  aoJjtion  ia  ambiguoua  btetute  tha 
complex  tarot  of  tha  anargy  apaetrum  of  tha  object  include  uoth  tha  taroi  of  tha  original 
ooiact  and  their  invaraai  Tha.gojact  oan  ba  raconatructad  uniquely  if  ita  tarot  can  ba 
idantiflad  eorraetiy  and  Walker  '  ahowad  that  thit  oan  bt  dona  uaing  an  axponantlal 
fiittr  in  tha  imaga  domain. 

Tna  avartga  anargy  apaetrum  of  tha  imagt,  <  l'*’(u)  )  •  (  |r(u)  |^  >,  it  found  in  tha 

uauil  way,  A  ateond  anargy  tpaotrum  i  :''‘’(u)  >,  it  found  by  mulTlplylng  tach  framt  l(k) 
by  an  axponantial  ,  axp( -lanx ) ,  whara  'a*  la  ?  eonttant.  It  car.  ba  ahown  that 


I '  (u)  )  ■  0'  ‘  (u)  .  f  P'  ‘^’(u)  1  , 


"I 


n 
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»o  th»t  Iron  th»  oriflJnal  data  oan  axtraet  th*  ancrgy  apaetrum  of  tha  objaet  and  that  of 
th»  axpontntla) ly  flltarad  objact.  In  prlnclpla,  thaaa  togathar  contain  auffleiant 
information  to  raoonatrjct  tha  Objact  Intanalty.  In  practlca,  ona  haa  to  uaa  a  itodiflad 
varalon  of  tha  rianup  algorithm’  to  raeovar  tha  objact  map  and  no  aatiafactory  modi f teat  Ion 
vKiata . 

C  Creat-apaetfum  method 

,„lhl8  method  waa  flrat  propoaad  oy  Kno*  and  Thorapaon"(KT)  eorralatlon, 

1 ,aji) ,  and  Ita  Fourlar  tranaform  (with  roapact  to  ij,),  1  ^a^lnad  by 


1 lx,  ,du)  •  J  l*(x)  1(15  ♦  X,)  axplJtldu.x)  dx  191 


and 


•  llji,)  l*(a,  *  Ajil  .  CO) 


wh'Ta  la  a  fixed  amall  vector. 

Whan  applying  Fqt.(9l  and  (  lO)  to  apackla  Imaging,  an  avaraga  la  taltan  ovar  many  apackla 
fra.maa  and  it  can  bt  ahown  that' 


t  I  Ig,  ,4u)  ^ 

(ID 

where 

0*'*'^’ (u,  ,4u  1  • 

0(u, )  0*(u,  ♦  4u) 

(U) 

and 

('  (Hi  -4ji)  > 

••  <  P**’(.H,  )  >  . 

1  <  p(4^)  >  1’  , 

(IJ) 

whara  <  p)  la  tha  aaalng*ltmttad  tranafar  function  It  followa  from  fq  (U)  that  Uul 
mutt  be  taaa  than  tha  aaalng-limltad  angular  fraquancy  of  r  /I  for  dlf f ractlon>llinltad 
information  to  be  retained  In  fq.()1l.  It  alao  followa  that  tha  phaaa  of  tha  KT  tranafar 
function  la  aero  and  thua,  taking  the  argumtnt  of  Fq  (11), 

(u,  ,aul  «  '  ’’C 

whar*  la  the  phaaa  of  and  o^iu)  '.a  the  phaaa  of  0(u)  . 

Tnua  the  KT  technique  glvea  phaaa  dltfaratoa. information  In  the  apactrum  of  the  object 
and  the  phaaa  haa  to  be  found  by  Integration  *•  avoiding  phaaa  dlalocatlona. 

D.  Triple  correlation 

Thla  method  waa  firat  applied  to  tpaekla  Imaging  by  Welgglt’^*^^  and  la  daacrlbad  in 

Tha  triple  (auto-)  eorralatlon 

(IS) 


more  detail  alaawhare  In  thaaa  Procaadlnga  and  In  ravlewa  j’. 
and  ita  Fourier  tranaform,  tha  blapaetrum,  are  defined  by 


i‘^’()I,, 


m 

I  ■*(»)  i (4  ♦  ai>  ca  *  sj)  da 
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and 


1 '  ^  ( J  ■  .  J 


I  'dj' 


I  ^  fH)  ■*  Jij*  • 


(16) 


In  ipaekl*  iir*9in9i  tha  tv«i«9«  Bliptctrum  of  th*  iiii*9«  li  ralitad  to  that 
throu9h  •  bitpaetruir  tranifar  function  <  ^  (u^iUj) 


of  eha  objact 


(  :*^’(u,,uj)  ) 


(17) 


At  with  tha  KT  tacnniquai  tha  blapcotrum  trtntftr  function  it  non-taro  in  naqnituda  ind  hat 
taro  phaaa  up  to  tha  diffraction  .Uirit  of  tha  taiaaeepa.  Uain9  tnri  fact  and  taking  tha 
argumtnt  of  Fq.lIC),  wt  find  that 


» '  ^ '  *ili 'Hj  *  *  *o’-l’  *  *o*-1  *  iij’  *  '  *'*’ 

whara  9'^’  la  tha  biapaetrum  phaaa.  A  racuraiva  algorithm  can  bj  uaad  to  raconttruet  the 
phaaa  of  tha  ob>ct  tranaform,  9g(u),  from  tha  biapaetrum  phaaa. 

eoMPABiaoN  or  cpoia-aptcrrjM  anp  aiapEciauM 

If  rq.(l6)  for  tha  biapaetrum  of  a  alngla  framt  ia  writtan  with  •  du.  than,  for  a 
tingle  'plana'  of  tha  biapaetrum,  ‘ 


l'^’(u,,4u)  •  Ku,)  I*(j,  *  du)  Kdu) 

•  l''^'''’(u,,du)  Kdu)  ,  (19) 


and  thut  thli  plana  of  tha  biapaetrum  aquila  tha  ereaa-apactrum  multiplied  by  an  image- 
dependant  number  Kdu)  Taking  tha  argument  of  Fq  (i9), 

■  9'**'''*(u,,dji)  ♦  9U^)  (JO) 


Tha  phaee  of  the  biapaetrum  le  invariant  to  a  trantlation  of  the  imiga  (thia  ia  clear 
from  rqi.()5)  and  (I6)  I  wharaaa  that  of  the  oroea-apaotrum  ie  not  thift-invariant  .  Whan 
implamanting  the  KT  algorithm  in  practiea,  it  la  nacaeaary  to  centroid  each  frame  prior  to 
procaeaing  and  thit  ia  equivalant,  in  tha  nolca  aanaa,  of  calculating  a  aingla  plane  of  the 
biapaetrum. 

Fxaet  analytical  axprattione  Cor  tha  KT  and  biapaetrum  trantfar  functiona  are  difficult  to 
obtain  and  inataad  wa  have  ueed  Monte  Carlo  tachniquea  to  compare  them.  Figure  1  comparaa 
tha  two  tranafar  function!  for  equivalant  aubplanea  for  raaliatlc  atmoinhaclc  paramatare. 

Signal -to-Noiaa  Patioa 

Both  tha  eroae-apactrum  and  triple  correlation  are  complex  quantitiaa  and  thua  both  the 
iignal-to-nolae  ratio  of  tha  ngdulue,  SNK  ,  and  tha  phaaa  error  F  are  of  intaraat.  Tha 
method  of  Goodman  and  Balahar  '  for  finding  BNP  for  the  energy  apactrum  can  be  axtandad  to 
tha  croea-epactrum  and  biapaetrum  for  both  8NR  '"and  F  *'.  Figure  2  ahowa  INK  and  F  in 
tha  frequency  domain  for  an  aateroid  type  of  object  and  an  imaging  ayatam  with  a  unit 
tranafar  function  (1  a  no  atmoapharic  turbulanca),  for  a  mean  number  N  •  100  photona  per 
frama  for  equivalant  croae-apactrum  (KT)  and  biapaetrum  aactiona.  Note  that  the  8NR  ie 
lower  (and  (T  higher)  for  the  biapaetrum  eaaa,  confirming  tha  general  raault  that  arrore  are 
alwaya  higher  on  higher  momenta. 

Note  that,  in  practice. 


(21  ) 
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as  iisuired  by  Niatnson  and  Fapaliolioa^^  in  thair  analyain  of  tho  affact  of  photon  noita  on 
tna  KT  alqorlthtn. 

In  tha  praaance  of  atmoipharic  turbulanca,  tha  axpreaaiont  for  tha  aignajito-noiia  ratio 
p*r  frarrt  in  th«  cro«<*spcctrum  and  blspvctrunt  can  bacon^a  quite  complicated  At  hi9h 
light  levele,  then 

8NR 1  .  IjijI  ‘  ^o^' 

and  122) 

>1  .  Uu|  <  r^/x 


For  liu;  >>  r^/X,  SNf"*'''' 


0 


wharaaa  for  |ijj|  >>  r^/X 


SKR 


TCI 


2t‘^’(u,,Uj) 


(23) 


whart  T*''*  li  tha  normaliiad  ovarlap  araa  of  n  pupila. 

At  low  light  Itvala,  in  tha  important  ragion  14u|  and  lu.|  <  r  /X,  both  algnal-to-noiaa 
ration  ahow  a  linair  dapandanea  on  tha  avtiaga  numBar  of  pBfitona  par  apacltla,  joat  lilta  tha 
powar  apactrum  cata.  For  Uui  and  |u, |  >>  r./X,  tha  croai-apactrum  hat  a  ligntl-to-noiia 
of  ataantitlly  laro,  wharaat~tha  biapictrum  naintaina  a  imall  valua  which  it  ttrongly 
dapandant  on  tha  avartga  numbar  of  apac)(laa  par  frama. 
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OBJECT  RECONSTRUCTION  FROM  RHOTON*LtMlTED  CENTROEDED  DATA 
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L.O.  da  FMIm,  M>  NorihaoU,  BJ.  Bramaa  aad  J.O.  Dainty 


Opllea  SaetloB  *  Blackatt  Uborateiy 
ImpatUl  CaUaia  •  SWT  3B2  •  U.Ki 


ABSTRACT 

Caatroldlaf  U  InvaallfaUd  ai  a  ilmpla  and  eemputatlonalN  lut  tachalqua  of  Inan  taconatiuctlon,  at  low  Uflit  laval,  of 
a  raodomly  traoalaUui  imaca.  Tha  ^taetad  framaa  aia  aoitaa  by  tbalr  nuatbat  of  pootona,  caatroldad,  aaparatod  avataiaa 
parfonnad  and  tbau  compand  with  Iba  uiual  way  ol  aantnldlAi  framaa.  An  alfoilthm  for  nltlavlai  tba  pbaia  for  ooa> 
dlmanalooal  eanirutdad  Imagloi  la  pnaaatad  aad  compiilar  almulaiad  data  la  uaed  to  taat  tba  theory  and  tha  laeonatrucUon 
taehniqua. 


1  .  TNTRQPUCTIQN 

Tha  dlffraetlon-Uniltad  lma|a  naolutloa  of  •  |rottad«baaad  4  raatn  talaacopa  oparatlng  at  400  nm  la  approxlmataly  0.036* 
(are  aaeond).  Atmoapharlc  turbulanca,  howevar.  makaa  tha  lma|a  of  an  imraaolvad  atar  broadar  ^  1.0*  (one  aacond  of  are)  aad 
axhlblt  a  (ranular  atruetura  raaambllni  a  apaekla  pattaro. 

A  Bumbar  of  taehfiiquaa  have  baao  davalopad  to  re'.rlava  tha  dlirraetlon.Umltad  lma|a  by  racordlqj  ahert.axpoaun  framaa 
of  thaaa  apaekla  patlaroa.  Tba  SK1FT-AND>ADD  (SAA)  matbod**"**  la  osa  of  tha  propoaad  taehalquaa.  It  tallaa  oa  tha 
propoaltioB  that  aaeh  frama  eoaalata  of  naoy  dlatortad  r^Ueaa  of  tha  trua  Imaja  and  that  an  Improved  aatlmata  of  tba  lma|a 
eould  ba  found  by  aupatlmpoalnc  thaaa  dlatortad  rapUeaa.  Tba  auparpoaitloB  u  earrlad  out  by  eoBaldarlsf  that  tba  brlibtaat 
part  of  a  apaekla  imaie  la,  moat  Ukaly,  a  dlatortad  varaioB  of  tba  brl|blaat  part  of  tba  true  ImMO.  la  eaeb  frame,  tba  brlibtar 
apaeklaa  are  found  aad  by  tblftlng  and  auparlmpoalof  them  at  tba  caatia  of  tba  frama  (SUirT>AND*ADD)  tba  ImparfaetloBt 
of  tba  Individual  ImagM  tood  to  average  out. 

At  vary  low  light  Itvala,  for  example,  tha  mean  Bumbar  of  phetoBi/frama  7^  <  lOpbotona/framt,  however,  eaeb  of  tbeae 
dlatortad  Inutgea  may  have  only  a  few  pbetona  and  no  ‘‘bright*  tpeekla  eon  ba  ehoaaa  in  order  to  Implament  the  SAA  technique. 
Aa  a  ftrat  approach  on  bow  the  ahlft<and-add  (SAA)  ta^tnique  behavaa,  at  very  low  light  levala,  a  tbeoratlcal  atudy  la  praaentad 
foeualog  on  ceatroidlng  photoo-llinltod  data  amltiad  from  a  randointy  tranalntlng  ImAge  .  A  nlatloaahlp  batweeu  tbeae 
cantroldad  Imagea  and  tha  atationary  normallaad  Image  /(r)  aa  wall  aa  a  pbaaa  reeenatruetien  algorithm  la  praaentad  for  tha 
caaa  of  OBo-dlmanalonal  object.  Computer  aimulatad  photon  data  amlttad  from  a  binary  atar  lyatam  la  uaed  to  uaaai  tba  theory 
and  tha  raconatruetlea  algorithm. 

3  .  THEORY  OF  CENTROIDING  DATA  OF  RANDOMLY  TRANSLATING  IMAGES 

T!b  ‘froeia*  a  randomly  moving  image,  pbotoaa.  that  ara  all  auppoaod  to  be  eaiaaatad  from  tba  aama  randomly  tranalating 
Image,  an  deteetad  during  a  aariaa  of  abort  time  iotarvala(framaa).  To  retrieve  tba  atationary  Imago  one  ihould,  for  each  frame, 
•hlft  the  photon  vectori  by  the  *010001  that  tha  trua  centroid  of  the  Image  >*  dlaplncad  la  teapaet  to  tba  centra  of  tha 
fiama  and  average(add)  over  all  tba  framaa.  Tiia  true  centroid  and  hence  the  trua  ablfi  vactor  ct,  howevar,  la  unknown  and 
another  ahift  va.'tor  Rt  rouit  ba  detarmlnad  (tea  Fig.  1).  Aa  an  aatimator  for  tba  true  eantroid,  tba  centroid  vauior  of  the 
datactod  pbotona  can  be  evaluated,  all  tha  photon  veetora  an  ibaa  abiftad  by  tbla  aatimator  and  aa  average  image  of  many 
auch  framaa  la  formad. 

Tha  nlatlonihip  beiwaan  tbla  aatlmatad  Image  and  tba  normallaad  atationary  image  /(r)  can  ba  derived  aaaumlog  that  the 
N.photon  data  ift(x)  datactod  on  tba  frame  it  modaladl*^  aa  art  iohomogtntou.i  Peiaton  procota  mathama:  tenlly  datcribed 
la:  ^ 

de(x)- 23 

i-i 

Let  £t(u,Xi,...,Xjv)  br  tha  Tburier  tranafonn  of  da(x) 

.5a(u,x . .  ■  2  7^f(x  -  Xj)  tzp(-i2rn  ■  x)dx 

IV 

-£e*p(-i3ru.x,)  (i) 

1-1 
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Fl|ur«  1  •  Typical  hiat  ahowlai  th«  tru*  aaii  tba  aallmaiad  ihifl  vaclora  •  ct  and  Kt 
The  N  photon  controld  vector  Ro  of  thi  h'*  frame  (Flguta  1)  b  doftnod  u 

Rt-yE*/ 

SabatltutUf  In  tormi  of  th«  eontrold  ooordloatoi  r/  ■<  x/  •>  Rt  (Fl|.  1)  wo  have 


(9) 


Afu.Xi . Xw)  ■  up{-i3wu  •  R»)  2^  ixp(-i2iru  ■  r>) 

/■» 

■  o»p(-<2»u  •  R»)  3j(u,ri,...r//) 

wbaro  3f(u,ri,...r^)  ■  f{(u,X|,...x/ir)'ii.tho  Fburior  Iraaafom  of  tbo  eantroldod 
Tbo  relation  between  tbe  eentrolded  and  noo'oeotroldid  data  epectrum  b  tbereforei 


(4) 


ft(u.Xi . xs)  -  e»p(+<J»u  •R*)i5i(u,Xi,..jtrf) 

N 

m  e<p(-l-i3ru  •  R*)  eopf-ilru  •  Xf) 
/■I 


<  e»p(+i3r«  •  ^). . e*p(-i3*u  •  |l  -  ^|x;) . e»p(+i2*u  • 


(9) 


Due  to  tbe  dependence  of  f{(u,xi,...,x/v)  on  tbe  photon  coordlnaiaa  Xj  tbe  eaeamble  averafe  of  tbene  X't(u,X|,...,XAr) 
can  be  performed  taklnf  Into  account  that  the  x/e  are  Independeai  random  vanablee  and  that  eoneequently  tbe  J^nt  probability 
denelty,  p(xi . xiv),  of  deiectlni  N  photone  between  ceordlnaiee  . . .  and  X|  -t-dbi, . .  •f  ixn  b  the  product  of 


'r<r 
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Ih*  liidtpwulwit  pcob4bUity  dnulUM  of  hoviaf  on*  photon  dotoetod  botwoan  xi  nod  X|  4  dxi  tlmw  tha  probobUity  dooiity  of 
dotoelini  utoihor  photon  bitwou  xi  ond  xt  f  dxi  tod  lo  on.  Htaco 


p(xi,.,.,x/y)dxi  ...dxw  -jKxi).p(xj) . p(XA/)dxj  ...dxjy 

whom  p(x>)  U  oqutl  to  tht  normnllnod  inUatity  Zi(x/)  nt  oootdlsnto  xy  and  ftamo  k  l.i.i 

P<*i)  -  Jt(xy) 


(») 


(T) 


'  Lot  ey  rtprount  tht  trut  oontrold  ihlft  voetor  at  ftamo  k.  Tht  furotioa  Xt(xy)  at  framt  k  and  eoordlnato  xy  U  roUtod  to 
tho  otatlonary  normallitd  ima|t  /(xy)  (Fl|.  I)  by 


At  a  contoqutaro 


X*(xy)  -  /(xy  -«t) 


»i(Xy)-X»(Xy) 

-/(Xy-Ct) 


(8) 


(8) 


whom  p(xyiet)  It  a  norfflaUitd  probahUlly  doulty  that  dtpondt  both  on  xy  and  on  tho  random  truo  etntrald  vaetor  ct. 

Tht  avtrait  of  3t(u,Xi . xh)  ovtr  tht  tottmbit  of  framti  hat  tharcfora  lo  bt  dont  lo  two  tltpti  Aitl  avtra|ln|  ovtr 

tht  dtitelad  photon  euordloattt  xy  and  than  ovar  tha  trna  eanttold  vaetor  et.  Lait  rapratant  by  l^{(u,et)  tht  raault  afiar 
•.  arailni  ovtr  Xi,...,x/r. 


J5{(u.et)  ■  (i>I(u,x, . ^ 

■  f p(*i.«t)-..F(xw,«t)I^I(u,Xi,.M,xw)dXi....dXN 

.Jh..  -et) . /(Xw  -ct)5l(u,x . . . dXff 


(iO) 


Subilllulinf  (S)  bio  (10)  ooa  gat: 


5;(u,c»)  ■  ^  /  /(xi  -et)tip(+i2)ru  -  ^)<ixi . 

y-iJi, 

. ^  /(*/-ea)t»y>(-t2au-xy[l-l/Ari)dxy., 

••m 

. J* I{xn  -  Ct)«rp(+i2nu  •  ^)dxN 

Tha  Inlagralt  In  (11)  can  bt  tvaluattd  tatlly  giving: 

Jn.,  -  Ck)tlf{-i2wa  ■  Xy)dxy  ■  tap(-i2r J  •  oj)  7(tf) 


(11) 


(12) 
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«rht(«  7(d)  ii  ih*  Fburi«r  tru*ii»rm  of  /(r)  «ad  d  mb  aitlMr  b«  BfuBl  to  -yx/S  or  to  •fu  [1  -  1/7V)> 

Uilni  mult  (13)  IbtJ'*  Urm  of  tha  luffl  io  (11)  U; 

■  7(-u/A^)  a»j»(+i3»^  •  at) . Tijull  -  l//r))a«p(-i2»  -^-u  •  e*) .....  7(-u/.V)  a4rp(+i8»^  •  e*) 

■  7(u[l  -  \IN\) m^-i3»  u  •  e*)  [f(-u/jy)]^  *  |a*f(+i2r •  «*)] 

■  7(u[l  -  UN])  [7’(-u/.y)] ^  ‘  a*p(+i3»|i^j^  -  u  •  c») 

■  7(u|l  -  l/AT])  [7(-u/^)j^”‘  0») 

which  U  IttdapaadaBt  of  J  tad  at. 

Tha  avBluBtlon  of  ii*  cum  ovar  J  mo,  tbaralara,  ha  parfarmad  atral|htforwardly  |1v1b{i 

-  iy7(uli-l/Ar))  [r(-u/A^)]''"‘  (U) 

Tha  mult  io  (14)  U  ladapaaitot  of  at  ood  haneai 

^•(u).  (^l(u,ct))^^ 

•  Ar7(u(i-1/Ar))  (7(-u/y)]''"‘  (i») 

It  It  coovaolaot  to  loModuca  a  BormBlliad  ipacttum  9;^(u)  dafload  m 

tylaldlni  tha  lioportaot  mult 

0;y(u)  -  7(u(l-l/Ar))  [r(-tt/ltr)]''"  (1(1) 

Expmaloo  (14)  U  a  ralatlooihlp  batwaao  tha  oormalliad  ipaetruot  of  tha  ctatloDary  Imafa  7(u)  aad  tha  apactrw  $n(u) 
o(  tha  caotroldad  avara|a  of  lhaaa  iramaa  cootalalof  axactly  8  photeoi/frama.  Tha  lodepaodaBca  of  (la)  oo  ct  ihowi  that 
4^jv(u)  ii  a  traoalatlnf  lavariaat  quaatity  that  dapaadi  only  oo  tha  ipaelruin  T  aad  oa  tha  aumbar  of  pholooi/fraina  N.  Aa  a 
eoaiaqueoca  the  quaatity  ^iv(u)  U  the  Mroa  althar  for  tha  raadomly  movlat  Imaia  aa  wall  u  for  tha  atatlouary  Imafa. 

Equalioo  (16)  eao  ba  ra-axpretaad  la  objact  ipaca: 

Qff(t)  ■.  ^*{4fw(u)) 

[7(-u/JV)j''"') 

-^‘{7(^u))» 

sr-HH-u/N})  ........ 9  n-HT{-u/n))  (17) 

(^^•1)  l$rm$ 
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An  MtsUaMlen  of  tkf  obovo  twulu  for  tono  pwtlcuUr  vtluot  of  N  tio  lortrueilvo.  For  lulaneo  for  ■  3  E4>(10)  roduen  to 

«.(«)- r(^) 

-r(|)  ^(f) 

.|7(|)|‘  (i») 


which  1*  ibo  powor  ipoolrum  of  r(u}  whore  the  vtrlohU  u  li  Molod  by  *  foetor  of  2.  Fbuhot  trutformiof  9t(u}  yloldt  th« 
Mom  outooer^lloo  of  the  liu|o  /(r). 

IbrNwS 


-  r(fu)  r(-f)7(-^)  (13) 

The  obovo  oquUloa  U  rcUlca  to  the  IrlpU  Mrr«Ultoa<*-*>  for  mS  fuaetlou  u 

/<”(xi,x,)  -  /(*)/{x  +  x,);(x  +  x,)*t  (20) 

404  lu  FeutUr  trooiform  7)'*(tt,v)  kaowo  u  the  bUpoctrua 

H*)(u,v)  ■  f(u)r(v)f(-tt- v)  (31) 

It  ffluit  bo  poiotod  out  thot  for  u  ■  )u,  v  ■  -^u  or  u  ■  -)u,  v  ■  }u,  which  oro  roproMoioileu  of  Uoo*  In  tho  u,v  pinno, 
oquntlon  (31)  roduco*  to  (16). 

Another  Umitlnf  com  U  when  S  -fOt,  Rememberin|  thet  7(u)  ht*  been  nono*Ueed  to  unit  et  u  ■  0  It  eon  be  loen  thoti 

5tf(u)  ■  |Jt“li  - 1//'^))^(-^)  } 

-  {7(utl  -  1/Afl)}  X  Jlrn^  {^(-^)'''‘ } 

«  r(u)  («) 

which  U  the  normoUiod  etotlonoty  Imoje  ipeetrum. 

^  .  CENTROIDINQ  PHQTQN  DATA 


Experln  eoU  were  eorried  out  uelni  elmuloted  one  dlmnulonol  photon  dote  leoeroted  by  4  computer  where  the  deelred 
low*U(ht  level  le  ochleved  by  eelictlni  the  volue  of  the  Poleeoo  moon  (N)  In  the  Poleiou  dletrlbutlon  of  photon!  delwted  In 
eoch  frome  Intervel.  in  the  portlculnr  elmulellon  experiment  deecrlbod  In  IbU  poper,  we  Mt  7  ■  3  photone/frome.  FYomee 
with  0  or  1  photnn/frome  ore  dlereivded.  beeouee  they  do  not  eorry  ooy  In/ormotion  coocernlof  the  Inteneliy  dletrlbutlon  of 
(be  Imoie,  ood  4  totol  of  $0,084  (romee  of  4  tondomly  tronelotlni  ImoKe  (of  4  blnory  itor)  contolnlni  Af  2  u  photone/frome 
were  peoeroted.  Theee  fromee  cor.  bo  (rouped  In  eele  of  fiomoe  coiitolulue  4  number  of  photone  rouilni  from  NaS  up  to  Nal3 
photone/frome  In  this  coee.  The  oetuol  dletrlbutlon  for  eoch  volue  of  N le  ehown  In  Toble  1  for  the  ^oion-dou  ueed  In  thie 

Fbr  0^  portlculnr  eel  of  N-pholon  fromee,  one  ebould  centroid  the  correepondlni  Ms  fromee,  overo(e(odd)  ood  divide 
by  N  to  And  the  normolleed  Imoge  eetimoior  Qs(t)  ood  lu  Fbu/ler  troneform  0^(u)  i 


I 


FVomee  with  N  photooi/ftome) 


(23) 


ood 


Fc 


‘ 
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(M) 


F1|un  9  dmitu  9y(r)  f«r  N«2  up  Ut  NalS  phuMiu/friuiM.  n>r  Na2  li  m»  b«  Mta  ibu  Q%[r)  ncUbtU  th*  mim  ihftp* 
u  lb«  tutoeemUtlM  of  lb«  blAtfy  but  dUtrlbuua  In  u  *nuti«r  (tbrlalMd)  rtfloa  o{  tpae*.  It  can  ilto  b<  mm  that  tb«  avtraf* 
Qiu(r)  baceoM  aara  aad  nan  ttmllar  to  tbo  itatleatiy  laafo  a*  N  incfoaM*. 

FhaM  ceotalaluf  noro  pbotoai  an  bouad  t«  tarry  otort  lalbroiatioa  coaMralat  tbo  Unan  tbaa  (raata  wltb  a  Ittttr 
Buntbtr  of  pbotoat.  Too  Umltlac  taa«  art  tboto  with  ae  laforniatloo  toatoat  (NaO)  aae  IboM  wtlb  tbo  iouft  ItMlf  (.Y  oe). 
Fbr  Uutaaot  ftamtt  wltb  Na2  oaly  carry  lafonnatlea  cooMratag  tbt  atodulut  of  all  peulbla  dltlaactt  wltala  thi  atatloaaiy 
lataft  wbortat  framat  wltb  NwS  paotoai/framt  btildti  tbt  modulua  alio  carry  laforaiatloB  coactraini  tbt  rtlativt  peiitloBi 
of  polati  la  tbt  objoet  (Ivlat  tbtrtfert  iaformatlon  about  lit  itiucturo.  Tbli  fact  allowi  oat  to  um  ftamtt  wltb  at  lout  NaS 
pbotOBi/fraiM  to  rtaoaitruel  tbt  baaft  ibrougb  tbt  ttebalqut  of  trlplt  terrtlatioa. 

Attaptlaf  tbli  Idaa  oat  ibould  flvt  lacrtatad  wtfbl  to  tboM  framtt  eoalalalaf  mort  pbotoai,  A  poMlblt  way  of  Implw 
Dtatlai  tall  u  to  eoaitrucl  aa  lina|a  ttilmator  by  a<lJta|,  wltb  taual  wtlibt.  ibt  quantltlta  0/y(r}  l.t.  : 


“  {n..Wr.UwaUUi.QN(r)} 


Nbi* 
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of  pbotoai  par  frame  N  l.t.i 
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Figure  2  -  Averaged 
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Ceniroldtd  Frames  with  N-  8  ph/tr 


Imag*  own*  ooardlMi* 


Cantroidad  Framas  with  N-  9  ph/lr 
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Cantroldad  Framaa  with  N-  12  ph/lr 
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Figure  2  -  continued  ’ 
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(c)  W) 

Figure  3  -  Compiirison  boiv'ccn  the  object , 


^MtHi  •  ^wvIlllUrO  and  the  Reconstruction. 

UiUi  N"3  ph/fr  Md  dtpudof  f»m  d(l)  oat  ((it  Ih#  phu*  »l  ?■  luiof  0i(3)i  from  p(a)  od*  foU  d(4)  uilog  ©i(6)j  from 
p(4)  oM  (tti  d(8}  uitog  9i(13)  but  from  d(8)  tlili  rtcurruM  nlklloo  con  not  b*  oppllod  onymor*  bccouM  It  would  giv«  u« 
th«  pbuo  ftt  d(ld)i  which  dooi  not  oxlit,  through  tb«  um  of  tlio  phiM  6}(24)  which  oUo  dooc  nut  oxlit. 

For  N>>4  ph/fr  tbo  rocurrinc*  rcUtion  (33)  bocomott 

©«(,*) -P(3t)- 3  d(*)  (33) 

Dtpwtlng  (lom  p(t)  end  uilng  th«  com*  proeodut*  d«ierlb«d  for  tbo  enio  of  Nb3  ph/fr  on#  flndj  d(3),  d(3)  nod  d(3)- 

In  gonortl,  for  o  froquoney  bln  k,  wboro  k  ii  n  primo  numbor,  eno  con  only  roeootlruet  lu  phuo  by  uilng  u  ovorngo 
4F/r(Afk)  luch  IhM  Nmk-fl  phetoai/frnoM.  tn  tbU  pnrtIcuUr  oxomplo  to  roeonitruet  tbo  pbnioo  up  to  bio  15  ono  ooodi 
ovorign  of  CfnaiM  with  N  up  to  13  photoni/fromo.  The  niodulut  of  T{k)  li  found  from  tbo  Fhuilor  trooifom  of  tho  luto 
corroUtlon  functiou,  tho  powor  ipoelrum  |^(v)|  • 
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a  .  OQNCLUllQNa  ANP  FUTURE  WORK 

Tb«  roiulti  ihow  tbat  eoatroldlBg,  although  a  vory  almplo  aad  eomputailoaally  ImI  toobalquo,  can  bo  uaod  to  rotrlovu  aa 
Imago  of  a  raadomly  traatlatlag  objMt  at  vory  low  light  lovol. 

^turo  work  could  addrooi  tbo  foUowlagi 
(i)  bow  to  combloo  tbo  Ib  optimum  wtyi 
(III  oxtoad  tbo  pbaoo  rotrloval  algorithm  to  two  dlmootioao  aad 
(111)  apply  tbo  tocbaiquo  to  opoeklod  Imaglog. 
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ABSTRACT 


A  comparative  study  has  been  made  of  the  Knox*Thompson  method 
and  triple  correlation  techniques  as  applied  to  image  restoration.  Both 
photon  noise  degroded  and  atmospheric  turbulence  degraded  imaging  have 
been  considered.  The  slgnal-to-noiie  ratios  of  the  methods  have  been  sni- 
died  analytically  and  with  the  aid  of  computer  simulations.  The  ability  to 
retain  diffraction  limited  information  on  imaging  through  turbulence  is 
considered  in  terms  of  phase  closure  relationships.  On  tiie  basis  of  this 
work  it  is  found  that  both  image  restoration  techniques  are  effectively 
equivalent. 
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INTRODUCTION 


The  problem  of  imaging  through  turbulent  media  hai  been  itudied  for  many 
yean.  Of  particular  intereit  here  ii  the  imaging  of  objecta  through  atmoipheric  tur> 
bulence.  A  detailed  deaeriptlon  of  thii  effect  ii  given  by  Roddier  [1]. 

At  atmoipheric  condition!  are  continuously  varying  the  ultimate  resolution 
obtainable,  from  conventional  long  exposure  imaging  using  a  large  telescope,  it 
severely  limited  by  turbulence  effeeu.  However,  if  short  exposure,  narrow  optical 
bandwidth,  images  are  recorded  atmospheric  induced  wavefront  perturbations  may  be 
'froten'.  The  turbulence-imposed  resolution  limit  is  therefore  removed  allowing 
diffraction  limited  resolution  to  be  obtained. 

It  is  well  known  that  this  type  of  imaging  of  a  point  source  object  results  in  a 
speckled  image  due  to  interference  effects.  This  is  a  consequence  of  the  complex 
amplitude  in  the  pupil  of  the  imaging  optics  consisting  of  typically  many  uncoirelated 
regions,  'seeing  cells',  the  site  scale  [2],  r^,  of  these  regions  Increulng  as  atmoipheric 
conditions  Improve.  The  long  exposure  angular  resolution  limit  imposed  is  approxi¬ 
mately,  rjX,  which  is  typically  very  much  less  than  the  difftvction  limit,  D/X,  of  a 
telescope  of  diameter,  D, 

The  technique  of  speckle  interferometry  [3]  utilises  the  diffraction  limited  infor¬ 
mation  present  in  short  exposure  images  to  obtain  the  spatial  autocorrelation  of  astro¬ 
nomical  objecu  even  in  poor  'seeing',  at  very  low  light  levels  and  in  the  presence  of 
telescope  aberrations  [4].  During  the  past  15  years,  a  large  number  of  speckle  imaging 
techniques  have  been  proposed,  see  reference  (5)  for  further  details,  although  none  are 
in  routine  use.  Two  of  the  more  promising  techniques  are  the  metnod  of  Knox  and 
Thompson  (KT)  [6], [7]  and  a  method  first  suggested  by  Weigelt  et  al  t8]-[10]  which 
uses  the  concept  of  triple  correlation  (TCJ. 

Both  techniques  are  analogous  to  speckle  interferometry  as  they  involve  averaging 
correlation  functions  of  many  short  exposure  images.  Additionally,  as  a  consequence  of 
'phase  closure'  [11],  all  three  techniques  possess  transfer  functions  which  allow 
diffraction  limited  information  to  be'  obtained  even  if  is  very  small.  In  addition  to 
the  Fourier  modulus,  the  KT  and  TC  methods  allow  direct  recovery  of  the  Fourier 
phase  of  an  object’s  intensity  distribution  from  their  associated  ensemble  averaged 
correlation  functions.  If  the  spatial  autoconeletion  of  the  object  alone  is  known  then 
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tuch  phue  information  is  potentially  only  obtainable  through  indirect  methods  [12] 
and  by  making  assumptions  concerning  the  object's  intensity  distribution.  Fov  both  the 
KT  and  TC  methods,  rigorous  analyses  (13],  [14]  show  that  an  object's  intensity  distri* 
budon  is  uniquely  related,  apart  from  bivisl  ambiguides,  to  the  relevant  averaged 
quandty  in  each  case  (although  there  is  no  guaraiitee  that  current  algorithms  recon¬ 
struct  this  unique  intensity  distribudon,  pirticularly  in  the  ^msence  of  noise). 

Although  these  coneladon  techniques  make  available  diffracdon  limited  informa- 
don,  its  quality  must  be  assessed  before  it  can  be  used  effectively.  A  common  meuure 
of  quality  used  is  the  signal  to  noise  rado,  SNR,  at  each  spatial  frequency,  U,  in  the 
Fourier  transform,  <iS(u)>,  of  the  pardcular  average  coneladon  being  considered.  The 
SNR  measure  can  be  defined  as 

SNR  = 

Oj 

where  is  d>e  standard  devladon  of  the  signal  5(u).  This  allows  similar  techniques 
to  be  compared  and  addidonally  gives  the  reladve  weight  of  the  information  piesent  so 
that  it  can  be  used  opdmally. 

A  consequence  of  using  short  exposure  dmes  and  the  low  luminosity  of  many 
utrunomical  objects  of  interest,  is  that  very  few  photon  e.ents  are  recorder  in  indivi¬ 
dual  speckle  images.  Photon  noise  is  therefore  introduced  which  can  seriously  lower 
the  obtainable  SNR  and  as  a  result  degrade  reconsmicied  image  quality.  Despite  such 
noise,  the  correladon  techiuques  have  beer  succesfu’Jy  used  to  study  low  luninosiiy 
utronomical  objects,  for  example  see  referent'  [IS]. 

The  aim  of  this  paper  is  to  compare  and  relate  the  KT  and  TC  methods  with  par¬ 
ticular  reference  to  the  SNR’s  of  their  usosciated  averaged  correlations.  Ihc  evalua¬ 
tion  of  the  SNR's  requires  the  ensemble  average  of  many  different  image  correlations 
to  be  carried  out  over  both  the  Poisson  statistics  of  the  photon  events  and  the  statistics 
of  atmospheric  turbulence  effects  Fortunately,  recent  work  [16], [17]  has  assosciated 
the  much  used  idea  of  'phase  closure*  in  radio  astronomy  with  the  comparatively  new 
ideas  of  triple  correlations  and  'phase  closure’  pupil  plane  interferometry  as  used  in  the 
optical  and  infra-red  wavebands.  This  assosciation  can  be  extended  and  generalisud  to 
allow,  previously  very  difficult,  ensemble  averages  of  many  different  image  correla¬ 
tions  to  approximated  on  the  basis  of  phase  closure  relationships. 
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The  effect  of  photon  noiie  -m  tpeckle  interferomeoy  wat  originally  invettigated 
by  Goodman  and  Belther  [IS],  their  work  being  extended  by  Dainty  and  Greenaway 
119]  with  particular  reference  to  the  problem  of  photon  biai.  The  latter  it  a  conkC* 
qusnee  of  correlating  photon  evenu  with  thenuelvet  and  it  pretent  in  both  the  KT  and 
TC  average  tignali.  Evaluation  of  the  eniemble  average  lignalt  hat  been  facilitated  by 
the  developnient  of  computationally  very  efftcient  photon  coordinate  differencing  algo¬ 
rithm!,  [20],  together  with  modem  photon  detectort,  such  at  the  PAPA  detector  [21] 
which  Ktumt  a  atream  of  time  tagged  photon  coordinate!.  Ute  of  the  time  tag  infor¬ 
mation  enablei  increated  SNR  to  be  obtained  [28].  Previout  work,  [22]-[24],  on  the 
effect  of  photon  noite  and  atmospheric  turbulein'.e  hat  been  carried  out  for  the  KT  and 
TC  techniques.  However,  only  very  simplified  aporoxlmationt  to  ensemble  averages  of 
image  correlations  over  the  photon  noite  and  atmctpheric  statistics  have  been  obtained. 
Previout  studies  have  alto  only  been  concerned  w'th  very  specific  regions  in  Fourier 
kpace  of  the  ensemble  averaged  signals,  the  worV  on  the  TC  beln"  only  concerned 
with  regions  having  the  lowest  SNR  which  are  intact  not  used  in  practice  [2S].  The 
work  on  the  Kf  method  hu  only  looked  at  the  basic  two  dimensional  conflation  tech¬ 
nique  u  origii.u  iy  proposed  [6],  whereas  the  technique  can  be  completely  generalised 
to  give  a  four  dimensional  average  signal  that  contains  considerably  more  information 
and  it  analogous  to  the  four  dimensional  TC  signal. 

The  pretent  work  it  organised  u  follows.  In  section  2  the  mathematical  form  of 
the  TC,  the  KT  correlation  and  the  autoconelation  arc  defined  in  both  image  and 
Fourier  space.  A  formalism  for  the  TC  it  put  forward  which  allows  it  to  be  simply 
related,  in  image  space,  to  the  KT  correlation.  The  correlation  techniques  are  i.een 
diagramatically  and  matherrutically  to  be  intimately  related,  each  being  simply  correla¬ 
tions  of  images  with  themselves  weighted  by  complex  exponentials.  Additionally,  a 
further  technique  [26], [27],  known  as  the  phase  gradient  method,  is  included  which  is 
shown  to  be  a  sub-set  of  the  more  general  KT. 

In  section  3  the  implementation  of  the  correlation  techniques  at  low  light  level  it 
described.  Of  particular  importance  is  the  development  of  a  differencing  algorithm 
wliich  allows  individual  tecdont  or  sub-planes  of  the  Fourier  transform  of  the  TC, 
known  at  the  bispcctrum,  to  be  evaluated,  the  full  TC  differencing  algorithm  being 
reduced  to  a  weighted  KT  Implementation.  This  hat  the  advantage  that  only  the  higher 
SNR  regions  need  be  calculated  at  opposed  to  the  whole  of  the  bitpectrum  which  is 
evaluated  if  a  full  differencing  algorithm  is  implemented.  This  problem  can  be 


52 


-5- 


overcome  by  evalutting  the  required  bitpeetrum  planet  directly  in  Fourier  ipano.  Thii 
however  lequiret  each  frame  to  be  Fourier  tranifomed  which  becomet  lets  efficient 
relative  to  photon  differencing  at  low  light  leveli. 

The  effect  of  photon  noiie  on  SNR't  it  conildered  in  lection  4.  An  ettimate  for 
the  error  on  the  phaie  of  a  complex  lignal  it  propoied  and  iti  dependency  on  photon 
nolle  itudled  for  the  KT  and  TC  lignali.  The  phue  error  ettimate  it  ihorvn  to  be 
related  to  the  SNR.  Averaging  varioui  image  correlationi  over  the  photon  ataditlei  hat 
been  facilitated  by  a  computer  implemenution  of  a  generaUtadon  of  the  method  of 
Goodman  and  Beliher  [18].  It  hat  therefore  been  potiible  to  extend  the  work  on  pho¬ 
ton  nolte  dependency. 

Section  5  briefly  contidert  the  imaging  of  a  randomly  iranilatlng  object.  The 
necetiity  to  centroid  framei  before  carrying  out  the  KT  procettlng  it  explained  and 
alto  the  problemt  aiiociated  with  centroiding  imagei  in  the  pretence  of  photon  noiie. 

In  lecdon  6  the  ute  of  the  KT  and  TC  for  imaging  through  atn>otpheiic  tur¬ 
bulence  it  itudied.  The  diffbaedon  limited  Informadon  preient  in  the  tranifer  funedont 
of  tpnckle  interferometry,  the  KT  and  the  TC  it  teen  to  be  a  direct  coniequenee  of  the 
concept  of  phaie  cloture.  Extenilve  ute  hat  been  made  of  a  computer  limularion  of 
the  imaging  of  objecu  dirough  atmoipherlc  turbulence  uiing  a  large  teleieope.  The 
dependency  of  SNR't  end  tranifer  funedoni  on  mean  number  of  photon  eventi  per 
image,  17,  and  mean  number  of  ipecklet  per  image,  K,  it  itudied  analydcally  and  via 
Monte-Carlo  type  computer  limuladont. 

Finally,  in  lecdon  7  reiulti  are  compared  and  concluiioni  drawn  at  to  the  reladve 
meritt  of  the  ipeckle  imaging  technlquei  with  reipect  to  the  treat  itudied. 
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2.  Definitions  and  relationships 


Coailder  t  two  dimenilon&l  intenilty  diitiibution,  i(x),  and  iU  Fourier  tnmifom, 
/(u),  defined  by 

/(u)  ■  J  i(x)  txp(-ivjux)  dx  .  (1) 

Both  the  triple  correlation  (TC),  t8]*(10],  and  the  correlation  of  Knox  and 
Thompaon  (KT),  [6][71,  may  be  defined  in  imaje  apace  at  eorrelationt  of  i(x),  or  in 
Fourier  apace  u  producti  of  /(u).  The  TC  and  iu  Fourier  trantform,  the  biapectrum, 
are  defined  by 

■  J  <*(*)  K*+Ki)  f(K+*2)  dx  (21 

;™(ui,ua)  ■  /(u,)  7*(ui+uj)  /(uj) . 

In  an  analogout  way  the  KT  double  correlation  and  Itt  Fourier 
by 

<*^(Xj,Au)  ■  I  /*(x)  <(x+xj)  «p(2rt/Aiix>  dx 

/*^(ui^u)  -  /(Ui)  / (Ui+Au) .  (5) 

where  X;-Xj,+Xi,  and  Xj^Xjjj+Xaj,  are  two  dlmcniional  tpatiel  coordinate  vectort, 
Uj«U)^<-Uij„  U2*ni2,-)-U2;  and  Au^Au^i-t-AUj  are  two  dimenaional  tpatial  frequency  vec- 
ton,  and  the  tupencript  *  indlcatec  the  complex  conjugate.  The  x  and  y  tubteriptt 
denote  orthogonal  component  vectors. 

It  in  apparent  that  Eqt.(l)>(5)  are  four  dimenaional.  Thlt  fact  cautet  Implemenu- 
tion  problemt,  extentive  evaluation  time  and  data  atorage  requirementa,  particularly  if 
the  correladont  are  performed  utlng  digitized  imaget  on  a  computer.  The  buic  KT 
technique  doet,  however,  only  contitt  of  evaluating  the  three  aub-planea  in  Fourier 
apace  coireapondlng  to  Au^Aiijp  AiiaiAu^,  and  AuBAiiy-fAu^,  with  the  Au  vector  com- 
ponenta  having  conatant  valuea.  If  digitized  imagea  are  uaed  then  Au  normally  [7] 
corteaponda  to  the  fundamental  aampling  vector  interval.  Further,  different  valuea  of 
Au  may  be  uaed  and  a  number  of  aub-planea  can  be  evaluated.  Thia  ia  analogoua  to  the 


(3) 

.  maform  are  defined 

(4) 


3i 


-7- 


many  iub*planer  cotitained  in  the  four  dimenilonal  biipectrum,  we  Fig.l. 

The  relationihipi  between  the  two  methodi  become  more  apparent  If  Eq.(2)  ii 
rewritten,  making  the  lubatitution  uj  •>  Au.  A  single  plane  of  the  bispecmim  can 
therefore  be  described  by 

/^(Ui,Au)  ■  /(ui)  /*(u,+Au)  /(Au) .  (6) 

The  coireiponding  image  space  expression  is  obtained  by  inverse  Fourier  transforming 
yielding  the  double  correlation 

i^(Xi,Au)  ■  /(Au)  J  /*(*)  /(x+Xi)  dx 

■  /(Au) /*’'(Xi,Au)  .  (7) 


On  comparing  equations  (5)  and  (6)  it  can  be  seen  that  for  a  fixed  spatiil  fre¬ 
quency  difference,  Au,  the  defined  blspertrum  plane  corresponds  to  a  weighted  version 
of  the  KT  product  in  Fourier  space.  Ttie  complex  weighting  factor  is,  however,  signal 
dependent  and  as  such  posseives  a  number  of  important  properties  which  will  become 
evident  in  later  sections. 

In  order  to  understand  the  correlation  techniques  it  is  instructive  to  consider  the 
autocorrelation  process.  This  conresponds  to  the  double  correlation  represented  by 
Eq.(4)  with  the  substitution  Au  ■  0.  The  resulting  representations  in  the  image  and 
Fourier  domains  are 


(-^(x,)  -  /'^(x,.  0)  -  J  i*(x)  ((x+x,)  dx 


(8) 


A'(Ui)  ■  /^(Ui,  0)  “  /(u,)/*(u,)  . 


(9) 


The  autocorrelation  of  i(x)  is  the  cotrelation  of  i(x)  and  /(x)  multiplied  by  a  complex 
exponential  factor  with  aero  spatial  frequency,  see  Flg.2(a).  Unfortunately  the  complete 
spatial  symmetry  of  the  operstion  prevents  the  preieivation  of  Fourier  phase  informa¬ 
tion.  The  argument  oi  Eq.(9)  is  always  aero. 


0^(11,)  -  (>(Ui)  -  i)>{u,)  ■  0  , 


(10) 


where  (|)(u)  it  the  phase  of  /(u).  and  arg 


defines  the  phase  of  the  complex  numbei 
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modulo  2r  . 


Now  coniider  the  KT  correlation  of  Bq.(4).  Thli  deflne»  the  correlation  of  l(x) 
end  <(x)  multiplied  by  a  complex  exponential  factor  with  ipatlal  frequency  greater 
than  aero,  aee  Pi|.2(b).  For  lamptad  itnagea  thii  generally  correiponda  to  the  funda¬ 
mental  ipatial  ftequeney.  The  KT  correlation  enablea  Fourier  phue  difference  informa¬ 
tion  10  be  retained.  The  argument  of  Eq.(5), 

af||/*^(Ui,Au)'i ■  0*^(Ui,Au)  -  Wui)  -  $(ui+Au) ,  (1 1) 


glvei  the  phau  difference  between  two  apadal  frequenelet  aeparated  by  the  vector  Au. 

The  TC  of  Biq.(7)  almllatly  allowa  the  retention  of  phaie  difference  Information. 
However,  abaolutk'  apadal  poiition  it  indeterminable  due  to  the  property  of  the  com¬ 
plex  weighting  flaetbr  that  linear  phue  li  cancelled.  Additionally,  biipectrum  planet 
contain  phue  infonnab’w  encoded  by  correlating  i(x)  with  l(x)  multiplied  by  complex 
axponentlalt  with  higher  xpatial  frequMclet,  for  an  example  tee  Fig.2(c).  The  argu¬ 
ment  of  Bq.(3), 

«rrf/^(ui,ua)|  ■  0'®(Uj,U3)  •  ♦(uj)  -  <>(U|+Ua)  +  0(uj) ,  (12) 


allowt  the  phue  dlffcrtnce  between  two  tpatlal  frequencict  teparated  by  the  vector  Au 
to  be  determined  provided  the  phaie  at  frequency  Uj  ■  Au  it  already  known.  In  par- 
deular,  for  U} «  Au  and  Au  equal  to  th«  fttndamenul  lampling  frequency  in  the  cate 
of  dlicrete  dau,  the  TC  givei  identical  phue  difference  information  at  the  KT  apart 
from  the  phue  at  Au  .  Thii  it  the  llnerr  phew  which  cancel!  with  the  linear  phase  in 
the  KT  phaw  difference  to  produce  the  shift  invariant  property  of  the  bitpectrum, 

®^(Ui,Au)  ■  0*^(Ui4u)  +  «Au) .  (13) 


Another  procetilng  technique,  the  phaw-gradient  method  of  Aitken  et  al, 
[261[27],  it  very  closely  related  to  the  correlationi  Juit  dweribed.  Phase  information  it 
again  reutlned  but  in  the  form  of  tu  danvxtive,  not  explicit  phue  dlfferencei.  The 
phaw  derivative  information  it  contained  in  the  Imaginary  part  of  a  double  product  in 
Fourier  space, 


/^fui.e)-!  IMAO 


AuilAu,) 


X 
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-  (  1  /(U,)  p  lt>’(u,)  )  f  . 

n 


(14) 


where  IMA 


Indleatei  the 


Imtginery  part,  ^'(u)  it  the  phaie  derivative,  /(u)  it  the 


derivative  of  /(u)  and  a  ii  a  vector  gradient  in  image  apace.  The  reiult  can  be 
achieved  through  image  apace  in  a  number  of  waya.  For  example,  the  following  two 
aimple  correlationa  may  be  performed, 


l^(Xi,  +a)  ■  J  /‘(X)  i(x+xi)  t+a.x  +  b)  dx  and 


-0)  ■  J  (*(x)  l(x+xi)  [-a.x  +  b]  dx, 


(15) 


where  b  la  a  real  conatant,  and  then  combined  and  Fourier  tranaformed  to  yield  the 
reault  of  Be[.(14).  Thia  la  aimilar  to  the  optical  implementation  of  the  technique  uaed  by 
Aitken  et  al, 


;^°(u,,a)-  f7|l'*®(xi, -a)  ■  -  FT|i^°(x,. +a)  J 


•  (  [/(ui)r(ui)]’-/(ui)/*(u,)  )~ 


■  (  IMAG 


-  (  I /(ui)  I*  V(u,)  ) -^  ,  (16) 

Tv 

The  expreaaiona  of  Eq.(15)  correipond  to  correlationa  of  ((Xj)  with  /(Xj)  multiplied  by 
linear  ramp  fUncUona  of  oppoilte  gradient,  aee  Fig.2(d). 

Thia  technique  la  analogoua  to  the  KT  with  a  Au  value  auch  that  there  it  only  a 
linear  fluctuation  of  phate  over  that  dittince.  Thia  cotretponda  to  the  complex 
exponential  in  the  KT  correlation  of  Bq.(4)  having  a  apttial  frequency,  Au,  auch  that 
the  image  being  correlated  only  extendt  over  the  linear  region  of  the  exponential.  The 
imaginary  part  of  the  KT  product  in  Fourier  apace  therefore  giver  the  phate  difference 
between  tpatial  frequenciet  Au  apart.  Obviously  division  by  Au  now  yields,  apart 
frx>m  a  constant  factor,  the  weighted  phase  gradient  conuined  in  the  result  of  Eq.(14). 
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IMAO  j/CUi)  /*(Ui) 


Au 


/MA<W/”‘(u,Au) 


>1 


Jau^ 


i-/WAO  j/(u,)/*(Ui+Au)i 

I  Jau-»o 

/MA{j|;(u,)  [  ;(ui)  +  /(u,)  Au  ]*| 


Au 

s:  I  /(Ut)  I*  V(Ui)  . 


(17) 
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3.  Implementatton  of  correlation  techniques  at  low  light  level 


In  s  number  of  litustioni,  for  eximpte  'ftmtng  out’  turbulence  effectt  on  Imtg- 
ing  through  the  etmotphere,  it  li  dMireble  to  record  miny  ihort  exposure  Images.  At  a 
result,  very  ftw  photon  events  an  deteeted  in  Individual  images.  However,  it  can  be 
shown  that  apart  from  nmovsble  photon  biu  terms,  one  of  the  properdet  common  to 
all  the  correlation  techniques  it  that  the  ensemble  average  over  many  low  light  level 
frames  it  equivalent  to  tlie  high  light  level  oonelatlon.  At  such  low  photon  levels  it  it 
more  efficient  computationally  to  perform  the  correlations  in  image  space  via  photon 
diffenneing  algorithms  rather  than  in  Fourier  space  via  fast  Pouiisr  transform  routines 
and  multiple  products. 


Using  the  notation  of  Dainty  and  Greenaway,  [19],  the  pth  photon  limited  image 
ia  npnsented  by  a  sum  of  Dirac  delta  functions  positioned  at  the  coordinates  of  the  Np 
photon  events, 


(18) 


snd  its  Fourier  transform 

DpW  -  k)  exp(“2t0ix),  (19) 

lm\ 

where  Xpi,  is  the  position  vector  of  the  kth  photon  in  the  pth  frame. 

The  photon  differencing  algorithm  tor  tlie  KT  can  be  written  is,  [20], 

i  HXpi^  -  X,*,  +  X,)  exp(2JVAuXp*,).  (20) 


Th'Ji  the  KT  correlation  for  each  image  it  carried  out  by  placing  the  complex  number 
txp(2ifJAmpi^)  at  position  x  ■  Xp^-Xpp^  whenever  this  vector  coordinate  difference 
occurs  between  two  photon  events.  The  tame  algorithm  can  be  used  to  obtain  the  auto¬ 
correlation  but  with  the  substitution  Au«0. 

There  are  many  approaches  to  obtaining  the  TC,  [28],  but  a  differencing  algo¬ 
rithm  bated  on  Eq.(7)  is  most  appropriately  tuted  here.  That  it,  each  bitpectrum  plane 
corresponding  to  a  different  Au  value  it  evaluated  indepeildently  using  the  following 
expression 
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i^(«i,Au)  ■  exp(“2:VAux^*) 

N  N 

■  U  8(*;rt,  -  +  *l)  #Xp(-2jVAUXp*,).  (21) 

Thus  th«  eoirelttion  It  reduced  to  the  limpter  KT  differencing  algorithm  followed  by  e 
tingle  complex  muldplicition  end  hence  lequirei  very  little  tdi  tionel  compuution 
time. 

‘Dte  phtie  gradient  conelationi  of  equation  (IS)  can  be  evaluated  using 
/'^(Xidt)  •  titi  "  Vi 

Jt|aU|pl 


Photon  blM  termi 

The  eniemble  avenging  of  photon  limited  correlttioni  unfortunately  introduces 
photon  biu  termi.  Such  teiini  cotreipond  to  tlie  conelation  of  photon  eventi  with 
themielvei  and  do  not  contribute  any  uiefbi  information  to  the  avenge.  The  elimina¬ 
tion  of  biai  terma  ii  limply  eairied  out  by  not  Including  in  the  dlffetencing  elgorithmi 
any  diffeiencei  between  a  photon  event  and  itself.  Therefore  Eqi.(20)  and  (22)  should 
include  the  proviso  that  itiedEj  and  equation  (21)  is  modified  tc  nad 

/^(Xi,Au)»  £  fl(xi-x^  exp(-27^Aux^ 

I  6(V.  ■  V,  +  *i  )  exp(-2n/Aux,*,)  (23) 

It  is  important  to  note,  ftom  the  standpoint  of  photon  noise  analysis,  that  the  bias  tenns 
of  eech  frame  [19],  are  removed  during  processing  and  not  the  ensemble  average  bias 
teimi. 

The  unbiased  quantities  averaged  in  each  technique  may  be  stated  in  Fourier 
space  as  follows. 

(a)  The  autocorrelation  : 

Dp(u,)  !3;(ui)  -  D^(0).  ■  (24) 

(b)  The  Knox-Thompson 
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/),(Ui)  n;(u,+Au)  -  o;(Au).  (23) 

fo)  The  triple  coireUtion  for  each  blipectrum  plane  ujaAu 

^p(Ui)  i>;(Uj+Au)  £)p(Au)  -  |D,(u,)|*  -  |Z)p(ui+Au)|*  -  |Z)p(Au)|»  +21Vp.  (26) 
(d)  The  phue  gradient  correlation 

[i5p(Ui)  D^(Ui-t-Au)  -  /)*(Au)  -  2)*(u,)  D/ii,+Au)  +  Z)^(Au)  ].  (27) 

The  variance  of  the  above  quantities  with  the  number  and  distribution  of  photon 
events  in  each  image  describes  the  effect  of  photon  noise  on  such  conelation  tech* 
niquei. 
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4.  Photon  noUe 


With  the  exception  of  the  autocorrelation  the  aim  of  the  teehniquei  under  etudy  ii 
primarily  to  preierve  aome  eidmate  of  the  Fourier  phaie.  Ai  the  averagea  are  com¬ 
plex  number!  both  the  variance  of  the  teal,  an\  and  imaginary  part,  of  the  Fourier 
tranifoim  are  anidled,  together  with  the  eovaiiuice  of  the  real  and  imaginary  parti, 
Cov(JJi).  Prevlouily,  [22], [23],  the  lignal  to  nolle  ratio,  SNR,  meaiure  uied 
correipondi  to, 


SNR..  — iJMSUUStif—  M 

^  ktandard  deviation  of  tlgnal 

.  VKI  (28) 

when  +  O/*,  and  M  ii  the  number  of  independent  naliiatloni  uied.  Thli 

meaiure  providei  a  naionable  indication  of  the  SNR;  however,  of  enncerr  here  li  the 
error  on  the  phue  of  the  complex  tlgnal.  An  eitimaie  of  the  phaie  error,  0b,  may  be 
bo  found  by  eoniidering  the  variance  of  the  lignal  in  a  direction  perpendicular  to  the 
mean  ilgnal. 

^  Vo;*  coi*  d  +  tin*  ♦  -  Cov(JJi)  lin  2d  \ 

tan  Oh  ■  - . . — .  (29) 

I  <S>  I 


when  <1  ii  the  argument  of  the  mean  complex  lignal,  <S>.  Provided  the  number  of 
realiiationi,  M,  uied  to  obtain  the  mean  ilgnai  li  large,  dien  then  the  umall  angle 
approximation  of  the  tangent  may  be  uwd  to  yield  the  phaie  enor  eitlmate  directly. 


0,-- 


'hj^  coi'^  ^  *  iltt^  ^  -  Cov(IJi)  iln  20  I 
I  <5>  I 


(30) 


The  o/,  Ob’  and  the  covariance,  Cov(lJ{),  of  all  the  unblued  eitLmated  quanti- 
tiee,  Eqi.(24)-(27).  have  been  evaluated  where  appropriate,  lee  Appendix(l).  The 
method  of  evaluation  ii  e  computer  Implementation  of  a  generaliiadon  of  the  method 
deicribed  by  Goodman  and  Beliher,  (18],  for  the  power  ipectrum.  The  evaluated  vari- 
oncei  aie  now  uied  to  compare  0,0  effect  of  photon  noiie  on  the  correlation  tech- 
nlquei. 
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The  effect  of  photon  noise  on  power  spectrum  estimation  has  been  adequately 
studied  by  other  authors,  [18],  [20],  Of  particular  interest  in  this  work  is  the  com¬ 
parison  of  photon  noise  sensitivity  of  phase  difference  estimates  obtained  using  the  KT 
and  the  TC  techniques.  It  is  apparent  from  the  relationships  described  in  section(2)  that 
the  KT  and  TC  sub-planes  may  be  directly  compared  as  they  effectively  yield  the  same 
Information.  The  additional  term  present  in  each  bispectrum  sub-plane  is  in  general 
known  from  sub-planes  defined  by  a  smaller  magnitude  frequency  vector,  uj.  How¬ 
ever,  the  bispectrum,  being  an  average  of  a  triple  product,  contains  noise  terms  to  the 
third  power  causing  the  variance  of  the  observed  signal  to  be  greater  than  that  of  the 
KT  double  product.  This  is  a  basic  property  of  high  order  correlkdons  and  suggests 
the  use  of  cotrelations  of  the  lowest  order  possible  which  allow  recovery  of  the  tame 
informadon.  For  example,  Plg.3(a)  shows  the  (heotedeal  error,  0b,  for  a  tingle  image, 
on  the  phase  difference  etdmatet  resulting  from  using  the  KT  plane  (UjyaO,  Au^bO) 
and  the  blepectrum  plane  (u^fOAUx,  Uiy«ujyB0).  The  results  are  for  an  'asteroid'  type 
object,  tee  Fig.3(e),  which  it  approximately  10  x  20  pixels  in  a  sampled  region  of  128 
X  128  pixels,  and  an  imaging  system  with  a  transfer  funcdon  of  value  unity  for  all  tpa- 
dal  frequenciei.  The  mean  number  of  photons  per  frame,  Ff,  it  10.  Crott-tecdonc 
through  these  planet  are  shown  in  the  graph  of  Flg.4(a),  corresponding  to  Au^bUjx^, 
1.0.  4  timet  the  frequency  sampling  interval.  The  signal  to  noise  ratio,  SNR„,  for  the 
tame  two  planet,  are  shown  in  Flg.3(b).  If  cross-secdont  through  these  plots  are  com¬ 
pared,  for  example  Fig.4(b)  shows  the  SNR  of  the  KT  section  AUj^  and  AUy,  Ui,b0) 
and  the  bispectrum  (u2yB4  and  U2yBU]yB0,  a  further  indication  of  the  advantage  of  the 
double  product  estimation  it  obtained. 

As  a  second  example,  consider  estimating  the  mean  number  of  photons  per  frame 
using  the  bispectrum  element  uibu2b0  and  the  KT  element  UjbAubO.  The  SNR  for  a 
tingle  frame  for  the  KT  is 

and  for  the  bispectrum  it 


It  it  apparent  that  the  KT  hat  a  greater  SNR  then  the  blspecuum  at  this  point  for  all 
values  of  f^.  In  particular,  for  the  limiting  cates  of  ^approaching  infinity  and  f7  very 
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much  leii  than  unity  the  following  retulti  are  obtained  uaing  the  expreiiions  for  the 
lignal  vaiiancei  given  in  Appendix(l]. 


Urn  SNRkj' 


2 


(33) 


and 


and 


Urn  SNRtc  ■ 

3 


SNRjf-^ 


/7«  1 


/7«  1 


(34) 

(35) 


(36) 


The  plots  of  9e  tnd  SNR„  shown  in  Figs.3  and  4  Indicate  a  close  relationship 
between  these  quantities  at  expected  from  Bqs.(28)  and  (29).  Typically, 


0B  ® 


1 

<^NR„' 


(37) 


however,  the  only  generalisation  that  can  be  made  provides  an  upper  bound  for  the 
phase  error  in  relation  to  SNR„. 


(38) 


Thi;:  inequality  assumes  nothing  about  the  statistics  of  the  fluctuating  signal. 
Nevertheless  it  is  a  useful  expression  at  it  allows  a  bound  on  the  phase  error  to  be 
defined  without  calculating  @£  which,  in  general,  it  a  far  more  difflculi  task  than 
evaluation  of  SNR„.  The  additional  VI  arises  in  this  expression  due  to  potential  corre¬ 
lation  between  the  real  and  imaginary  parts  of  the  complex  signal,  resulting  in  a  non¬ 
zero  covariance. 

Nltenton  and  Papalloliot,  [22],  assume  that  in  the  KT  method,  for  small  Au 
values,  o,„  is  the  modulus  of  a  complex  error  vector  which  can  possets  any  phase 
between  ±rt  with  equal  probability.  This  then  gives,  on  averaging  over  all  possible 
phases,  in  the  expression  for  the  phase  error  given  by  Eq.(37).  If  the  simple  imaging 
esse  previously  described  it  again  considered  then  the  validity  of  these  assumptions 
can  be  investigated.  Figure  3  shows  plots  of  the  ratio  of  Og,  defined  by  Eq.(30),  to 
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for  the  tame  KT  and  biipectrum  pianos  u  used  in  Fig.3.  Both  plots  show 
the  tame  important  features  and  to  only  one  of  them,  the  KT  plane  in  this  cate,  needs 
to  be  studied  in  detail.  It  should  be  noted  that  the  result  it  unity  in  general  with,  in 
particular,  the  exception  of  the  region  around  Au^^.  This  it  important  u  this  it  the 
main  region  used  for  the  reconstruction  of  the  object’s  Fourier  transform.  These  results 
are  straightforward  to  explain  if  the  Fourier  transform  of  a  tingle  realiution  it 
r  ^presented  by  D^(u),  defined  by  Eq.(19).  This  it  a  summation  over  a  number  of  unit 
Ttoduli  complex  exponendalt  resulting  in  a  complex  exponential  at  each  ipatial  fre¬ 
quency.  u,  which  pctsetset  a  phase  that  it  effectively  uniformly  distributed  modulo 
2n,  conflated  over  a  short  range  of  spatial  frequencies  and  uncoirelated  between  reali¬ 
sations,  the  corrrladon  over  spatial  fiequenciet  increuing  if  frames  contain  very  few 
photon  events.  On  this  basis  it  is  easily  teen  that  away  from  the  axes,  where  the  KT 
signal  it  a  product  of  Fourier  transform  values  for  well  separated  spatial  frequencies 
(large  Au),  then  the  real  and  imaginary  parts  of  the  signal  an  essentially  uncorrelated 
random  variables  with  equal  variances.  Thit  gives  rite  to  the  tame  result  at  that  of 
Nitenton  and  Papalioliot.  However,  the  regions  close  to  the  axis,  (small  au) 
correspond  to  where  the  Nitenton  and  Papalioliot  attumptiont  are  made  whereas  it  it 
exactly  ''.'ese  regions  for  which  the  utumptiont  break  down.  The  phue  of  the  signal 
in  each  frame  near  the  axis  tends  to  be  small,  tee  Eq.(ll),  giving  rite  to  a  large  real 
and  small  imaginary  part.  The  real  and  Imaginary  parts  therefore  tend  to  be  correlited 
and  a  non-zero  covariance  results.  Additionally  the  real  variance  is  large  and  th  <  .1;'- 
ginary  verlance  small.  In  general,  the  phase,  of  the  averaged  signal  it  also  tm."  .id 
so  the  following  expression  for  the  phase  enor  resulu. 


_ 1_ 

i<S>|  ^ 


(39) 


As  can  be  seen  in  Fig.5(b),  this  is  smaller  than  that  predicted  by  V{'J2SNR„)  which 
contains  the  Irjrge  real  variance  contribution.  Unfortunately,  the  nature  of  sampling 
leads  to  occatssional  discontinuities  in  the  Fourier  phase  of  the  object  being  imaged. 
This  can  consequenUy  give  a  large  signal  phase,  0,  which  increases  the  contributions 
of  the  real  variance  and  the  covariance  to  the  phase  error,  see  Eq.(30).  The  result  is 
that  Eq.(38)  tends  to  the  equality.  This  also  can  be  seen  in  Fig.5(h)  where  large  values 
for  the  plotted  rsiio  .ire  visible  near  the  axis.  7'he  implication  of  this  is  that  a  larger 
number  of  frames  are  required  in  these  high  phase  error  regions  than  that  predicted  by 
the  SNR„  value.  Ihe  number  of  frames  required  could  in  fact  be  more  than  double  the 
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don  of  the  expected  error  on  the  phase  of  the  mean  signal,  again  indicating  the  useful 
neat  of  evaluating  the  SNR^  as  opposed  to  the  more  complicated  Sg. 
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5.  Imaging  a  randomly  translating  object 


The  KT  Imaging  technique  it,  unlike  the  TC,  dependent  on  image  ihifU  teaultiiig 
in  lignal  degradation.  Consider  trying  to  image,  using  the  KT  technique,  an  object 
whose  image  is  raitdomly  translating.  Such  translations  may,  for  example,  be  a  conse¬ 
quence  of  the  first  order  degrading  effect  of  atmospheric  turbulence.  The  cause  of 
inr,age  shift  is  wavefront  tilt  relative  to  the  optical  axis  of  the  imaging  system,  which, 
as  in  die  case  of  atmospheric  turbulence,  is  cffoctively  independent  of  other  wavefront 
perturbations.  Suppose  such  tilts  result  in  a  Gaussian  probability  distribution,  p(a) ,  of 
image  shifts,  a. 


where  is  the  standard  deviation  of  the  shifts.  An  ensemble  averaged  transfer  func¬ 
tion  may  therefore  be  defined  which  reduces  the  KT  signal  by  an  amount  dependent  on 
Au  and  9^  . 

<T/u)>  ■  Jexp{27^ua)  exp(2j^(u-i-Ati)a)  p(a)  da  (41) 

■  exp(-27t*Au  V) 

As  an  example  consider  a  frequency  sampling  interval  Au  of  0.5  arcsecondt~‘  and  a 
standard  deviation  of  image  shift  of  0.5  arcseconds.  Such  numbers  are  reasonably  typi¬ 
cal  in  the  astronomical  case  and  result  in  the  signal  being  reduced  by  approximately 
0.75.  This  effectively  reduces  the  improved  SNR  obtained  by  averaging  the  KT  dou¬ 
ble  product  as  opposed  to  the  equivalent  bispectrum  plane. 

In  the  absence  of  photon  noise  centroiding  images  before  correlating  makes  the 
KT  method  a  shift  invariant  process  and  so  removes  the  above  effect.  Assuming  sam¬ 
pled  data  and  Au  equal  to  the  spatial  frequency  sampling  interval.  The  Fourier  spec¬ 
trum  of  the  centrolded  image,  /(u),  can  be  expressed  using  Fourier  shift  theorem  as 


/(u)  ■  f(u 


au 


tl/(Au)l, 

and  the  KT  average  performed  with  centroided  frames  becomes 


(42) 


/(u)  Tfu+Au) 


'  '1,TOU>|J 


u  ..  u+Au 

|/{Au)| 
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(43) 


The  tdvintage  of  performing  a  double  product  for  the  KT  hat  been  Icit  and  the  corre¬ 
lation  hu  now  become  a  triple  product  In  Fourier  ipace,  equivalent  with  reipect  to 
third  order  nolle  termi  to  the  bitpeotrum. 


At  low  light  leveli  the  problem  li  complicated  by  the  pretence  of  photon  nolle 
which  prevents  the  accurate  determination  of  the  image  centroid.  The  KT  average 
obtained  on  averaging  a  leriet  of  low  light  level  randomly  moving  imaget,  centrolding 
each  image  before  proceiiing,  can  be  ihown  to  give,  aee  Appendix  (2). 


<  %  (Wp-1 )  /(u+~)  >  (44) 

where  Np  ii  the  number  of  photon  event!  in  the  pth  image  and  /(u)  it  the  Fourier 
traniform  of  the  object  notmallied  by  HO).  The  eniemble  average  over  the  Poliion 
itaditici  of  the  number  of  photon  eventi  in  each  image  ii  left  unreiolved  due  to  the 
intriniic  object  dependency  of  the  above  reiult.  Thli  reiult  differ*  from  the  deiired 
reiult, 


/(u)  /“(u+Au)  >,  (43) 

became  of  the  pretence  of  the  fiequency  ihlfl  vector,  The  latter  li,  however, 

negligibly  imall  for  all  imaget  except  thoie  containing  a  few  photoni.  In  particular,  for 
framei  containing  two  photont  the  autocorrelation  reiulti  and  no  phaie  information  it 
preierved.  The  centrolding  error  effect  obviouily  increaaei  with  Au.  Thii  remit  it 
important  at  it  luggesti  that  the  technique  cannot  be  uied  in  iti  preient  form  when 
only  a  few  photon  evenu  per  frame  are  recorded.  Thlt  it  in  contratt  to  the  ihlft  invari¬ 
ant  TC  which  allows  objecti  to  be  Imaged  reasonably  well  at  very  low  photon  levelt 
[29].[301. 
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6.  Imaging  through  atmospheric  turbulence 


Both  the  KT  vid  the  TC  techniques  have  been  used  luecesafully  for  imaging 
Mtronomical  objecu.  Their  lucceit  it  bated  upon  the  technique!  overcoming  the 
degrading  propertiei  of  atmoiplieric  turbulence.  Many  images  are  recorded,  with  expo¬ 
sure  times  short  enough  to  fteeie  atmospheric  perturbatlont,  and  processed  to  yield 
ensemble  average  transfer  functions  which  contain  diffiractlon  limited  information. 
Obtainable  retoiutlont  are  therefore  iildmately  limited  by  telescope  eharaoterlttlct. 

The  complex  mathematical  nature  of  the  transfer  functions  concerned  mtket  a 
complete  analytical  study  extremely  dlfAeult.  We  have  therefore  used  Monte  Carlo 
computer  simulations  to  obtain  a  reaionabte  idea  of  the  effects  of  the  sutistical  fluctua¬ 
tions  of  atmorpheric  turbulence.  The  wavefront  phue  penurbadont  induced  by  the  tur¬ 
bulence  have  been  considered  to  have  a  Qausiian  comladon  fUncdon  with  a  standard 
deviadon  of  value  being  the  Pried  [2]  parameter  describing  the  atmospheric 

condidons  being  studied.  Wavefront  amplitude  perturbadons  an  not  pnsent  in  the 
slmuladons.  The  simuladons  allow  any  form  of  telescope  pupil  function  to  be  used  and 
the  imaging  of  arbitrary  object  intensity  distribudons.  The  introducdon  of  photon  noise 
is  also  possible. 

Inldally,  consider  a  flmher  slmpiiflcadon  of  the  problem  bued  on  the  work  of 
Roddier  [16].  Suppose  the  atmosphem  is  modelled  by  a  dlsttlbudon  of  'seeing'  cells 
of  diameter  in  angular  spadal  fnqueucy  space.  The  wavefront  phase  is  usumed 
constant  across  each  seeing  cell  and  the  complex  amplitude  is  assumed  uncorrelated 
between  cells.  Using  this  basic  model  die  processing  techniques  may  be  simply 
explained. 

Consider  two  'seeing'  cells  separated  by  a  vector  in  the  teletcope  pupil,  Xu, 
where  X  is  the  mean  wavelength  and  u  an  angular  spadal  frequency  vector.  If  a  point 
source  is  imaged  through  the  telescope  using  a  pupil  funedon  conslsdng  of  two  aper¬ 
tures,  corresponding  to  the  two  'seeing'  cells,  then  a  fringe  pattern  is  produced  with  a 
narrow  spadal  frequency  bandwidth.  The  nutjor  component,  at  frequency  u,  /(u),  is 
produced  by  contribudons  from  all  pairs  of  points  with  a  separadun  Xu,  one  point  in 
each  apenure,  sec  Pig.6(a).  If  the  major  component  is  averaged  over  many  frames  then 
the  result  for  frequencies  greater  then  tends  to  aero.  This  is  due  to  the  phase 
difference,  modulo  2n  between  the  two  apertures  being  uniformly  distributed 
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between  ±n  with  uro  mean,  tee  Pig.7(a).  Thii  retulta  in  the  Fourier  component  per¬ 
forming  a  random  walk  in  the  complex  plane  and  averaging  to  aero. 


o'‘r|/(u)j»  0(u)  + 


</(u)>  ■  0  ;  u  >  rjk 


(46) 

(47) 


where  6(u)  la  the  Fourier  phaae  at  ftequency  u,  and  <  >  indicatea  an  entemble  average 
over  many  fnmea. 

Coniider  now  the  autocorrelation  technique  which  in  Fourier  apace  coireipondi  to 
eatimating  the  power  apeetrum.  The  nu\jor  Fourier  component  of  the  fhnge  pattern  la 
averaged  at  a  product  with  lu  complex  coqjttgate  and  ao  the  atmotpherlc  phue  contri¬ 
bution  ia  eliminated  and  the  averaged  tlgnal  ia  non-zeio,  tee  Fig.6(b). 

arg|/(u)  /(-u)j  -  0(u)  +  -  C>a  +  ^)(-u)  -  +  d>e  ■  0  (48) 

and 


<  /(u)  /(-u)  >  ■  <  l/(u)P  >  *»  0.  (49) 

Unfortunately  phaae  Information  it  not  preaerved.  Tliia  it  the  timpleat  form  of  the 
'phaae  cloture'  [11]  idea  uted  in  radio  attronomy. 

The  KT  technique  ia  a  email  modification  of  the  above.  The  m^jor  Fourier  com¬ 
ponent  of  the  fringe  pattern  it  sveraged  with  a  component  at  a  frequency  dlaplaced  by 
a  email  vector  Au.  Provided  the  vector  ditplacement  Au  doet  not  force  the  vector 
difference  -u-Au  to  be  outtide  the  tpadal  frequency  bandwidth  of  the  fringe  pattern 
then  Fourier  phaae  difference  information  it  preaerved  in  the  averaged  tlgnal,  tee 
Fig.7(b).  Again  the  atmotpherlc  phaae  effectively  forma  a  doted  loop,  tee  Fig. 6(c). 

flrgj/(u)  /(-u  -  Au)|  -  6(u)  +  -  O3  -(•  4>(-u-Au)  -  (50) 


■  6(u)  -  6{u+Au) 


and 


<  /(u)  /(-u-Au)  >  »•  0  ;  Au  <  r/k. 
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SuppoM  thil  tyitem  of  two  aporturei  it  extended  to  three  and  the  Au  value  ii 
nude  greater  than  r^/X  .  tee  Plg.6(d). 

Now 


•  * 

arg  I(u)  /(-u  -  Au)  ■  ■ 

k  < 


0(u)  +  “  Ofl  -  ^(-u-Au)  - 


•  0(u)  -  ^(u+Au)  -  +  Oc  (52) 

<  /(u)  /(hi-Au)  >  -  0  :  Au  <  rA  (33) 

It  can  be  leen  that  the  atmoipheiic  phaie  contribution  li  not  eloied.  The  reiulting 
phase  beeomei  unifomily  distributed,  see  Fig.7(c),  and  the  signal  performs  a  random 
walk  and  averagei  to  aero.  Thus  the  KT  technique  la  limited  to  frequency  differences 
Au  <  Kg/k.  If  however  the  bispectrum  average  is  performed  then  the  phase  it  agsin 
closed  and  Fourier  phue  difference  Information  is  preserved,  sea  Fig.7(d)  and  Fig.7(e). 

arg|/(u)  /(Au)  /(-u-Au)  |  ■  0(u)  +  ®x  ~  +  ♦(^u)  (34) 


+  -  Off  +  0(-u>Au)  -  ®x  + 


■  0(u)  -  ^(u-fAu)  *  0(Au) 

Thus  the  bitpeetrum,  unlike  the  KT.  can  ooutin  phase  difference  information  for  phase 
differences,  Au,  greater  than  r^/k.  It  it  apparent  that  the  phue  cloture  requirement  can 
be  restated  in  terms  of  the  spatial  frequency  vectors  which  define  a  particulur  technique 
in  Fourier  space:  the  vector  sum  of  the  spatial  frequency  vectors  mutt  be  lest  than  r^/k 
for  the  ensemble  average  of  the  signal  to  be  non-zero,  the  only  exception  being  if  fre¬ 
quency  vecton  can  be  grouped  together  such  that  all  the  groups  independently  have  a 
sum  that  it  lets  than  r^/X. 

Using  this  simplified  model  it  can  be  teen  how  the  techniques  are  related  through 
the  phase  cloture  concept  and  how  they  allow  preservation  of  Fourier  phase  informa¬ 
tion.  However,  the  use  of  a  telescope  pupil  consisting  of  a  few  tub-apertures  only 
allows  a  limited  number  of  spatial  frequencies  to  be  studied.  This  may  be  desirable  in 
certain  circumstances  but  in  general  the  fuU  telescope  aperture  it  used  to  make  use  of 
aU  the  available  light  flux  and  to  allow  many  spatial  frequencies  to  be  studied. 
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SuppoM  a  lariei  of  ihort  axpoiure  imaiei  are  recorded.  The  pth  image  may  be 
tepreiented,  using  the  incoherent  imaging  equation  in  Fourier  apace,  u 

/,(u)  ■  0(u)  Sp[\i)  (55j 

where  Sp(u)  ia  the  combined  tranifer  function  of  the  teleicope-atmoiphere  imaging 
ayitem.  The  ensemble  averaged  transfer  fbnctloni  of  the  KT  and  TC  techniquei  are 
defined  by 

7<*^(Ui,du)  •  <  5^(Ui)  5^(Ui+Au)  >  (36) 

and 

7<«^(Ui,Uj,)  -  <  5,(u,)  5;(u,‘Hi,)  S/uj)  >  (57) 

respectively. 

The  effect  of  using  the  full  aperture  on  these  transfer  functions  is,  firstly,  to  intro¬ 
duce  information  on  all  spatial  frequencies  out  to  the  telescope  diffraction  limit  and, 
secondly,  to  lower  the  transfer  flinction  for  those  frequencies  which  were  present  with 
the  three  sub-aperture  pupil  function.  The  reduction  of  the  average  transfer  functions  is 
due  to  contributions  from  combinations  of  sub-apertures,  present  in  the  full  pupil,  for 
which  phase  closure  Ic  not  realised.  For  example,  signal  contributions  from  the  two 
frequency  vectors,  u  and  uMu  in  Fig, 6(c),  but  with  both  of  the  vectors  lying  between 
different  pairs  of  sub-apertures  thus  preventing  phase  closure. 

Although  both  transfer  functions  are  four  dimensional  a  reuonable  indication  of 
their  behaviour  may  be  obtained  by  studying  two  dimensional  sub-planes  through  them 
coirespondlng  to  uiBUg^.  with  AuaAUy  for  the  KT  and  for  the  bispectrum. 

Fig.8  shows  plots  of  these  planes,  obtained  using  the  Monte  Carlo  computer  simula¬ 
tion,  allowing  a  comparison  of  their  buic  properties.  Both  possets  a  high  central 
region  for  frequencies  Uk,  Au,  and  Uj,  being  lest  than  the  'teeing  cutoff  of  approxi¬ 
mately  rg/\,  Additionally  both  planet  contain  a  high  region  around  the  axis  U|^. 
This  region  allows  information  about  frequency  differences  of  up  to  approximately  r^/X 
to  be  recovered  with  high  signal-to-noite-ratio.  It  should  be  noted  that  this  region  it 
wider  in  the  KT  cate.  This  it  due  to  the  third  term  present  in  the  blapecirum  triple  pro¬ 
duct  which  reduces  the  bitpeetrum  modulus,  falling  off  in  a  similar  manner  to  the 
short-exposure  transfer  fimctlon.  Crott-tectio.it  through  these  planet,  see  Flg.9,  show 
this  fact.  However,  the  benefit  of  this  additional  term  can  also  be  teen.  At  the  fre¬ 
quency  difference  vectors,  Au,  and  increase,  ihe  KT  signal,  in  contrast  to  the 
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biipeccrum.  reduces  In  value  and  becomes  more  randcm  due  to  the  lack  of  exact  phue 
closure.  The  remaining  high  regions  present  in  the  bispectrum  are  redundant  due  to  the 
twelve  fold  symmetry  of  the  function. 

Of  pardeulu  interest,  is  the  number  of  different  frequency  difference  vectors,  du, 
and  ujji,  high  signal  to  noise  information  can  be  acquired.  In  general  the  greater  the 
number  the  better  the  resulting  reconstructed  object  intensity  dlstribution[14],[23]. 
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Signal  to  nolle  ratio  (SNK) 


Again  the  problem  arliei  of  detcribing  the  ttatlitlct  of  complex  lignali.  The  uie 
of  the  phue  eiror  meaiuie,  6g  (Eq.(30)),  potei  many  problems  in  an  analytic  study 
due  to  the  complexity  of  the  signal  variances  required.  Consequently,  the  SNR  measure 
5NR„  (Bq.(28))  is  used,  as  in  section  4,  to  compare  analytically  the  KT  and  TC  tech¬ 
niques.  The  phase  error  measure,  Bg,  hu  hovhwver  been  studied  compuutionally,  using 
the  Monte  Carlo  simulation  previously  described  and  the  following  results  comparing 
Bg  and  the  SNK^  when  imaging  through  atmospheric  turbulence  were  obtained. 

Figure  (10)  shows  plots  of  the  ratio  of  the  phase  error  Bg  to  l/(V£SAfRJ  for  the 
sub  planes  of  the  KT  signal  and  the  biipectrum  defined  by  U|aui„  Au-Au,. 
These  plots  are  the  result  of  compuudonaily  simulating  the  effect  of  imaging  the 
'asteroid*  shown  in  Pig.3(e)  with  a  2  metre  telescope  and  0.7  urc  second  'seeing'.  Five 
thousand  independent  realisations  of  the  atmosphere  were  used  to  obtain  the  ensemble 
averages.  These  plou  should  be  directly  compared  with  those  in  Flg.S.  Although  no 
photon  noise  is  present  the  atmospheric  turbulence  induced  signal  fluctuations  result  in 
the  same  overall  effect;  the  variances  of  the  real  and  imaginary  part  of  the  signal  tend 
to  be  independent  and  equal,  with  zero  covariance,  away  from  the  axes,  whereas  near 
the  axes  the  covariance  tends  to  be  larger  with  a  small  imaginary  variance  and  a  large 
real  variance.  Again  the  SNR^  appears  to  give  a  reiionable  indication  of  the  phase 
error,  however  the  fluctuation  of  the  ratio  about  unity  near  the  axes  is  still  present,  a 
consequence  of  the  non-zero  covariance. 


The  signal  to  noise  ritloi  of  the  KT  ana  TC  in  Fourier  space  arr  defined  oy 

where 


and 


-  <  1  V**l)  P  >  -  1  <  Vui)  f;(u,+Au)  >  p. 


^(TQ 


(59) 


where 


-  •>7- 


<fhv)  ■  <  I  Vmi)  ''JCUl+Ui)  /pCUj)  I*  >  - 1  <  /,(U|)  /p(U,-i-U2)  IpiUi)  >  p, 

and  M  it  the  number  of  independent  i.nages. 

In  order  to  compare  the  SNR't  the  average*  muit  be  evaluated  for  a  range  of  fre< 
quuncy  vector*.  Based  on  the  modej  of  /loddier.  [16^  the  ^oKowlng  results  have  been 
obtained.  The  mean  ntunber  of  speckles  per  fmme,  is  assumed  to  be  much  greater 
than  unity  and  Kf  is  2.3  0fr^)\  D  bring  the  telescopu  diameter. 

The  approximadons  reqi^ired  to  etdmate  the  effect  of  atmospheric  turbulence  on 
ensemble  averages  are  repetitive,  and  tedious.  Consequently,  a  typical  quindtyt  the 
ensemble  average  mean  square  of  the  bispectrum,  is  treated  in  detail  as  a  worked 
example,  ue  in  Appendix  (3).  Ihe  effect  on  imaging  a  point  source  ac  high  light  level 
is  considared  first  The  Fourier  spectrum  of  the  pth  image,  /p(U])>  therefore  reduces  to 
that  of  the  point  spread  funcdon  of  the  telescope-atmospheru  combination,  5p(\i'h 
For  U),  Mj,  Au  >  rjkt 

<  5,(11,)  s;(u,+Uj)  5,(Uj)  >  *  7<»)(u„U2)  (60) 

<  |5,(Ui)  5^(0, +U2)  5,(U2)P  > 

e  -i-  T<*>(Ui)  7<*>(U2)  7<*^(Uj-t-U2)  (61) 

<  J,(u,)  S;(Ui+Au)  >bO  (62) 


and  for  U2,  Au  <  r/k, 


<  S,(ui)  5;(Ui+U2)  5,(U2)  >  -  7^’(U,)  <  SpiUi)  >  (63) 

<  5,(u,)  5;fu,+Au)  >  ■=  -^  T<^>(u,)  (63) 

<  |5,(u,)  5;(u,+U2)  5,(U2)P  >  “  ^  |7^^(u,)P  I  <  5,(U2)  >  P  (64) 

<  |5,(u,)  5;(u,+Au)P  >  =  -^  ir«^(u,)P  (63) 
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when  7^)(u)  is  the  normaliietl  overlap  intetral  of  two  pupils  separated  by  u.  and 
7^’^(U|,U2)  is  die  normalised  overlap  integral  of  three  pupils  separated  by  ui  and  U2. 
Consequently,  the  SNR's  an 


SNR<^ts 


2  7<^>(ui.U2) _ 

[7<J)(u,)  Ti*)(ui+U2) 


“ii  «2  >  rA 


a  ‘'OT  M2  <  r^/X 


(66) 


and 


B  0  Ml,  dM  >  rjx 

B  Va7  Am  <;  r^X.  (67) 

For  U],  Au  <  fg/},  at  high  light  level  both  techniques  have  sudstics  which,  within  the 
approximadons  made  here,  nflect  the  negadve  exponendal  behaviour  of  the  SKR  of 
the  power  spectrum[18].  For  frequencies  outside  this  range,  for  which  no  power  spec¬ 
trum  equivalent  exists,  the  KT  has  approximately  zero  SNR  whenas  the  blspectrum 
SNR  has  a  dependency.  This  nsult  disagrees  with  that  obtained  by  Wimitzer 
[23]  which  shows  no  such  dependency.  The  SNR  of  the  informadon  for  U2  >  is 
thenfore  typically  of  the  order  of  lower  than  that  for  M2  <  VX.  If  individual 
images  am  consldend  to  be  independent  then  in  older  to  obtain  equivalent  SNR  in 
both  nglons  a  muldpllcative  factor  of  order  W,  mon  frames  an  nquind  when 
U2  >  VX, 

Figure  11  shows  plots  of  a  single  plane  of  the  KT  and  bispectrum 
(U|*Ui;p  U2BU2J,  AubAu^  )  SNR’s  at  high  light  level  obtained  using  the  computer 
simuladon.  Again  the  high  ngion  around  should  be  compared  between  the  two 
plots.  The  width  of  this  ngion  gives  an  indicadon  of  the  number  of  frequency 
dlffennce  vecion  high  SNR  infoimadon  is  available  for  object  reconsirucdon.  On 
comparing  cross-secdons  for  consunt  AUj,  and  U2x  values  it  it  found  that  both  tech¬ 
niques  have  very  similar  SNR's  if  diffennce  vectors  an  lets  than  approximately  rJK, 
tee  Fig.  12.  Beyond  this  value  the  closun  phase  property  of  the  blspectrum  maintains 
its  SNR  above  that  of  the  KT  which  tends  to  fluctuate  wildly.  The  high  SNR  ngions 
around  the  axes  an  constricted  near  the  origin  for  both  techniques.  This  is  particularly 
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nodCBble  in  ihe  biipectrum  where  distinct  regions  of  low  SNR  exist,  corresponding  to 
positions  where  two  of  the  terms  in  the  bispectrum  tiiple  product  are  identicel,  hence 
increasing  the  variance  of  the  bispectrum  phase.  These  regions  are  clearly  visible  in 
Fig.  11,  for  example  around  the  line  defined  by  Uu  "Uir  In  contrast  to  the  transfer 
functions,  the  SNR's  for  the  cross*sections  shown,  with  <  r^,  maintain  a 

roughly  constant  value  for  values  of  Uij  up  to  approximately  ryX  of  the  telescope 
imposed  cutoff. 

Unfortunately,  the  need  to  observe  astronomical  objects  with  low  luminosity  and 
the  short  exposure  requirement  of  'freezing*  the  atmosphere  result  in  low  numbers  of 
photon  events  being  recorded  in  individual  images.  For  such  low  light  level  images  the 
effect  of  photon  noise  becomes  important.  The  photon  noise  dependent  variances,  see 
Appendix  (1),  therefore  need  to  be  averaged  over  the  atmospheric  statistics  in  various 
regimes.  These  regimes  being  defined  in  general  by  the  ratio,  K,  of  the  average 
number  of  photons,  fJ,  to  the  average  number  of  speckles  per  frame,  and  by  the 
value  of  the  fhequeney  vectors  U|,  U],  Au.  Again  using  the  model  of  Roddier[16], 
averaging  over  the  atmospheric  statistics  yields  the  following  approximations  for  the 
required  statistical  moments. 

For  M„  uj,  Alt  >  r^, 


fr, 


+  7<«(U,+U2)  ^  ^ 

fTj  rr, 

7^2)  (ii.) 

+  ^  7<«(Uj)  ^  f7>  1, 


(63) 


0^  =  ^  M>l,»-«l 


(69) 
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*  0- 


(70) 


For  Mj,  <  rJK, 


a  /7»  +  |<  S.(iij)  >|»  +  7<«(ui)  +  ■—  7<*)(u,+u,) 

+  7<>>(u,)  |<  S/Uj)  >1*  +  ^  7<*)(u,+U2)  |<  5,(U2)  >1* 

«,  "i 


+  ~  7<*)(ui)  7<2>(ui+uj)  l<  5p(uj)  >P  /7>  1 


(71) 


^  l<  5p(U2)  >P  /ir>  1.  H-c  1 


(72) 


^KT)  a  aJ®  +  ~  7<*>(Ui)  +  ^  T<^)(U,+AU) 


+  7<*)(ui)  7<*)(Ui+Au)  ff>  1 


(73) 


(74) 


The  mean  signali  are 

7^^1(Ui.U2)  and  7^*^(ui,Au) 

for  the  biipectnim  and  KT  reipectlvely. 

The  SNR'i  of  the  two  techniques  may  now  be  compared  for  the  various  regimes. 

For  44].  442,  ^  ^ 

SNR(^0  s  Vw 

=  ^  Kc  1  (75) 
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The  result  is  more  complicated  for  larger  n  but  the  general  form  corresponds  to  a  gra¬ 
dual  falling  off  of  the  dependence  on  .V  due  to  the  variance  becoming  more  dependent 
on  //  u  succesive  terms  in  Eq.(68}  take  precedence.  Eventually,  for  H'>  1,  the  SNR  is 
simple  dependent  on  as  the  dependent  term  of  Eq.(68)  dominates  the  variance. 
When  is  appro-ximately  equal  to  fT,  then  there  is  roughly  a  linear  dependence  of  the 
on  the  ^  dependent  term  dominates  the  variance.  The  KT  has  approxi¬ 
mately  zero  SNR  in  the  equivalent  region,  Au  >  r^/X, 


For  Mj,  du  <  and  1, 

+  Sp(U2)  >1* 


(76) 


and 


SSR^KT)  -  R  7<A71(u  du)  ^/W.  (77) 

K, 

In  particular  Eq.(76)  reduces  to  a  linear  depndence  on  IT,  equivalent  to  the  power 
f,peetrum,  for  ><2  ■  0.  u  expected.  Again  for  larger  IT  both  SNR’s  are  far  more  compU- 
cated  but,  u  before,  their  dependence  on  ^gradually  falls  off  as  tlie  terms  in  the  vari¬ 
ances  dependent  on  higher  powers  of  R  dominate.  Both  SNR's  tend  to  the  high  light 
level  value  of  unity,  at  for  the  power  spectrum,  for  U2  ■  du  aO. 

A  Monte  Carlo  simulation  hu  been  carried  out,  tee  Fig.  13,  which  reasonably 
confimu  the  ^and  dependencies  of  the  above  SNR  results. 

The  ensemble  averaged  variances  derived  can  easily  be  incorporated  into  phase 
reconstruction  algorithms  (28]  to  weight  optimally  the  information  present  in  the  aver¬ 
age  signals.  The  expression  for  the  variance  of  the  bispectrum  giA’en  in  Eq.(68)  is  par- 
dculariy  useful  at  a  weighting  function  at  it  can  be  evaluated  from  knowledge  cf  N 
and  the  ensemble  average  image  spatial  power  spectrum  at  various  spatial  ftequencies. 
This  weighting  function  will,  however,  suffer  from  the  fact  that  it  it  independent  of  the 
phase  of  the  bispectrum  at  each  point.  Two  possible  approaches  to  overcoming  this  ar: 
the  following.  Pintly,  a  larger  weighting  function  can  be  given  to  those  points  for 
which  the  bispectrum  phase  it  not  close  to  zero  on  the  basis  of  the  phase  error  tending 
to  the  equality  in  E.,  38).  Alternatively,  the  variances  and  covariance  of  the  real  and 
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inuginuy  pirti  of  the  complex  lignel,  lee  Appendix  (1),  could  be  approximated  in  a 
limilar  manner  to  the  variance  of  the  modulus  and  the  full  expression  for  the  phase 
error  defined  in  Eq.(30)  evaluated. 
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CONCLUSION 


It  it  evident  that  the  KT  end  TC  methodi  of  Image  proceiiing  are  very  closely 
related  from  both  the  mathematical  and  implementation  viewpoints.  The  extension  of 
the  KT  method  to  a  mote  general  form,  whereby  the  ^quency  difference  vector,  Au, 
can  take  on  a  wide  range  of  values  a.novea  the  TC  advantage  of  a  having  a  redun¬ 
dancy  of  phue  information  that  can  increue  the  SNR.  It  is  necessary  to  centroid 
images  before  processing  using  the  KT  method  and  this  removes  the  signal-to-noise- 
ratio,  SNR,  gain  over  the  TC  which  should  result  as  a  consequence  of  being  a  lower 
order  correlation.  Additionally,  the  need  to  centroid  and  the  development  of  an 
efficient  photon  differencing  TC  algorithm  result  in  Implementation  times  being 
equivalent,  despite  the  KT  method  being  a  Iowa-  order  correlation. 

A  more  appropriate  measure  of  the  phue  error,  0£,  of  a  complex  signal  hu  been 
suggested.  This  has  been  compared  with  the  traditional  SNR  bued  measure  whereby 
the  phue  error  is  approximated  by  l/('^SiV/f).  It  is  found  that  this  gives  a  reasonable 
indication  of  the  phase  error;  however,  in  certain  circumstances,  for  smell  frequency 
difference  vectors,  the  value  predicted  by  this  means  can  be  considerably  in  error.  It  is 
seen  that  the  SNR  can  be  used  to  provide  an  upper  bound  on  the  phase  error. 


Speckle  Interferometry,  the  KT  correlation  and  the  TC  have  been  seen  to  be  con. 
nected  through  phase  closure  relationships.  By  considering  phase  closure  in  terms  of 
spatial  frequency  vectors  forming  closed  loops  in  the  telescope  pupil,  it  hu  been  possi. 
ble  to  evaluate  many  ensemble  averages  of  image  correlations  over  the  statistics  of 
atmospheric  turbulence.  A  computer  implemenution  of  a  generalisation  of  the  tech* 
nique  of  Goodman  and  Belsher[18]  has  allowed  photon  noise  dependent  SNR's  to  be 
investigated.  Using  these  developments,  together  with  a  computer  simulation  of  the 
effecu  of  imaging  through  turbulence,  the  SNR's  of  the  TC  and  KT  methods  have 
been  compared  with,  in  particular,  their  dependency  on  and  57! 

It  is  found  that  the  the  highest  SNR  regions  of  both  the  TC  and  KT  average  sig. 
nals  are  of  equivalent  value  and  extent.  They  correspond  to  using  frequency  difference 
vectors,  au,  U2,  that  have  magnitudes  that  are  lest  than,  or  of  the  order  of,  rj\.  In  this 


81 


-  34- 


region  both  SNR't  are  linearly  dependent  on  the  number  of  photoni  per  ipeckle  at 
very  low  light  levelt  and  approximately  unity  ai  high  light  level.  They  are  reflecting 
their  cloae  uioiciation  with  the  negative  exponential  itatiatlci  of  the  power  spectrum. 
Outalde  these  high  SNR  regions  the  KT  signal  it  iniigniflcant  and  the  bitpectrum  SNR 
ihowt  a  dependency  at  veiy  low  light  levelt  and  a  1/V^  at  high  level.  Again 

thit  it  tuggetting  the  effective  equtvalance  of  the  two  techniquet  at  for  typical  values 
of  W,,  of  the  order  of  1000.  the  additional  information  present  in  the  bitpectrum  would 
require  many  more  flramet  to  be  procetted  to  bt  of  any  use.  Retultt  from  a  Monte- 
Carlo  type  computer  limulation,  which  enabled  the  ensemble  average  over  the  atmos¬ 
pheric  ttatitiict  to  be  evaluated  within  the  conitraintt  of  the  model  used,  agree  well 
with  the  resulu  obtained  from  the  analytic  study  bated  on  phue  cloture  reladont. 

This  work  wu  supported  by  the  UK  Science  and  Engineering  Research  Council 
(OR/D  92332)  and  the  US  Army  (DAIA  45-8S-C-0028). 
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Appendices 


Appendix  (1) 


A  computer  implemenudon  of  t  genenliittion  of  Ooodtntn  end  Belsher'i 
method  [18]  for  tvereging  ligntl*  in  Fourier  ipece  in  the  preienee  of  photon  noite  hu 
been  uied  to  yield  the  following  reiulu.  Coniider  inuging,  with  t  meen  number  of 
photon  events  per  frame  of  an  intensity  distribution  which  has  a  high  light  level 
image  with  the  Fourier  transform  /(u). 

The  KT  signal  of  the  pth  frame  it,  in  Fourier  space, 

Z)<*”(ui,Au)  -  D;,(u,)  D;(Ui+Au)  -  I>p(Au),  (Al.l) 

resulting  in  the  following  expressions  for  the  second  order  ttatittict.  The  second 
moment  of  the  teal  part  of  the  signal  it  given  by 

1/4  <  [  Z)j*^(Ui,Au)  +  ]*  >  (A1.2) 

and  averaging  over  the  photon  noite  sutittict  gives 

IM  [  /(-U|)  /(-U|)  /(ui+Au)  /(ui+Au)  *  ee 
*  2  /(-U|-Au)  /(-U|)  ;(U|)  /(U|+Au)  ,9" 

+  /(-V|)  /(-Ml)  /(2ui+2Au)  ^  +  ec 
+  2/(-u,)'(u,)^ 

+  2  /(-Ui)  /(Au)  /(Ui't’Au)  R*  *  cc 
+  2  /(-Uj-Au)  /(-Ui)  /(2ui+Au)  fP  +  cc 
+  /(-Ui-Au)  /(-U|-Au)  /C2ui)  R*  *  cc 
+  2  /(-Ui-Au)  /(U|+Au)  R* 

+  /(-Au)  /(-Au)  R**cc 
+  /(-2ui)  /(2ui+2Au)  P  +  rc 
+  2  /{-2a.-Au)  /(2ui+Au)  R* 

*2R‘]. 

where  cc  indicates  the  complex  conjugate.  It  should  be  noted  that  each  cc  is  the  com* 
plex  conjugate  of  the  whole  of  the  term  to  its  left,  including  the  sign. 


83 


-36- 


The  lecond  moment  of  the  imigintty  part  of  the  ilgnel  it  given  by 

-1/4  <  [  )*  >.  (A1.3) 

and  averaging  over  the  photon  noiie  itatittict  gives 

-1/4  [/(-U|)  /(-Ui)  /(UffAu)  /(Ui+Au)  +  ec 

-  2  /(-Ui-Au)  /(-Uil  /(Ui)  i(tti+Au)  ff* 

+  I{-U\)  /(-Ui)  /(2M1+2AU)  FP  +  ce 

-2/(-tti)/(u,)^ 

2  /(-U|)  /(Au)  /(Ui-t'Au)  ^  -f  cc 

-  2  7(-U|-Au)  /(-Ui)  /(lui-i-Au)  P-*-  ee 

-  /(-Ui-Au)  /(-Ui-Au)  /(2ui)  P  *  ct 

-  2  /(-U|-AU)  /(u,+Au)  P 
>  /(-Au)  /(-Au)  P-¥cc 

+  /(-2ui)  /(2ui+2Au)  P*cc 
■  2  /(~2U|-Au)  /(2u|-i’Au)  P 

-2/7’]. 

The  Aret  ,'oint  moment  of  tlie  real  and  imaginary  paru  of  the  lignal  it  given  by 

-y/4  <  [  -  0),”^*(Ui,Au)  ] 

.[  Dj*^(Uj,Au)  +  DJP^’(Ui,Au)  ]>.  (A1.4) 

and  avenging  over  the  photon  nolie  itatiitici  gives 

•JiA  [  /(-U|-Au;  /(-U|-Att)  /(Ui)  Hvt)P  -  cc 

-  /(-Uj-Au)  /(-Ui-Ah)  /(2ui)  ^  -  cc 
4  2  /(-Ui-Au)  /(-Au)  /(U|)  P  -  cc 

+  /(-2ui-2Au)  /(Ui)  /(Ui)  P  -  cc 
f  /(-Au)  /(-Au)  P-  cc 
+  /(-2ui-2Au)  /(2ui)  ^  -  cc  ]. 
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The  TC  ligntl  of  the  pth  fnme  is,  in  Fourier  space, 

f?^(Uj,Uj)  ■  Z?^(Ui+Ui)  £),(U2) 

-  |D/u,)l^  -  |Dp(u,+uj)|*  -  !Z),(U2)I*  +  2Np. 

resulting  in  the  following  expressions  for  the  second  order  sudslics. 
The  second  moment  of  the  real  part  of  the  signal  is  given  by 

1/4  <  [  iDj^(Ui,Uj)  +  :)^*(Ui.U2)  ]*  > 
and  averaging  over  the  photon  noise  sutlitlcs  gives 

1/4  [  /<-Ui)  /(-Ui)  /(-Ui)  A-ua)  /(2ui+2ua)  fJ*  *ec 
+  2  A-ui)  A-Ui)  A-ua)  /(Ui)  Aui+uj)  f^*ec 
+  A“«i)  A-U|)  /(-2ua)  /(Ui+ua)  Aui+ua)  f?  *  ec 
+  2  A-Ui)  A-ua)  /(“vO  /(«»)  /(“i+Ua)  l^*e£ 

+  2  A-Ui)  A-U;)  /(Ma)  AUi)  ^ 

+  2  A-Mi-Ui)  /(“U|)  /(-uj)  Aui)  Aui+2ua)  ^  +  ec 
+  2  /(-M|-Ua)  A-Ui)  A-Ma)  AU|+Ua)  Aui+«a)  ^  +  ee 
+  2  A-UpUj)  A-Ui)  A-Ui)  Aua)  /(2U)+ua)  P*ce 

*  2  A-Uru»)  A-ui)  Aui)  Aui+ua)  P 

+  2  A-ui-ua)  /(“Ui)  /(ua)  Aoi-ua)  /(ui*ua)  f^*ce 
+  /(-ui-ua)  A-ui-ua)  Aua)  Aua)  A2ui)  *  ec 

*  2  A-Ui-Ua)  /(-Ua)  /(ua)  Aui+ua)  P 

*  /(-Ui)  A-Ui)  A-2ua)  A2ui+2ua)  *  cc 
+  2/(-ui)A-ui)/{u,)Aui)/7' 

+  4  A-ui)  A-Ua)  Aua)  Aui)  ^ 

+  2  A-U|)  A-2ua)  /(ua)  /<Ui+Ua)  +  cc 
+  2A-0i)AU|)^ 

+  2  A-U|)  Aua)  Aui-Ua)  f^*cc 
+  2  A--Mi-ui)  A-Ui)  /(-Ua)  A2ui+2ua)  *  cc 
+  4  A-Ui-Ua)  /(-ui)  /(ui)  /(Ui+Ua) 

+  2  /(-U|-Ua)  /(-U|)  /(ui-ua)  /(ui+2ua)  *cc 
+  2  /(-U|-ua)  A-Ui)  /(Ua)  /(2ui)  fP*ce 
+  2  /(-U|-Ua)  A-Ui)  /(2ui+ua)  *  cc 


(A1.5) 


(A1.6) 
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■f  2  /(-ui-uj)  /(-Uj+ttj)  /(-Uj)  /(2U|+Uj)  /J*  +  ee 
+  2  A-urUj)  /(“U|-*'Ui)  /(U|>U})  /(uj+uj)  ^ 

+  2  /(-urOiJ  /(-urui)  /(U)-Hii)  /(ui+uj) 

+  A-Ui-tn)  /(-UrW)  AJuj)  /(2ui)  fJ*  *ee 
+  2  /(-urUj)  A-’Jj)  /(Li+2ui)  /?*  +  «< 

+  4  /(-uj)  /(uij  Aui+Ma)  ^ 

+  2  A-ui-uj)  /(tti+ui) 

+  2  /(-ur2u|)  A-Ui)  Aui)  Aui+2mj)  ^ 

+  2  /(-Ur2u»)  A-ui)  Aui)  A3u,+U2)  ^  +  ee 
+  2/(-u,)A-u»)Amj)Auj)A^ 

+  2  A-uj)  Aiii)  ^ 

*  A-2ui)  A-uj)  A-uj)  A2ui+2ui) 

+  2  A-2urUj)  A-uj)  Auj)  A2ui+uj)  ^ 

*  A-Ui)  A-Ui)  A2ui)  fP  *ee 

f  2  A“«i)  A-ui)  Aui+Uj)  fp  +  ee 

*  2  A-ui+uj)  Aui-uj)  ^ 

+  A-Ui-Ui)  A2Mi*2ua)  fP  *ce 

+  2  A"Ur2uj)  A-Ui+Uj)  /(2ui+ua)  fP  ^cc 

*  2  A-Ui-luj)  AU|+2ui)  fP 
+  A-ui)  A-Uj)  A2ua)  fP  +  cc 

+  A-2ui)  A-2ua)  A2ui+2ua)  FP  ■*ec 
+  2  A“2uruj)  A2ui+uj)  FP 
*2FP 

*  A-tti)  A-«ii)  A-ua)  A-ua)  Aui+ua)  /(Ui+«i)  ^  +  cc 
+  2  A-Ui-Ua)  A-U|)  A-Ua)  Aual  /(Ui)  Aui*«a)  FP  ]■ 

The  lecond  moment  of  the  Imaginuy  put  of  the  lignal  It  given  by 

-1/4  <  [  Z?^^(U,.U2)  -  0^^‘(U„U2)  ]*  > 
tmd  averaging  over  the  phcton  nolle  tutistici  gives 

-1/4  ( A-ui)  A-U|)  /(-ua)  A-Ua)  A2ui+2ua)  FP  +  ce 
+  2  A-Ui)  A-U|)  A-ut)  Aui)  Aui+ua)  FP  +  ec 
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+  A-Ui)  /(-ui)  /(-2u,)  /(Ujt-uj)  /(u,+u,)  f^  +  ee 
+  2  /(-ui)  /(-Ui)  /(-Oi)  /(Uj)  7(U|-t-U|)  ^  cc 

-  2  /(-tti)  /(-uj)  /(uj)  /(ui)  fr 

-  2  /(-urM»)  /(“«,)  7(-Ui)  /(ui)  /{U|42u»)  ^  +  <c 

+  2  /(-U|-Ui)  /(-U|)  /(-Uj)  7(ui+ui)  /(U|+U})  ^  +  ft 

-  2  /(-u,-ut)  /(Mil)  A-u,)  /(u^  /(2ui+ui)  ^  +  cc 

-  2  /(-U|)  /(u,)  /(ui-t-m)  >?» 

-  2  /(-U1-U2)  /(-U|)  /(Uj)  /(Ui-Ua)  /(UftUj)  fP  +  cc 
+  /(-uru»)  A-ui-uj)  Auj)  /(uj)  /(2tii)  ^  +  cc 

-  2  /(-u,-Uj3  A-tti)  /(Ui)  /(u,t-uj)  fP 

+  /<-Mi)  /(-Ui)  /(-2ua)  /<2ui+2ui)  fP cc 
2  /(-ui)  /(-U|)  /(ui)  /(U|)  fP 

+  4/(-u,);(-ui)/(u,)/(ui)^ 

+  2  /(-u,)  /(-2ua)  /(Ui)  /(ui+g,)  fP  *  ee 

-  2  /(-U|)  /'U|)  fP 

-  2  /(-Ui)  /(Ui)  /(ui-ui)  ^  +  cf 

+  2  A-Ui-Ui)  /(-oi)  A-Ui)  /(2ui+2uj)  ^  +  cf 

*  4  A-uro>)  A-Ui)  Atti)  /(u.+uj)  FP  +  cf 

-  2  A-U)-Ui)  A-Ui)  Aui-uj)  /(ui+2iij)  FP*cc 

♦  2  A-ttrui)  A-U|)  Auj)  A2u|)  fP*ce 

-  2  A-Ur«;)  A-«i)  /{2u,+Ui)  ^  +  cf 

-  2  A-urUj)  /(-Mj+u,)  /(-uj)  A2ui+uj)  FP  cc 

-  2  /(-ui-ui)  A-Ui+u,)  Aui-uj)  Aui+Ui)  FP 
•r  2  A-U|-Ui}  /(-U|~Ui)  AU|-f|l]}  /(Uf)>U|)  fP 
+  A-Ui-Uj)  /(-U,-Ui)  A2ua)  A2ui)  FP*ec 

-  2  A-Ui-ui)  /(-Ui)  /(U|+2ui)  fP  *  cc 
+  4  A-ui-ui)  A-uj)  Au,)  /(u,*ui)  FP 

-  2  /(— U|-Ui)  Aui^U])  FP  *  ce 

-  2  /(-ui-2ui)  A-ui)  /(ui)  /(u.+2ua)  FP 

-  2  /{-ui-2u,)  A-ui)  /(u,)  /(2u,•^u,)  ^  +  ce 
••■2/(-ui)A-u,)/(Ui)/(u,)^ 

-  2  A-ui)  /(Ui)  )9* 

+  /(-2u,)  A-u,)  A-u,)  /(2u,+2ui)  FP*cc 

-  2  /(-2ui-u,)  A-Ui)  /(Ui)  A2ui+Ui) 

+  /(-Ui)A-u,)/(2u,)/^  +  fc 
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+  2  /(-ui)  /(-uj)  /(ui+ui)  fP  +  cc 

-  2  /(-uj+uj)  /(ui-uj)  ^ 

+  /(-U|-nj)  /(-Ui-Uj)  /(2ui+2uj)  fP^ec 

-  2  /(-bi-2uj)  /(-«i+uj)  /(2ui+ui)  IP  *  ce 

-  2  /(-ttfluj)  /(U)+2ui)  A/* 

+  /(Hij)  /(2uj)  ^  +  ee 

+  /(-2ui)  /(-2uj)  /(2ui+2ui)  fP  *  ec 

-2/(-2urUi)A2u,+Uj)^ 

-2^ 

+  /(-Ui)  /(-U|)  /(-Uj)  /(-Uj)  /<Ui+Ui)  /(Ui+Uj)  ^  « 

-  2  /(-Ui-Uj)  /(-Ul)  /(-ttl)  Aui)  /(U,)  /(U,-MI,)  ^  ]. 

The  flrtt  Joint  moment  of  the  reel  and  imaginary  parti  of  the  lignal  it  given 


-JIA  <  r  D^(u,.ui)  -  Z)^‘(U,.U2)  ] 

,[  1>  (A1.8) 

and  averaging  over  the  photon  noiie  iiatiitici  givei 

-JIA  [  -  ;(-Ui)  /(-Ml)  /(-uj)  /(2ut+2ua)  fP-ce 

-  2  ;(-ui)  A-ui)  /(-ui)  /(Ui)  /(tti+ui)  ^  -  ee 

-  /(-a,)  /(-ui)  /(-2uj)  /(u,+uj)  /(ui+ui)  fP-ec 

-  2  A-ui)  /(-Uj)  /(-uj)  /(Uj)  /(Ui+Uj)  -  ec 

-  2  /(-Ui-Uj)  /(-U|)  /(-Uj)  /(Ui+Uj)  /(Uj+Uj)  fP  -  CC 
+  /(-Ui-Uj)  /(-U,-Uj)  /(Uj)  /(Uj)  /CUi)  ^  -  CC 

-  /(-U,)  /(-U|)  /(-2uj)  /{2u,+2uj)  ^  -  re 

-  2  /(-ui)  /(-2uj)  /(Uj)  /(Ui+Uj)  ^  -  ee 

-  2  /(-Uj-Uj)  /{-U|)  /(-Uj)  /(2ui+2uj)  -  cc 

+  2  /(-ui-uj)  /(-U|)  /(Ui)  /(2ui)  -  cc 

+  /(-Uj-Uj)  /(-tti-iij)  /(2uj)  /(2ui)  IP  ~ce 

-  /(-2ui)  /(-uj)  /(-uj)  /(2ui+2uj)  FP  -cc 
+  /(-Ui)  /(-Ui)  /(2u,)  FP-ce 

+  2  /(-U|)  /(-Uj)  /(U|+uj}  A**  -  cc 

-  /(-Ui-Uj)  /(-Ui-Uj)  /(2U|+2uj)  ^  -  cc 
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Appendix  (2) 


Centrolding  fnunei  before  performing  the  KT  technique  reiults  in  eirora  being 
introduced  into  the  eniemble  avenge  when  photon  noiie  it  preient. 


The  pth  low  light  level  image  li  deicribed  by  Eq.(18)  and  iti  Fourier  trantform 
by  Eq,(19).  The  following  expreiiion  maybe  written  relating  the  thllt,  X^,  from  the 
image  centroid  of  the  photon  event  diitribution.  llte  image  centroid  being  the  'centre 
of  gravity'  of  the  image  photon  event  diitribution. 


(A2.1) 


where  Xg  it  any  point  in  the  frame  used  at  a  reference,  Hence, 

I  *k 

(A2.2) 

The  KT  lignal  in  Fourier  space  it 

flpCUi)  Z);(u,+Au)  -  o;(Au),  (A2.3) 

and  if  Dp(,v)  it  centroided,  using  the  Fourier  shift  theorem,  to  give, 

£)p(u)  exp(-2rt/XpU),  (A2.4) 

then  the  resulting  KT  signal  formed  is 

Z)p(Ui)  exp(-2?t/XpUi)  Dp(ui+Au)  exp(2>t/Xp(Ui+Au) 

-  f3p(Au)  exp(2tt/XpAu).  (A2.5) 

Using  the  method  of  Goodman  ind  Beliher  (18]  the  ensemble  average  over  the  photon 
event  diitribution  yields 

<  Np  (\p-i)  /(u+-^)  r(u+Au-^)  > 

/Vp  Np 
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where  /(u)  It  the  object  Fourier  tnntform,  /(u),  normalised  by  f(0)  and  <  >  indicates 
in  eniembie  average  over  the  Poisson  itatiitict  of  the  photon  events,  Obviously  the 
bias  terms  cincel  as  desired.  However,  the  resulting  signal  differs  from  the  required 
result, 

<  Aui+Au)  >.  (A2.7) 

obtained  fn>m  a  series  of  photon  limited  images  of  a  stationary  object  which  is  not 
centroided  before  performing  the  KT  average.  The  average  over  the  bias  term  alone 
yields  a  result  which  is  completeiy  general  for  any  frequency  Au.  Consequently  the 
same  result  maybe  used  to  describe  the  resulting  image  obtained  by  summing  a  series 
of  centroided  images  which  contain  photon  noise. 

<  £(0,(u)exp(-2j^/XpU)]  >  -  <  N.  (A2.8) 

/* 

where  £[]  indicates  the  expectation  value  of  the  signal  due  to  the  photon  event  distri¬ 
bution. 
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Appendix  (3) 


Roddier[16]  hai  developed  a  mooel  which  may  be  extended,  in  conjunction  with 
tome  tlmple  aitumptioni  conceminj  the  ttatistics  of  atmotphere  induced  wavefront 
perturbadont,  to  approximate  the  effect  of  atmospheric  turbulence  on  ensemble  aver- 
aset  of  intensity  distribution  correlations  in  fieneral. 

Imaging  through  atmospheric  turbulence  with  a  large  tingle  aperture  telescope  it 
modelled  by  describing  the  telescope  pupil  function,  at  a  distribution  of  'seeing 
cells'  of  diameter  rg.  The  wavefront  phase  it  assumed  constant  across  each  cell  and 
the  complex  amplitude  in  the  pupil,  ACx),  is  attumed  uncorrelated  between  cells.  The 
complex  amplitude  modulus  it  attumed  to  be  unity.  The  Fourier  transform  of  the 
intensity  distribution  of  an  image  of  a  point  source  is  tlrerefore  described  by, 

J  Afg(x)  A/J(x+?iu)  A(x)  ,4*(x+Xu)  dx 

- .  (A3.1) 

f  |W,(x)l*dx 

puplt 

This  integral  may  be  approximated,  using  the  simple  model,  by  a  sum  of  complex 
exponentials  corresponding  to  the  complex  amplitude  difference  between  pairs  of  'tee¬ 
ing  cells'. 

SpW  -  T  cxp[/(<5,-4.;)]  (A3.2) 

where  it  the  wavefront  phaie  difference  between  the  sth  pair  of  'seeing  cells' 

separated  by  ku.  There  are  approximately  A'(u)  such  pairs  in  the  pupil. 

In  order  to  study  the  signal  to  noise  ratios  of  the  KT  and  TC  techniques  a  number 
of  different  correlationt  need  to  be  averaged  over  the  atmospheric  ttatistics.  Using  the 
above  model  a  particular  correlation  it  treated  to  try  to  indicate  the  form  of  the 
approximations  required  in  general.  These  and  similar  approximations  ate  used  to  aver¬ 
age  all  the  correlationt  needed  for  the  investigation  allowing  the  form  of  the  dominant 
terms  to  be  determined. 

Consider,  as  an  example,  averaging  the  modulus  squared  of  Uie  bispectrum  of  a 
point  source,  for  frequencies  u,,  uj  >r^K  .  This  term  it  of  particular  importance  as  it 
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tllowi  the  variance  of  the  biipectrum  to  be  evaluated. 

<  |5/Ui)  S;(ui+U2)  5,(U2)P  > 

,  «(u,)  W(u,+uJ 

•  <  III  “  «*♦) + A®,  -  ^r)  +  -  <»,)] 

W(0)’  y.1  ,.1  «■! 

W(“l)  WtUi-MIJ  N{ui) 

•  2  t  t  expW®,-  -  O  •)  +  i(«>/  -  ®,)  +  X®/  -  ®/)]  >  (A?.3) 

^■l  rml  «■! 

In  general,  thii  expreuion  conreipondi  to  aununing  over  contilbulioni  from  lix 
different  pain  of  'teeing  celli’  and  for  tuch  ccmbinationi  the  entemble  average  tends 
to  zero.  However,  certain  other  combinadoni,  for  which  'seeing  cells'  occur  in  more 
than  one  pair,  result  in  'phase  closure’  and  a  non>zero  ensemble  average  contribution. 
Such  combinations  all  result  in  the  phase  closure  relationship 

+  +  (Oi-®/)  +  (®i“®,)  -  0.  (A3, 4) 

which  may  be  realised  in  a  number  of  different  ways.  To  see  how  these  arise  and  to 
obtain  their  contribudion  to  the  ensemble  average  it  it  inttrucdve  to  consider  imaging 
through  a  telescope  pupil  with  an  atmosphere  induced  complex  amplitude  dittiibudon 
which  it  delta  correlated. 

Q(x)-<A(x‘)  4V+x)>‘-5(*)  (A3.5) 

This  corresponds  to  making  N(0)  approach  infinity  so  that  the  lummadon  approxima* 
don  of  eqn.  (A3.3)  again  becomes  an  Integral  expression. 

<  l3/ui)  5;(ui+«s)  > 

4JJM  W,(Xi)  /fJCxi+Xu,)  /C(*i->.(ui+uj))  W,(x>)  Wl(xj+^Uj) 

.  H,(.x^+h^^)  //^(xj)  «,(xr-k(U|+M,))  w;(x,) 


<  A{xi)  4'(X|-t-Xui)  A(Xa)  4*(Xj-X(0;*U;))  4(xa)  4’(xj-t-Xuj)  A’(X4)  /4(X4+Xui) 

4*(Xj)  4(X5-\(Ui+uiy)  A’lXt)  4(x«+Xuj)  >  dX)  4X)  4xj  0x4  dxf  4x«, 
where  the  normalisadon  factor  Y  is  defined  by 

-[/II  W,(xi)l’  |«,(x,)|»  <fx,  rfx,  dx,  ]*  (A3.6) 
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By  ivenging  over  the  complex  amplitude  ttatiitici  the  phue  cloture  non^zero 
ensemble  average  contribution  terms  can  be  approximated.  Since  the  complex  ampli¬ 
tude  it  assumed  to  possets  unit  modulus  iu  higher  moments  cannot  uc  generated  using 
Reed's  theorem  [31]  which  it  only  applicable  if  the  complex  amplitude  it  a  gaustian 
random  variable.  Instead  the  high  order  moment  it  broken  down  into  combinations  of 
lower  order  moments  on  the  buit  of  phue  closure  contideradont.  This  break  down  of 
the  higher  order  moments  it  perhaps  more  euily  seen  if  the  phase  cloture  relationship 
it  thought  of  in  terms  of  spatial  frequency  vectors  forming  closed  loops  in  the  tele¬ 
scope  pupil,  tee  Fig.  14. 

(a)  Consider  terms  for  which 

(«>,  -  «»i)  +  («;■  -  0,0  ■  0* 

(O,  -  <)  +  (O;'  -  0,0  -  0, 

(0,-0;)  +  (0;-0y)-0, 

These  correspond  to  the  following  non-tero  contribution  to  th.e  ensemble  average  over 
the  complex  amplitude  statittict,  tee  Pig.l4(b}. 

<  4'(X4)  >  <  4(iu+Xu,)  4*(X|+Xu,)  >  <  A(xi)  4‘(xj)  > 

,  <  4(Xr-X<U|+Ui))  4'(*»“X(Ui+Ua))  >  <  4(Xa)  A’(x«)  >  <  A(x,+Xu»)  A’tXa-rXua)  > 

-  8(X4-x,)  8(xrXi)  8(xr-x,)  (A3.7) 

Resulting  in  the  following  contribution  to  the  average  modulus  squared  of  the  bitpec- 
trum  which  is  non-zero  for  all  spatial  frequencies  out  to  the  diffraction  limit  of  the 
telescope.  For  an  unapodited  aperture  this  it  simplified  to  a  product  of  terms  depen¬ 
dent  on  the  normalised  diffraction  limited  transfer  function,  7^*(u). 

y”  J  1  J  |/f,(xi-X(ui+U;))P  W,(Xj)|*  U/,(X|+Xuj)|*  dx,  dx,  dxy 

-  A  TO’iui)  TO'Cu,)  7<>J(u,+ua)  (A3.8) 

where  the  mean  number  of  speckles  per  frame,  n,  is  defined  by 

J  \HM\*  4x 

- 

J  IC„(x)p  dx 


94 


.47  - 


(b)  Coniider  tarms  for  which 

<♦,  -  ■  0, 

(♦,  -  ♦i)  +  -  0,0  +  -  •!)  +  <©y  -  ©y)  ■  0, 

These  cotreipond  to  two  terms,  see  Fig,  14(c),  that  have  identical  non*  zero  contribu¬ 
tions  to  the  ensemble  average  over  the  complex  amplitude  statistics,  provided  the  pupil 
is  symmetric. 

For  example, 

<  4(Xi)  /l‘(X4)  >  <  >t(iu+3tui)  4‘(K,+Xtt,)  >  <  A(Xi)  A*(x«)  > 

.  <  A(Kii  X*(Kr-^Ui+Ui))  >  <  A(*j-X(Ui+Uj))  A*(xj+Xui)  >  <  A(XrrXuj)  A‘(X5)  > 

-  d(x4-X|)  8(x<-xi)  8(X|-X|-Xui)  8{xr-x,-X(ui+uj))  8(xi-xi-X(ui+tt>>-Xui).  (A3.9) 

results  in  the  following  contribution  to  the  average  modulus  squared  of  the  bispectrum 
which  is  non-zero  for  all  spatial  frequencies  out  to  the  diffhiedon  limit  of  the  tele¬ 
scope.  This  simplifies  for  the  case  of  an  unapodised  aperture  to  a  product  of  terms 
dependent  on  the  normalised  diffraction  limited  transfer  function,  7^>(u),  and  the  nor¬ 
malised  overlap  integral,  T^)(U],u,tui),  of  four  pupils  with  separations  fTom  an  arbitrary 
origin  of  zero,  Xuj,  X(U|^])  and  X(u|4-Ui>fXui. 

1  f  W,(x,)|«  W,(x,+Xu,)l»  As, 

Ani 

^  IW,(xj-Xu,)|’  W,(xrX(ui+uj))|*  |ff,(Xj+Xttj)|’  d*i 

■  7^’(ui)  T*’(Ui,ui+Uj)  (A3. 10) 

A) 

The  symmetry  of  the  bispectrum  squared  modulus  results  in  further  terms  for  the 
same  phase  closure  realisation,  namely 

^  7®(uj)  r<'>(ui,U|+uj)  and  -2- T®>(ui+aj)  T^'^tui.uj).  (A3.ll) 

"1 

(c)  Consider  terms  for  which 

(©,  -  ©i)  +  {©,  -  ©;>  +  (©,  -  ©;> .  0. 

(©^.  -  ©^.)  +  (©;.  -  ©^,)  +  (©;.  -  ©^.) .  0. 
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Thete  correspond  to  four  terms  that  have  identical  non>zero  contributions  to  the  en¬ 
semble  average  over  the  complex  amplitude  statistics.  For  example, 

<  A(xj)  A*(xi+Xui)  >  <  ,S(Xi)  M*(Xj-^U|+Ma))  >  <  /i(Xi)  > 

<  4(X4+XUi)  A‘(X^  >  <  4(*5-X(U|+Uj))  ,S‘(X4)  >  <  4(X«+)4Ij)  A*(Xj)  > 

-  e(xi-Xj+Jui|)  8(Xi-xi-Vui+ua))  8(xt-xa+Xua)  8(Xf-X4-Xui) 

.&(X4-XfrX(U|+Ua))  8(xr-X4-Xua)  (A3. 12) 

Another  example  is  shown  dlagrammatically  in  Pig.  14(d)  Again  a  diffraction  limited 
contribution  results  which,  after  combining  all  terms,  is  slmplifled  by  considering  an 
unapodlsed  aperture,  to  yield 

[  Y  l«,(X|+X(ui+ua))l‘  4x, 

■ -4  I<’)(tt„Ui)  7<»)(tt,.u,)  (A3. 13) 

The  normalised  diffraction  limited  bispectrai  transfer  function  7^’(U|,U})  it  the  overlap 
integral  of  three  pupils  with  separations  from  an  arbitrary  origin  of  zero,  Xui  and 

XCui+ua). 

The  cloture  phase  realisation  being  considered  in  tl-it  section  hat  tome  redun¬ 
dancy  with  the  previous  realisations,  (a)  and  (b).  For  example,  in  eqn.  (A3. 10)  if 
xit'X.U)  >  X]  then  the  phase  closure  conditions  of  (c)  are  realised.  The  redundancy  mutt 
therefore  be  eliminated  by  subtracting  the  relevent  duplicated  terms.  In  this  example 
this  corresponds  to  tubtraedng 

7<*l(u,,ui+u,). 

Compensating  for  all  the  redundancy  requires  subtracting  the  following. 

^  [  f'*’(U}.Ui+u,)  +  r^’to^ui-ruj)  +  7<*>{u,.ttj)  ]  * -^  T^^Ui.uj)  (A3.14) 

Figures  14(e)  and  14(f)  shown  examples  of  phase  closure  relations  which  shown 
redundancy. 

(d)  Finally,  consider  realising  the  cloture  phase  of  eqn.  (A3.4)  without  any  of  the 
equality  conditions  of  (a)  -  (c)  being  true.  Averaging  over  the  complex  amplitudes 
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yield!  40  termi  which  have  the  following,  identical,  non-aero  contribution!  to  the 
eniemble  average  provided  the  pupil  ia  lymmetric. 


I«,(*i)l'  l«.(*i+Xu,)r  l«.<*,+X(u,+u,))l*  4*1 


r 


(A3. 15) 


The  eniemble  averaged  niodului  iquared  of  the  biipecirum  therefore  coniiitt  of  a 
lum  of  many  termi  with,  in  particular  differing  dependenciei  on  K,,  Howevv,  it  hai 
been  aaiumed  that  K,  li  very  much  greater  then  unity  and  to  In  this  limit  h  rmi  which 
are  inveriely  proportional  to  the  loweit  power  of  n,  will  euily  dominate.  In  thii  exam¬ 
ple  the  contribution  of  the  term  defined  by  Eq.(A3.g)  will  dominate  all  othen. 

Thii  approach  hai  been  applied  to  all  the  eniemble  averagei  over  the  atmoipheric 
itaditini  required  for  the  SNR  study. 
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Figure  captions 


Figure  (1) ;  Correipondlng  lub-planes  of  (t)  the  blspectrum  (  UibUu,  U2"U2,  )  and  (b) 
the  KT  (  u,»ui„  AubAUx  )  four  dimeniional  transfer  functions  of  a  circular  aperture 
lens. 

Figure  (2)  :  Pictorial  representation  of  correlation  techniques  in  the  image  plane,  (a) 
The  autocorrelation,  (b)  the  Knox-Thompson  correlation  for  a  single  plane,  AuaAU;c 
(c)  a  tingle  blspectrum  plane,  U}b2Aub2AUj,  represented  at  a  double  correlation  in  the 
image  plane,  and  (d)  coTrelttlons  required  for  phase-gradient  technique,  a«a,.  The 
asterisk  (*)  Indicates  the  correlation  of  the  two  dimensional  images  and  multiplication 
by  a  function  it  indicated  by  the  overlayed  dashed  axes  and  1-D  function,  The  func¬ 
tions  are  constant  in  the  orthogonal  direction.  The  real  and  imaginary  pans  of  the 
correlations  are  shown  separately. 

Figure  (3)  t  Plots  (a)  shows  the  theoretical  urror,  ®|.  for  a  tingle  image,  on  the  phase 

difference  estimates  retuldng  from  using  thu  bispectrum  plane  (U2(^U;t) 

and  the  KT  plane  (Uiy«0,  AUj,»0).  Thu  signal  to  noise  ratio,  5NR„,  for  the  same  two 

planes,  are  shown  in  (b).  Plots  (a)  and  (b)  are  on  a  natural  logarithmic  scale.  The 

results  are  for  an  'asteroid'  type  object  which  it  shown  in  (c),  together  with  its  Fourier 

modulus. 

Figure  (4)  :  Graph  (a)  compares  the  phase  error  on  corresponding  cross  sections  of  the 
bispectrum  (  Ui«U](,  U2'«U2x  )  and  the  Knox-Thompson  (  Uj»U|,,  Au»Au„ )  sub-planes 
for  ujh^Auk*  4  frequency  sampling  intervals.  Graph  (h)  compares  the  SNR^  of  i'’* 
tame  two  crots-tections  of  the  KT  and  blspectrum  sub-planet. 

Figure  (5)  ;  Plots  of  the  ratio  of  ®b,  defined  by  Eq.(30),  to  l/(V55iVi?„)  for  (a)  the 
blspectrum  plane  defined  by  uj^bAUji,  U|^bu2],b0  ,  and  (b)  the  KT  plane  defined  by 
U]^bO,  AUj,aO.  Departure  of  the  ratio  from  unity  only  occun  near  the  axes. 
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Figure  (6)  :  Diagrammatic  representation  of  systems  of  pupil  sub-ap>ertures,  of  diameter 
showing  how  ‘phase  closure*  is  realised.  The  phase  closure  requirement  can  be 
redefined  in  terms  of  spatial  frequency  vectors  fonrung  a  closed  loop,  (a)  Vectors  con¬ 
tributing  to  Fourier  spatial  frequency  u.  (b)  Phase  closure  is  achieved  in  the  power 
spectrum  by  taking  pairs  of  vectors  of  opposite  sign,  -u  and  -Hi,  to  form  a  closed  loop, 
(c)  Approximate  phase  closure  is  achieved  in  the  Knox  and  Thompson  method  by  tak¬ 
ing  two  vectors,  u  and  u-fiku,  and  assuming  that  the  pupil  phase  is  constant  over  Au. 
The  bispectrum  adds  a  third  vector  Au  to  explicitly  form  phase  closure,  (d)  The  third 
vector,  Au,  is  essential  for  phue  closure  when  >  r^.  The  KT  method  fails  with 
this  arrangement. 

Figure  (7)  ;  The  graphs  illustrate  the  retention  of  phase  information  in  the  KT  and  TC 
techniques.  The  resulu  have  been  obuined  using  a  computer  simulation  of  imaging  a 
point  source  through  a  2  metre  telescope  and  approximately  0.7  are  second  seeing  at 
high  light  level.  The  point  source  is  centroided  and  so  has  xero  Fourier  phase  at  all 
spatial  frequencies.  The  phase  distribution  statistics  arc  obtained  using  an  ensemble  of 
SOOO  different  atmospheric  realisations.  The  fluctuation  of  phase  between  frames  is 
plotted  for  the  first  ISO  realisations  only.  The  telescope  pupil  consists  of  three  widely 
spaced  sub-apertures  of  diameter  Graphs  (a)  show  the  uniformly  distributed  Fourier 
phase  usosciated  with  spatial  fr:equencies  greater  than  r^X.  Graphs  (b)  and  fd)  show 
Fourier  phase  difference  information  is  retained  using  both  the  KT  and  TC  techniques 
for  Au  <  r^K.  Graph  (c)  shows  Fourier  phase  informitlon  is  lost  if  the  KT  is  used 
with  Au  >  r^X.  Graph  (e)  shows  phase  difference  information  is  again  retained,  even 
though  Au  >  rgiK,  due  to  the  phase  closure  property  of  the  bispectrum. 

Figure  (8)  :  Plots  showing  the  moduli  of  corresponding  sub-p'anes  of  (a)  the  btspec- 
trum  (U|BUij,U2«U}y)  and  (b)  the  KT  (U|«U]|,AubAUx  )  four  dimensional  transfer  func¬ 
tions  for  imaging  through  a  2  metre  telescope  with  approximately  0.7  arc  second  ’see¬ 
ing’.  The  data  was  obtained  using  a  Monte  Carlo  computer  simulation,  generating 
10000  independent  realisations  of  the  atmosphere  to  obtain  the  ensemble  average.  Both 
(a)  and  (b)  are  plotted  using  the  same  natural  logarithmic  scale. 

Figure  (9)  :  Graphs  (a)  and  (b)  respectively,  show  corresponding  cross  sections  through 
the  bispectrum  (  Uisui,i,  U2BU2(  )  and  the  Knox-Thompson  (  UjbUix.  Au»Au,  )  sub- 
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planes  of  their  atmosphere-telescope  combination  transfer  functions,  The  cross-sections 
coiretpond  to  frequency  difference  vectors,  U2x>  of  0,  2,  4,  6,  and  8  frequency 
sampling  intervals.  The  telescope  cutoff  occurs  at  50  and  r/K  at  approximately  4.5  fre¬ 
quency  sampling  intervals. 

Figure  (10) :  Plou  of  the  ratio  of  the  phase  error  @e  to  for  the  sub  planes 

of  the  KT  and  bispectrum  of  the  'asteroid'  type  object  shown  in  Fig.3(c)  defined  by 
UiBUix,  U2BU2JI,  AubAU;,.  These  plots  are  the  result  of  computationally  simulating  the 
effect  of  imaging  the  'asteroid'  shown  in  Fig.3(c)  with  a  2  metre  telescope  and  0.7  arc 
second  'seeing'.  Departure  from  unity  only  occurs  near  the  axes. 

Figure  (11)  :  plots  of  a  single  plane  of  the  KT,  (a),  and  bispectrum,  (b), 
(Ui"Ui„  UjBUjp  Au>AU;,  )  SMR's  at  high  light  level  obtained  using  the  computer 
simulation  to  obtain  ensemble  averages  over  10,000  independent  realisations  of  atmos¬ 
pheric  turbulence. 

Figure  (12)  :  Graphs  (a)  and  (b)  respectively,  show  corresponding  cross  sections 
through  the  Knox-Thompson  (  Ui-Uty,  AubAu^  )  and  the  bispectrum  (  Uiauj,,  U}>U}x 
)  sub-planes  of  their  SNR's  shown  in  Pig.ll.  The  cross-sections  correspond  to  fre¬ 
quency  difference  vectors,  U]],,  Au„,  of  0,  2,  4,  6,  and  8  frequency  sampling  Intervals. 
The  telescope  cutoff  occurs  at  50  and  at  approximately  4.5  frequency  sampling 
intervals. 

Figure  (13) :  The  graphs  in  this  figure  are  the  results  of  a  Monte  Carlo  computer  simu¬ 
lation  whereby  ensemble  averages  over  the  atmospheric  sutistlcs  have  been  carried  out 
on  the  photon  noise  dependent  variances  defined  in  Appeiidix(l).  Graph  (l  shows  the 
dependency  of  the  SNR  of  the  bispectrum  on  W,  and  R.  The  relative  SNR  It  the  ratio 
of  the  SNR  for  the  bispectrum  point  defined  by  mi"0.4£)A,,  m2«0  at  r,  =  0.1 85D  to 
that  at  b  0.0950.  This  ratio  of  corresponds  to  a  ratio  of  0.25  in  the  mean  number 
of  speckles  per  frame,  8).  The  plot  tends  to  4  at  very  low  light  level,  corresponding  to 
a  1/8",  depenoency,  see  Eqs.(76),  (77),  and  tends  to  unity  at  high  light  level, 
corresponding  to  a  SNR  which  it  approximately  Independent  of  K,,  tee  Eq.(66).  Graph 
(b)  shows  an  equivalent  plot  to  (a)  for  the  bispectrum  point  defined  by 
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U]b0.60/X,  u^.2DI\,  The  plot  tends  to  a  value  of  16  at  very  low  light  level, 
corretponding  to  a  \IW,^  dependency,  see  Bq.(?5)  and  tends  to  2  at  high  light  level, 
conretponding  to  a  1/V^  dependency,  see  Eq(66).  Graph  (c)  shows  the  approximate 
linear  dependence  upon  N  of  the  SNR  of  the  same  bispectrum  point  when  the  mean 
number  of  photons  per  speckle  is  approximately  one.  The  mean  number  of  speckles 
per  image,  R*,,  in  this  example  is  approximately  1000.  Graph  (d)  shows  the  approxi¬ 
mate  F!^  dependency  of  the  same  bispectrum  point  when  the  mean  number  of  photons 
per  speckle  is  very  much  less  than  unity.  See  Eq.(75).  The  upper  line  is  a  depen¬ 
dent  reference.  Graph  (e)  shows  the  linear  dependence  on  N  at  very  low  light  level  of 
the  SNR's  of  the  KT  and  TC  signals  defined  by  u^aCAD/\,  ugaO  AuaO.  See 
Eqs.(76),(77).  The  upper  line  is  a  linear  reference.  The  KT  signal  has  the  slightly 
higher  SNR,  however  the  l*nes  appear  to  be  superimposed. 

Figure  (14)  ;  Phase  closure  relationships  can  be  viewed  in  terms  of  spatial  frequency 
vectors  forming  closed  loops  in  the  telescope  pupil.  Diagram  (a)  shows  a  representa¬ 
tion  of  the  6  vectors  which  need  to  be  considered  when  analysing  the  modulus  squared 
of  the  bispectrum  at  a  point  defined  by  U|,  u^.  Diagram  (b)  shows  how  the  vectors 
can  be  selected  in  cot\jugatc  pairs  to  cancel  the  phase.  Each  pair  being  analogous  to 
the  power  spectrum  at  that  particular  frequency.  Diagrams  (c)  and  (d)  show  other  pos¬ 
sible  combinations  of  vectors  which  result  in  phase  closure.  Diagrams  (e)  and  (f)  show 
redundent  configurations  between  those  represented  in  (b)-(d). 
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ABSTRACT 

th«  triple  conalation  taehnique  of  tpaekle  Imaginf  hu  boon  Inves¬ 
tigated  in  the  lovr  light  level  regime.  Three  practical  methods  for  digital 
triple  eonelation  processing  have  been  developed.  Bach  of  these 
methods  hu  successfully  reconstructed  imagu  hrom  computer  simulated 
speckle  interferognms.  In  this  paper  wo  describe  and  compare  the  three 
algorithms. 
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1.  Introduction. 

The  angular  reiolution  of  large  ground  bated  teletcopei  li  leverely  limited  by 
itmotpheric  turbulence.  In  order  to  overcome  thii  limiution,  lonte  form  of  ipeckle 
inteiferometry  teohnique  may  be  uied.  Theie  teehniquet  rely  on  taking  a  leriei  of 
meaiurementt.  with  luf&eiently  ihort  expolure  timet  to  freeee  the  effbctt  of  atmoi* 
pherlc  turbulence.  The  retulting  dau  it  then  procetted  to  obtain  the  required 
diffttcdon  limited  information. 

Uiing  the  traditional  interferometiy  technique,  m  tuggetted  by  Labeyrie’  ^  we 
may  obtain  modulut  information  on  the  object  of  intereit  However,  with  the  advent 
of  reatonably  cheap  digital  computing,  it  it  now  becoming  potilble  to  uie  more  com¬ 
plex  methodt  tuch  at  the  triple  correlation^  ^  ^  (in  the  following  referred  to  at  TC)  ind 
Knox-Thompion^  ^  technique!.  With  thete  technique!  we  can  obtain  the  phaee  at  well 
u  the  modului,  and  thut  reconitruct  an  image  of  the  object 


2.  Triple  correlation  procettlng 

The  triple  correlation  /^^(X|,X2)  and  itt  Fourier  traniform,  the  blipecirum  /^^(u.v) 
are  defined  at  followi; 


J  i*(x)i(X'»-Xj)i(x-fX2)dk . 


(1) 


/^\u,v)  ■  /(u)/(v),'*(u+v) .  (2) 

where  i(x)  repteientt  the  tlgnal  and  /(u)  itt  Fourier  trantform. 
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The  buio  methodoloiiei  of  TC  ind  jutoconelation  procttilns  in  very  limilir.  In 
both  cues  we  lum  the  conelition  fUnctiont  of  eieh  fnme  over  the  complete  dau  net. 
From  Eq.  (1)  It  it  appinnt  that  the  TC  of  a  two  dlmeniional  lignal  li  a  foitr  dlmen- 
ilonal  fbncdon  which  leadi  to  difficuldei  in  oonq^uudon  and  manipulation  of  the  triple 
correlation.  In  praodee  we  thenfon  (voceii  only  a  imall  lub^wt  of  the  TC  or  bltpec* 
trum. 

In  the  following  we  will  addreai  ipeciflcally  the  cue  of  photon  limited  interfero* 
grama,  in  pardeutar  the  regime  of  leti  than  one  photon,  on  average,  per  ipeckle.  In  thit 
ngime  data  ii  moit  conveniently  rapnaented  ai  a  teriei  of  time  tagged  photon  coordi¬ 
nate!.  Bach  photo-event  it  thui  npniented  by  a  delu  fUncdon  in  x,y,t  ipace,  le 

8(  x-x*)8(  <“'*)  • 

Integradng  over  a  imall  time  interval  we  obtain  a  dme  frame  containing,  lay,  N  pho- 
toni  which  ii  repreiented  ai^  ; 

ImS 

I(x)  -  1«(  *-*a) .  (3) 

It!  Fourier  traniform  ii, 

/(u)- s' 8(x-x*)  #■*”**,  (4) 

M) 

where  the  product  UX;^  repreienU  the  dot  product  of  the  two  vecton,  (  u  x^^ ).  In  Eqi. 
(3)  and  (4)  X|  it  the  locadon  of  the  kf*  photo-event. 

Most  current  imaging  photon  counting  detectort  are  "dme  framed  devicei"  in 
which  typically  the  light  lignal  it  pre-amplifled  by  an  image  inteniifler,  then  integrated 
for  a  "dme  frame"  of  leveral  tent  of  millitecondt  on  a  film  or  C.C.O.  Due  to  the  lim¬ 
ited  frame  rate  Achievable  by  luch  detecton,  temporal  reioludon  it  it  reladvely  poor. 
To  gain  maximum  SNR  (lignal  to  noiie)  the  frame  rate  muit  be  acyuited,  at  recording 
time,  to  a  rate  dependent  upon  the  coherence  dme  of  the  atmoipheric  turbulence. 

Future  generadoni  of  imaging  photon  counteri .  for  example  the  IPD  *  and 
PAPA  *  detecton,  will  be  of  the  dme  tagging  type.  With  theie  devicei,  individual 
photo-evenu  are  recorded  with  a  temporal  reioludon  typically  of  the  order  of  a  few 
micro-iecondi.  Thii  gain  in  temportU  reioludon  allowi  ui  to  make  better  uie  of  the 
limited  number  of  photo-evenu  from  a  faint  lource.  For  initance,  the  procetiing 
flame  rate  could  be  dynamically  varied  when  proceiilng  the  data,  in  an  attempt  to 
obtain  the  optimum  trade-off  between  ipeckle  contrut  and  lignal  Intenilty. 
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Probably  tha  n.^it  Important  advantage  of  time  framed  data,  i«  that  Instead 
of  dividing  the  dam  stream  into  time  frames,  we  can  look  at  it  through  a  moving  time 
window.  One  finds  that  when  using  the  photon  differencing  algorithm  in  a  moving 
time  window,  it  li  potiible  to  compute  3  times  as  many  diflbrencei  or  blvectori,  than 
it  possible  when  the  data  it  processed  in  frames,  tee  Appendix  and  Pig.  1.  The  SNR 
of  the  bispeetrum  along  its  axea,  at  the  tow  photon  levels  of  interest,  scales  roughly  at 
the  mean  number  of  photons  per  flrame  Equivalently  we  could  say  that  the  SNR 
scales  u  where  is  the  tnean  number  of  bivectors  calculable,  per  frame.  Thus 
increasing  the  number  of  biveetors  calculated  from  a  data>ttream  by  a  factor  of  three, 
will  increase  the  SNR  by  tha  order  of  3‘^  or  1.44. 

2.1  Reconstruction  from  the  bispeetrum. 

The  reconstruction  algorithms  used  in  the  following  sections  are  based  on  a  recur¬ 
sive  technique  .  The  phue  of  the  bispeetrum  lepresentadon  given  in  Eq.  (2)  it: 

3(i)(u,v)  ■  3(u)  +  P(v)  -  3(u+v),  (3) 

where  it  the  phue  of  the  bispeetrum  and  3('")  the  phase  of  the  signal 

transform. 

In  the  cue  of  a  one  dimensional  signal  if  we  put  v  ■  1  we  get  the  following 
equations: 


foru"l 

P(2)-3(i)f  P(i)-3"^(i.i). 

(6) 

foru«2 

P(3)-p(l)-t-P(2)-3W(1.2), 

(7) 

for  14*3 

3(4)  ■  3(1)  +  P(3)  -  3^^(1.3)  • 

(8) 

for  umn 

3(rH-l)  -  3(1)  -»■  3(«)  -  P^’ld.n) . 

(9) 

We  are  flnee  to  choose  the  first  two  phase  values,  since  the  phase  at  zero  fre¬ 
quency  is  irrelevant,  and  the  first  frequency  phue  merely  determines  the  position  of 
the  object  Arbitrarily  setting  the  phue  at  the  fint  frequency  3(1)  to  uro,  we  are  thus 
able  to  compute  the  phue  at  3(2),  using  Eq.  (6) .  Similarly  Eqs.  (7)  to  (9)  allow  us  to 
compute  all  lemainlng  phaw  valuu.  This  reconstruedon  technique  loavea  much  scope 
for  improvement,  since  due  to  its  recursive  nature  it  is  fairly  sensldve  to  noiu.  Unfor¬ 
tunately  a  leut  squares  solution  suffers  from  the  fact  that  (he  bispeetrum  phue  is 
determined  only  modulo  2x  .  When  noise  ia  present  this  can  lead  to  2n  mismatches  in 
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th<  multlpls  phtH  oitimatM  for  each  ftequency.  We  avoid  this  problem  In  the  recur¬ 
sive  algorithm  by  multiplying  and  dividing  complex  exponendali,  rather  than  adding 
and  subtracting  phases. 

3.  Algorithms. 

3.1  Direct  calculation  of  the  Bispectrum. 

With  this  method  we  ealoulate  direetly.  from  the  Fourier  transforms  of  each 
frame,  selected  points  in  the  biapeotrum.  To  do  this  we  make  use  of  the  triple  product 
reladonahip  repnaented  by  Bq.  (2). 

The  processing  steps  are: 

a:  Using  traditional  speckle  Interferometry  techniques,  estimate  the  modulus  of  the 
object 

b:  Decide  which  bispectrum  values  ate  to  be  calculated.  For  maximum  flexibility 
the  coordinates  of  the  bispectrum  points  to  be  calculated  are  formed  into  a  list 
Cl  Calculate  the  Fourier  transform  of  each  frame.  From  the  Fourier  transforms  calcu¬ 
late  the  mean  value  of  each  of  the  listed  bispectrum  pointa,  over  the  data  set. 
d;  Correct  the  bispectrum  for  the  effecu  of  photon  biu 

e:  Using  the  list  of  biapectrum  pointa  as  a  guide,  reconstruct  an  estimate  of  the 

object  transform. 

Since  we  are  only  calculating  a  lirrdted  subset  of  the  total  bispectrum,  we  obvi¬ 
ously  wish  to  select  those  parts  with  the  maximum  SNR.  Even  given  a  knowledge  of 
the  object  spectrum,  calculation  of  the  bispectrum  SNR,  allowing  for  atmospheric  and 
Poisson  noise  is  difflcult^^.  In  general  however  we  can  say  that  given  biipecmim  terms 
of  the  following  form; 

/<*>(u,dv)  ■  /(u)/(dv)/*(u+dv) . 

then  terms  for  which  4v  is  leas  than  Frled't  coherence  parameter  rg/X  ,  will  have  the 
highest  SNR. 

With  photon  limited  data,  due  to  poor  signal  to  noise,  we  do  not  expect  to  be  able 
to  reconstruct  the  object  transform  right  out  to  the  diffraction  limit  of  the  telescope.  In 
order  to  conserve  computation  resources,  we  estimate  the  likely  frequency  cut-off  of 
the  phase  reconstruction,  and  compute  bispectrum  points  only  out  to  this  limit.  Fig  2 
shows  experimentally  determined  curvet  of  cut-off  frequency  vs  number  of  data 
frames,  at  various  photon  levels,  for  a  2m  telescope  with  rg/X  of  approx  18  cm.  These 
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wen  derived  by  attemptins  to  reconitruct  the  phaie  of  a  point  rource,  from  simulated 
ipeokle  data,  using  the  direct  biipectrum  method.  For  more  complex  objects  the 
curves  can  be  used  at  a  upper  bound  on  the  reconstruction  cut-off  frequency.  By  cut¬ 
off  frequency  we  mean,  the  frequency  the  maximum  frequency  to  which  phase  values 
can  be  calculated  with  reasonable  accuracy,  for  example  within  n/4 

3.2  The  Radon  Transform. 

The  ctntral-tUce  or  projeetion-tlice  theorem  lUtei  that  if  we  form  a  one  dimen¬ 
sional  projection  of  a  tvuo  dimensional  fitnetion,  than  the  Fourier  transform  of  that  pro¬ 
jection  Is  a  central-slice  ( ie  passing  through  the  origin)  of  die  Fourier  transform  of  the 
object.  We  can  make  use  of  this  theorem  to  reduce  the  processing  of  a  two  dimen¬ 
sional  ftinction  to  the  processing  of  a  numb«  of  one  dimensional  functions.  We  And 
that  with  the  reladvaly  simple  objects  which  we  have  studied,  a  fairly  small  number  of 
projection  angles  (eg  36)  is  suffleient  for  a  reasonable  reconstruction. 

The  steps  in  processing  aie;- 

a;  Project  the  two  dimensional  data  onto  a  number  of  one  dimensional  lines  at  vari¬ 
ous  angles.  This  is  done  by  rotating  the  coordinate  system  by  an  angle  a  using 
the  rotation  matrix; 

cos(o)  -t-  sin(a)l 
-  tln(o)  cos(a)  J 

Ignoring  the  y  coordinate  effects  the  projecdon. 

b:  For  each  angle  of  projecdon,  compute  the  average  triple  comladon  or  bispectrum 
over  the  data  set  At  this  step  we  may  use  a  photon  differencing  algorithm  to 
compute  the  triple  correlation  (Fig.  2). 

c;  From  the  mean  blspectrum  at  each  angle,  reconstruct  the  correspondl.ng  central 
slice  of  the  object  Fourier  transform. 

d:  Combine  the  various  central  slices  into  an  estimate  of  the  object  Fourier 
transform,  thus  arrive  at  an  estimate  for  the  object  If  possible  object  domain  con- 
stralnu  such  as  posidvlty,  should  be  included  in  this  step. 

In  praedee  the  object  modulus  is  computed  separately  using  the  standard  speckle 
interferometry  techniques.  The  various  projection  triple  correladons  can  be  computed 
numerically,  possibly  dividing  the  calciiladons  among  several  machines.  Alternatively 
since  we  are  woridng  with  one  dii-tensional  fitnedons  a  hybrid  opto-electronic  proces¬ 
sor  could  be  used  to  compute  the  correlations. 
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It  should  be  noted  that  when  reconstructing  from  the  blspectmm  one  is  free  to 
chose  the  values  of  the  transform  phase  at  the  xsroth  and  first  frequencies.  To  fint 
order  it  is  equivalent  to  being  free  to  specify  the  phue  gradient  at  the  origin.  Since 
the  Fourier  traniform  of  the  object  is  anolyde,  it  follows  that  the  phase  gradient  at  any 
point  must  be  Independent  of  direction,  Thus  provided  we  make  appropriate  choices 
for  the  tot  bvo  frequencies  of  ewh  recmisttuction  ( eg  both  uro ),  we  can  guarantee 
there  will  be  no  phase  erroti  between  the  slices.  In  the  image  domain  this  is  equivalent 
to  making  the  centres  of  gravity  of  each  projection  coincident. 


3.3  Sub  plane  calculation. 

Both  of  the  previous  methods  have  Integrated  the  blspectmm  the  series  of  data- 
frames.  It  it  of  course  possible  to  integrate  the  TC  directly.  At  low  light  levels  this  is 
most  conveniently  done  by  using  a  photon  differencing  algorithm. 

To  tee  how  the  photon  differencing  algoilthm  it  formulated,  substitute  Eq.  (3) 
representing  a  photon  limited  image  into  the  triple  correlation  Bq.  (1): 

m 

/frx  f  V  V 

■  d*  Z  8(x-x*)  Z8(k+X|-*,)  Z8(x+*2-xJ  .  (10) 

J  Ml  hi  iml 

M  m*k 

iml 

The  conditions  M  and  are  necestaty  to  eliminate  the  photon  bias  terms .  The 

bias  terms  are  described  in  detail  by  Wimitzer 
The  Integrand  it  non  zero  iff; 


X|  ■  X+Xi 


X,  •  X|ffX. 


X«  -  X+X,  =e  X^  -  X*+X2  . 

Using  equations  (12)  and  (13),  equation  (10)  can  be  rewritten  as  follows. 

iV  N 

■  Z  Z  Z  8(*1-  (*I  -**))  8(xa- 
M)  t>l  M-1 
Ini  mni 


Equation  (U)  clearly  repretenu  a  photon  differenking  process,  see  Fig  3 
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Unfortunately  due  to  the  very  large  lize  of  the  TC  for  any  reaionably  dimen- 
lioned  image,  it  ii  not  practical  to  implement  the  photon  difTerencing  TC  u  it  itandi. 
For  inatance  a  128  iquare  Image  would  normally  generate  a  256k10*  element  TC, 
which  would  require  about  10b  of  itore,  or  170Mb  if  we  made  ftill  uie  of  the  tym- 
metty  of  the  TC. 

Fortunately,  at  the  coit  of  inereuing  the  computational  burden,  (3  integer  opt.  -* 
8  floating  point  opt.  )  we  can  compute  a  tingle  two  dlmentional  tub-plane  of  the 
bitpecirum,  uiing  a  modified  photon  differencing  technique.  Bttendally  we  do  thit  by 
computing  two  dimentiont,  of  the  four  dimentional  Fourier  trantform,  which 
trantformi  the  TC  to  the  blipectnim,  for  tome  conttant  Fourier  vector  V.  At  already 
pointed  out,  the  regiont  eloteit  to  the  axet  of  the  bltpectrum  in  general  have  the 
higheat  SNR,  If  we  compute  thete  high  SNR.  regiont  by  the  photon  differencing  tub- 
plane  method,  the  algorithm  tumt  out,  to  be  timilar  to  the  photon  differencing  method 
of  generating  the  Knox-Thompton  trantform 

Uting  Eq.  (14}  the  bltpectrum  may  be  written  at: 


I^\u,v) 


M  m 


-J«iu  ti.-Vtly  X| 


X  (*/  -**))  *<*2-  (*«-**))• 

ini  ial  mol 


where  to  aid  clarity,  we  have  omitted  the  condidont  that  eliminate  the  photon 


If  we  now  put  V  equal  to  a  conttant,  tay  v  ■  V,  we  may  carry  out  the  X2  integra¬ 
tion  forthwith. 

Since  X]  it  non  xero  only  at  X2  ■  (x^-Xj^  ;• 

P\u,V)  -  I  «<xi  *'  zz'z  -**))  •  (16) 

i-0  m-0 

Separating  out  the  lummation  over  m  :• 


/<^)(u,V) 


m 

f  *'  [S'  1  fz  f  5{*1-  (»/  -*»)) 

^  l_»nO  j(jUI  M 


*(*"'V**) 


131 


-  8  - 


Excluding  the  photon  noiie  blei  lenns  thii  becomet; 


/<^)(u.V) 


£(l-8.,*)(l-8,„j)a<“*’‘'^'-^ 

1 

(18) 


X 


i  Z  «(»!-(*/ 

Im  Ml 


where  bu  1*  Icronecker  delti. 

The  quantity  iniide  the  lecond  tet  of  large  braeketa  can  be  calculated  using  a 
photon  differencing  technique,  where  each  vector  ii  given  u  complex  weighting 

After  finding  the  mean  value  of  thii  quantity  over  the  data  set,  the 


integration  over  X(  (  Fourier  transform)  it  carried  out,  yielding  the  biipectrum  subplane 
/<^>(u,V) . 

We  could  of  course  put  u  and  v  equal  to  a  constant,  this  would  generate  a  single 
bispectum  value.  It  it  however  more  efficient  to  generate  a  whole  plane  at  a  time,  thus 
replacing  two  complex  multiplication!  per  bivector  per  point  with  a  final  Fourier 
transform  after  the  avenging  it  complete. 

For  a  two  dimensional  signal  the  minimum  requirement  for  reconstructing  an 
image  Is  to  have  calculated  the  two  bispeetrum  tubplanes  V  ■  (0,1)  and  V  -  (1,0)  .  In 
this  minimal  configuration  the  method  it  vety  closely  related  to  the  Knox-Thompson 
method. 


4.  RetulU. 

Examples  of  an  object  reconstructed  by  all  three  techniques  are  shown  in  Figs  4 
CO  7.  Each  of  the  reconstructions  represents  roughly  the  same  amount  of  computing 
time,  taking  about  12  seconds  per  data  frame  on  a  Sun  3/160  workstation  (Motorola 
68020  +  68881  at  16.7  MHz,  C  language).  The  direct  method  clearly  produces  the 
best  reconstruction.  However  unlike  the  other  two  cates,  the  program  implementing 
this  method  makes  use  of  an  array  processor,  which  gives  about  a  uix  fold  increase  in 
floating  point  speed.  The  relatively  poor  performance  of  the  Radon  transform,  can 
probably  be  put  down  to  Insufficient  numbers  of  slices,  and  a  very  poor  recombination 
toutine. 

If  one  is  fonunate  enough  to  be  using  time-tagged  data,  then  the  method  of 
choice  would  clearly  be  a  photon  differencing  method.  The  reason  being  that  this 
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method  enables  full  use  of  the  factor  of  1.44  improvement  in  SNR,  which  comes  as  a 
result  of  time*tagginit,  as  described  in  the  Appendix. 

With  conventional  time-framed  data  the  direct  method,  which  is  applicable  in 
both  low  and  high  light  level  regimes  would  seem  to  be  most  appropriate.  However  of 
all  three  algorithms  the  Radon  transform  method  is  the  simplest  to  code,  and  needs 
minimal  amounts  of  memory,  it  would  probably  be  appropriate  for  a  network  of  small 
computers.  Furthermore  it  can  be  implemented  as  a  series  of  photon  differencing 
processes,  and  thus  benefit  from  the  time-tag  Jng  SNR  improvement. 
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Appendix  A. 


In  the  following  we  ihow  that  the  the  SNR  obtainable  by  calculating  photon 
differencing  corretationi  by  integrating  through  a  moving  box-car  window  is 
sighifloantly  impraved  over  that  obtainable  using  a  set  of  discrete  time  windows. 

The  photon  dif'erencing  spatial  auto-correlation  algorithm  speeides  that  we  form 
a  histogram  of  all  tlie  spatial  vector  co-ordinate  differences  betwetm  ali  permuutions  of 
two  non-ldendcal  photons  firom  a  data  ftime.  The  analogous  procedure  for  the  triple 
correlation  it  to  fbrm  a  histogram  of  all  spadal  bi-vector  (tee  Pig.  3)  differences, 
between  all  permutadons  of  three  non-identical  photons  (tom  a  dau-fhune. 

First  let.  us  consider  a  date-stream  consisting  of  a  set  of  evenu  spaced  regularly  in 
time  u  shown  in  Pig.  8.  The  stream  is  divided  up  into  a4Jacent  time  frames  labelled 
(1-4)  of  length  T,  with  N  photons  per  ftwne. 

If  we  have  calculated  the  auto-correlation  for  frame  1,  how  many  additional 
difference  vectors  does  frame  2  contribute?  Using  the  tradidonal  dme-ftome  approach 
there  an  N(N-l)/2  independent  diffrrence  vectors.  In  the  case  of  triple  correlation 
there  an  N(N-l)(N-2)/d  independent  bi-vectort  and  for  k"*  order  correladon  we  get 


Nl 


or  independent  (k-l)*^  order  difference  vectors. 


(N-k)l  kl 

Using  a  moving  window  approach  we  select  each  photon  in  frame  2  in  turn  and 
difference  it  with  the  photons  in  the  preceding  period  of  length  T.  In  the  case  of  the 
auto-correlation,  w«  would  get  (N-1)  Independent  difference  vectors  per  additional 
photon,  giving  a  total  conttlbudon  of  N(N-1)  independent  difference  vectors  from  tlie 
frame.  For  triple  coneladon  we  get  N(N-l}(N-2)/2  independent  bl-vectors  ind  in  the 


case  of  a  k"*  order  correlation  we  get 


Nl 


Independent  (k-l)'*’  order 


(N^)l  (k-l)l 

difference  vectors.  Thus  for  the  auto-correlation  we  get  twice  as  many  difference  vec¬ 
tors,  for  the  triple  correlsrion  three  timer  as  many  bivectors,  and  for  a  k'*'  order  copt- 
ladon  k  times  as  many  difference  vectors. 

Figure  9  illustnstea  a  strasm  of  photons  obeying  uniform  (ie  constant  rate  )  Poll- 
son  sutisdes  in  dme.  For  a  frame  conttining  N  photons  for  a  k*"  (k  $  N)  order  correla¬ 
tion  w*e  get  "-"1^  Independent  (k-1)'*’  order  difference  vectors,  is  in  the  regularly 
spaced  case  above.  Averaging  this  over  uniform  Poisson  stadsdes  gives  e  mein  of 

Rk  o<k-l) 

- ...r. - * - tN 


~  independent  difference  vectors  per  frame,  or 


kl 
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A  buic  property  of  i  uniform  Foiiton  proeeti  ii  that  itt  imdidcs  ve  independent 
of  when  we  ohooM  to  itirt  meuurini  the  proceii.  Thui  if  we  meuure  the  number  of 
photoni  in  a  time  window  of  length  T,  the  diitrlbution  obierved  will  be  independent  of 
the  potitioni  of  the  windowi,  provided  they  are  uniformly  diitributed  throughout  time. 
In  particular  we  may  chnoie  to  observe  the  proceii  through  windows  whose  start  time 
is  itself  a  Poisson  process. 


Using  the  time  vtdndoVv  apixoaeh  on  a  Poisson  process,  we  move  our  window 
along  the  time  axis  undi  itt  leading  edge  is  coincident  with  a  photon.  The  resulting  set 
of  time  frames  will  be  Poisson  distributed  in  time  and  u  pointed  out  above  this  will 
not  effect  the  observed  photon  stuUtdes.  Thus  the  expeetadon  value  of  the  number  of 
photons  in  each  frame,  not  including  the  photon  on  the  leading  edge  is  simply  R 
With  N  photons  in  a  partieu.'ar  frame,  diffstencing  the  leading  photon  with  the  N  pho- 


(N-lc+l)I  (k-l)l 


difrisience  vectors.  When  we  average  this  quantity  over  a 


Poisson  process  we  get  ■  difAsrenee  vectors  per  obsen-ed  photon.  Thus  eompar' 
ing  this  to  the  mean  number  of  independent  difference  vectors  P"  photon 


given  above,  we  get  the  tame  factor  k  increase  in  the  number  of  difference  vector?  at 
we  did  with  a  regular  event  spacing. 


The  SNR  in  the  powe^tpectIum  at  low  photon  levels  scales  at  Ff  ^nd  similarly 
the  SNR  of  those  portions  of  the  bitpectrum  close  to  the  ixet  (close  to  the  power* 
specttum  terms),  where  the  SNR  it  highest,  alto  tcalo  u  IT  Assuming  an  analo* 
gout  result  holds  for  higher  order  spectra,  the  SNR  of  a  k"*  spectrum  ;.cale8  at 
where  \  it  the  number  of  k'*'  order  difference  vectors  per  observed  photon.  Since 
when  calculating  a  k***  order  corrsladon  using  a  condnuoui  dme  window  gives  us  k 
timet  more  difference  vectors,  we  can  say  that  the  SNK  improvement  in  the  order 
’'.yectrum  it  of  the  order  of  This  improvement  in  SNR  is  for  the  auto- 
correladon  and  ’VI  for  the  triple  correladon. 


Figure  1; 

Uiing  ■  condnuoui  tlnu  window  of  length  T,  we  gain  the  correlttion  term 
tetween  the  photoni  C,p^,  which  would  not  erriie  when  ueing  ■  more  conven¬ 
tional  time  framed  data-atream. 

Figujo  2:, 

Cut-off  frequency  of  TC  raconitruction  vt  No.  of  fhunM  .fqr  tqeonitructlon  of  a 
point  lource  uiing  the  direct  method.  The  diameter  of  the  teleicope  is  2m, 
equivalent  to  SO  hint.  Frieda  parameter  la  wt  to  18cm.  The  cut  off  la  arbitrarily 
placed  at  the  radiui  ‘  wham  the  p:.««e  varianoe  of  the  leconatruetion  becomes 
.  greater  than /4.  Teleicope  cut-off  SO  bins. 

Figure  3: 

Above  figure  la  a  graphical  repreienution  of  the  formation  of  a  bl-vector  from 
three  photo-event  cb-oidinatea.  One  of  the  blvectora  implied  by  eqn  (15)  la 
drawn,  both  of  the  component  vecton  are  ahown  u  arrows  between  two  photo¬ 
events.  The  vector  corresponding  to  8(x(X/~X;())  la  shown  pointing  from  photon 
k  to  photon  1,  similarly  for  the  vector  correiponding  to  8(x(x^-X;^)).  Theie  two 
vectors  would  be  catenated  to  form  a  four  vector  corresponding  to  a  triple  corre¬ 
lation  element. 

Figure  4: 

Oifftaction  limited  object  seen  through  a  2^1  telescope 
Figure  S; 

Reconstruction  using  direct  method  with  20(XX)  fitmes  at  200  photoni  per  frame. 
Figure  ti: 

Reconstruction  using  radon  transform  with  23000  fhunes  at  200  photons  per 
frame,  projections  every  4  degrees. 

Figure  7: 

Reconiuuction  using  sub-plane  method  with  12300  frames  at  200  photons  per 
frame,  4  planes  used 
Figure  8; 

A  regular  temporal  process,  divided  into  discrete  time  frames. 

Figure  9; 

A  temporal  Poisson  process,  divided  into  discrete  time  frames. 
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Dlfftnctlon  limited  object  teen  through  a  2m  telescope.  ! 
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Object 
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RMOiiitructlon  uilng  direct  method. 


20000  frames  at  200  photons  per  frame. 

Reconstructed  Object 


Reconstructed  Phase 


Reconitruction  using  ths  Radon  method. 


2SOOO  frunes  at  200  photons  per  frame»  Projections  every  4  degrees. 
Reconstructed  Object 
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Reconstructed  Phase 
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Reconstruction  using  the  Sub*plane  method. 


12500  frames  at  200  photons  per  frame,  4  planes  used. 
Reconstructed  Object 
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Reconstructed  Phase 


Fig? 


U4 


REF  9 


OBJTCCT  RECONSTRUCTION  FROM  PHOTON-LIMITED  CENTROIDED 
DATA  OF  RANDOMLY  TRANSLATING  IMAGES 

L.C.  de  FVeitae  and  J.C.  Dainty 
Blackett  Laboratory 

Imperial  College  •  London  SWT  2BZ  -UK. 

ABSTRACT 

Centroiding  it  investigated  at  a  limple  and  computationelly  fn  >t  technique  of  object  re¬ 
construction,  at  low  light  level,  from  randomly  translating  images.  A  relationship  between  the 
spectrum  of  the  average  N- photon  centroided  frame  and  the  object  spectrum  is  presented  as 
well  as  a,;  algorithm  for  retrieving  the  phase  in  the  case  of  one-dlmen'lonal  objects.  Computer 
simulated  data  is  used  to  test  the  relationship  and  the  reconstruction  algorithm. 


I  .  INTRODUCTION 


A  number  of  techniquei  li&ve  been  developed  to  retrieve  dlt&action-Umited  Mtronomicftl  im¬ 
ages  by  recording  short-exposure  frames  of  speckle  patterns.  The  SHIFT-AND-AOD  (SAA) 
method^‘~*^  is  one  of  the  proposed  techniques.  It  relics  on  the  proposition  that  each  frame 
consists  of  mar.y  distorted  replicas  of  the  true  image  and  that  an  improved  estimate  of  the 
image  can  be  found  by  superlmpoeing  these  distorted  replicas.  At  very  low  light  level  each  of 
these  distorted  images  may  have  only  a  few  photons  and  no  “bri]t;ht'*  speckle  can  be  chosen 
in  order  to  implement  the  SAA  technique. 

Ceuttoiding  is  investigated  as  a  new  way  of  circumventing  this  limitation  and  for  this 
purpose  a  relationship  between  the  N-photon  averaged  centroided  images  and  the  stationary 
normalised  image  /(r)  as  well  as  a  phase  reconstruction  algorithm  is  prtaented  for  the  special 
case  of  randomly  translating  images  of  one-dim  'nsional  objects. 

2  -  THEORY  OF  CENTROIDINQ  DATA  Hi..'  RANDOMLY  TRANSLATING 

IMAGES. 

To  “freese"  a  randomly  moving  image,  photons,  that  are  all  supposed  to  be  emanated 
from  th  e  same  ranuomly  translating  image,  are  detected  during  a  series  of  short  time  inter- 
vals(ff»tnea).  To  retrieve  the  stationary  image  one  should,  for  each  frame  k,  shift  the  photon 
vectors  by  the  amount  that  the  true  centroid  of  the  Image  is  displaced  in  respect  to  the 
centre  of  the  frame  and  average(add)  over  all  the  frames.  The  true  shift  vector  cs  is  unknown 
and  as  an  estimator,  the  centroid  vector  Rs  ■  of  the  detected  photons  is 

evaluated  and  tlie  photon  vectors  are  shifted  by  this  estimator  before  averaging  many  such 
frames. 

The  relationship  between  the  centroidod  and  the  non  centroided  spectrum  of  the  N- 
photon  data  di,{x)  detected  on  the  k*^  frame  is; 

-  exp{+i2iru-R*)Z>*(u,xi,...,xw)  (1) 

and  the  average  of  DJ(u,Xi,  ...,Xftr)  over  all  poMible  sets  of  detected  N-photon  coordinates 
is  performed  using  the  tvormalised  object  intensity  ae  the  probability  deuslty  distribution. 
Thus  Eq.(l)  leads  to  a  quantity  ^^(u)  9  D°(u)/N  related  to  the  normalized  object  intensity 
spectrom^*^; 

Qn(u)  - /(ull-l/AT))  [/(-u/AT)]''’’  (2) 

Expression  (2)  is  a  translation  invariant  relationship  between  the  normalized  spectrum 
of  the  stationary  image  7(u)  and  the  spectrum  $^(u)  of  the  centroided  averege  of  those 
frames  containing  exactly  N  photons  per  frame.  For  instance  for  «  2  Eq.(2)  reduces  to 

2 


U8 


^](u)  a  p(^)|  power  ipectrvun  of  /(u)  waere  ttie  vbriftblc  u  ii  icaled  by  a 

factor  of  2. 

FbrN-3 

5,(u)-  /(|tt)  (3) 

which  is  an  expreuion  related  to  the  btipectrum^^^  when  it  is  evaluated  at  frequencies  u  -» 

V  -^u  or  when  u  -•  -^u,  v  •>«  both  solutions  being  representations  of  lines  in 
the  blspeetnun  u,  v  plane.  It  can  also  be  seen  that  in  the  limiting  case  of  AT  -*  oo  thai. 
Qu(u)  -*  /(u)  as  one  would  expected. 

3  •  OBJECT  RETRIEVAL  FROM  ONB>DIMBNSIONAL  CENTROIDED  IM¬ 
AGES  •  QN(r). 

In  Eq.(2)  If  one  changes  u/Af  -»  u  and  consider  the  ease  of  real  one  dimensional  Intensity 
(/(-u)  -  /•(«))  one  hu: 

(4) 

FVom  Eq.  (4)  follows  a  recurrence  relation  linking  phases  of  the  object,  at  the  discrete 
frequencies  k  and  [Af  -  IJ*  with  the  phase  6w(A?ifc)  of  the  average  quantity  ^/v(Afife); 

-  ^((Af  -  l]k)  -(N-  im)  (6) 

For  N  m  2  photons/friune  63(2!;)  ■  ^(k)  —  ^(k)  m  0  and  therefore  no  information 
about  the  phase  can  be  retrieved.  But  for  Nb3  photons/ frame  Eq.(5)  becomes: 

e,(3k)  -  ^(2k)  -  2^(k)  (6) 

and  hence  from  phase  at  frequency  k  one  can  reconstruct  the  phase  at  2k  up  to  3k  <  (£/2)-l 
where  C  is  the  actual  number  of  bins  used  to  sample  Qs(r). 

The  method  can  be  better  understood  through  an  example  in  which  a  one  dimensional 
image  is  sampled  at  3?  points  and  therefore  ^jv(Afk)  is  determined  only  at  32  frequency  bins. 
Due  to  the  fact  that  7{—u)  >■  r’Cu)  one  has  to  find  the  phases  only  in  half  the  total  number 
of  bins  and  then  reverse  their  sign.  In  this  particular  example  one  has,  therefore,  to  consider 
only  16  frequency  bitui  (Figure  1). 

For  kwO  the  phase  d(0)  ■  0.  As  QN(Nk)  is  shift-invariant  one  can  always  determine 
/(k)  apart  from  a  constant  phase  shifting  factor.  By  appropriately  choosing  this  factor  one 
can  always  set  the  phase  at  ^(1)  >1  0. 


3 
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In  general,  for  a  frequency  hin  k,  where  k  ie  o  prinae  uuuber,  one  can  only  reconatruct  iti 
phase  by  using  an  average  Q;v(Afib)  such  that  N^k+l  photons  per  frame.  In  this  particular 
example  to  reconstruct  the  phaees  up  to  bin  Ifi  one  needs  averages  of  h'ames  with  N  up 
to  16  photons  pc*  frame.  The  modulus  of  T(k)  it  found  from  the  Fourier  transform  of  the 
auto-correlation  function,  the  power  spectrum  |f^(u)|  • 

Experiments  were  carried  out  using  simulated  one  dimensional  photon  data  emitted  from 
the  object  shown  in  Fig.  2(a)  and  in  the  particular  simulation  experiment  described  in  this 
paper,  we  set  the  Poisson  distribution  mean  value  TT  ••  S  photons/frame.  Frames  v.'ith  0  or 
1  photons/frame  rjre  disregarded,  because  they  do  not  carry  any  information  concerning  the 
intensity  distribution  of  the  image,  and  a  total  of  S0,084  frames  of  a  randomly  translating 
image  (of  a  binary  star)  containing  N  ^  3  photons/frame  were  generated.  Theee  frames  can 
be  grouped  in  sets  of  frames  containing  a  number  of  photons  ranging  from  Nb3  up  to  NwlS 
photons/frame  in  this  case. 

Figure  3  depicts  Qjv(r)  for  Nw2  and  NaS  photons/frame,  hbr  Naa2  it  can  be  seen  that 
Qt(T)  exhibits  the  same  shape  as  the  autocorrelation  of  the  binary  but  distributed  in  a  smaller 
region  of  space. 

Figure  2(b)  shows  the  reconstructed  image  obtained  using  this  method  with  C  ■  128 
bins  and  using  frames  with  N  only  up  to  13  photons/frame  i.e.  not  all  the  frequencies  were 
reconstructed.  Despite  the  fret  that  no  other  image  proceeslng  technique  has  been  used,  e  g. 
enforcing  positivity.  Fig.  2(b),  nevertheless,  gives  a  good  estimate  of  the  image  and  in  thin 
particular  case  also  a  good  estimate  of  the  relative  brightness  of  the  stare. 
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Fig.  1  •  Pictorial  raprctantatlon  of  tho  phat*  raconstruction  algorithm  for  Nag 
and  Na4  photoni/frama 


Fig.  3  •  (a)  •  Stationary  Object  Intanaltyi  (b)  •  Raconatruetton  with  N  ranging 
fkrom  3  up  to  13  photona/firama 


Fig.  3  •  (a)  -  Centroldad  framaa  with  Na3  photona/firama  and  (b)  -  Centrolded 
framaa  with  Nag  photoni/frame 
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